Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_007793 1,693,048 (T)7→8 coding (1124/1125 nt) SAUSA300_RS08410 ← coproporphyrinogen III oxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0077931,693,0411.T97.3% 102.4 / ‑3.5 37intergenic (‑94/+6)SAUSA300_RS08405/SAUSA300_RS08410HrcA family transcriptional regulator/coproporphyrinogen III oxidase
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base T (1/35);  total (1/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.34e-01

CACTCACTTATCTCAAGTGCTAATTATAATTTATCAAATTGGTCAAAGTAAGTCAATGTTAAAGACTCGAAATTTCAA‑‑‑‑‑TTTT‑TTTAATCATTTATTAGGAAAGCTTCAAAAACCTCATTACCTATGACTTTCCCTCTATTTGTAAGTGCAATCACATCGTTCTTTT  >  NC_007793/1692964‑1693129
                                                                              |                                                                                             
cactcactTATCTCAAGTGCTAATTATAATTTATCAAATTGGTCAAAGTAAGTCAATGTTAAAGACTCGAAATTTCAAT‑‑‑‑TTTT‑TTTAATCATTTATTAgg                                                                     <  1:5255856/100‑1 (MQ=255)
           ctcAAGTGCTAATTATAATTTATCAAATTGGTCAAAGTAAGTCAATGTTAAAGACTCGAAATTTCAAT‑‑‑‑TTTT‑TTTAATCATTTATTAGGAAAGCTTCaag                                                          <  1:2947108/100‑2 (MQ=255)
             cAAGTGCTAATTATAATTTATCAAATTGGTCAAAGTAAGTCAATGTTAAAGACTCGAAATTTCAATTTTTTTTT‑TTTAATCATTTATTAGGAAAGCTTCa                                                            >  1:6595374/1‑100 (MQ=255)
                    tAATTATAATTTATCAAATTGGTCAAAGTAAGTCAATGTTAAAGACTCGAAATTTCAAT‑‑‑‑TTTT‑TTTAATCATTTATTAGGAAGGCTTCAAAAACCTCAtt                                                 <  1:5213420/100‑1 (MQ=255)
                        tataATTTATCAAATTGGTCAAAGTAAGTCAATGTTAAAGACTCGAAATTTCAAT‑‑‑‑TTTT‑TTTAATCATTTATTAGGAAAGCTTCAAAACCCTCATTACCt                                             <  1:3077356/100‑1 (MQ=255)
                          taATTTATCAAATTGGTCAAAGTAAGTCAATGTTAAAGACTCGAAATTTCAAT‑‑‑‑TTTT‑TTTAATCATTTATTAGGAAAGCTTCAAAAACCTCATTACCTAt                                           <  1:1253410/100‑1 (MQ=255)
                           aaTTTATCAAATTGGTCAAAGTAAGTCAATGTTAAAGACTCGAAATTTCAAT‑‑‑‑TTTT‑TTTAATCATTTATTAGGAAAGCTTCAAAAACCTCATTACCTagg                                          <  1:1958975/100‑3 (MQ=255)
                           aaTTTATCAAATTGGTCAAAGTAAGTCAATGTTAAAGACTCGAAATTTCAA‑‑‑‑‑TTTTATTTAATCATTTATTAGGAAAGCTTCAAAAACCTCATTACCTATg                                          <  1:6386888/100‑1 (MQ=255)
                                             aagtaagtCAATGTTAAAGCCTCGAAATTTCAAT‑‑‑‑TTTT‑TTTAATCATTTATTAGGAAAGCTTCACAACCCTCATTACCTATGACTTTCCCTCTATTTGTa                        <  1:5046961/100‑1 (MQ=255)
                                             aagtaagtCAATGTTAAAGACTCGAAATTTCAATT‑‑‑TTTT‑TTTAATCATTTATTAGGAAAGCTTCAAAAACCTCATTACCTATGACTTTCCCTCTacggca                         <  1:3116038/100‑6 (MQ=255)
                                             aagtaagtCAATGTTAAAGACTCGAAATTTCAAT‑‑‑‑TTTT‑TTTAATCATTTATTAGGAAAGCTTCAAAAACCTCATTACCTATGACTTTCCCTCTATTTGTa                        <  1:292885/100‑1 (MQ=255)
                                             aagtaagtCAATGTTAAAGACTCGAAATTTCAAT‑‑‑‑TTTT‑TTTAATCATTTATTAGGAAAGCTTCAAAAACCTCATTACCTATGACTTTCCCTCTATTTGTa                        <  1:3217249/100‑1 (MQ=255)
                                             aagtaagtCAATGTTAAAGACTCGAAATTTCAAT‑‑‑‑TTTT‑TTTAATCATTTATTAGGAAAGCTTCAAAAACCTCATTACCTATGACTTTCCCTCTATTTGTa                        <  1:3319835/100‑1 (MQ=255)
                                             aagtaagtCAATGTTAAAGACTCGAAATTTCAAT‑‑‑‑TTTT‑TTTAATCATTTATTAGGAAAGCTTCAAAAACCTCATTACCTATGACTTTCCCTCTATTTGTa                        <  1:589915/100‑1 (MQ=255)
                                             aagtaagtCAATGTTAAAGACTCGAAATTTCAAT‑‑‑‑TTTT‑TTTAATCATTTATTAGGAAAGCTTCAAAAACCTCATTACCTATGACTTTCCCTCTATTTGTa                        <  1:1884689/100‑1 (MQ=255)
                                               gtaagtCAATGTTAAAGACTCGAAATTTCAAT‑‑‑‑TTTT‑TTTAATCATTTATTAGGAAAGCTTCAAAAACCTCATTACCTATGACTTTCCCTCTATTTGTAAg                      <  1:2072573/100‑1 (MQ=255)
                                                  agtCAATGTTAAAGACTCGAAATTTCAAT‑‑‑‑TTTT‑TTTAATCATTTATTAGGAAAGCTTCAAAAACCTCATTACCTATGCCTTTCCCTCTATTTGTAAGTGc                   <  1:3861186/100‑1 (MQ=255)
                                                  agtCAATGTTAAAGACTCGAAATTTCAAT‑‑‑‑TTTT‑TTTAATCATTTATTAGGAAAGCTTCAAAAACCTCATTACCTATGACTTTCCCTCTATTTGTAAGTGc                   <  1:3382707/100‑1 (MQ=255)
                                                  agtCAATGTTAAAGACTCGAAATTTCAAT‑‑‑‑TTTT‑TTTAATCATTTATTAGGAAAGCTTCAAAAACCTCATTACCTATGACTTTCCCTCTATTTGTAAGTGc                   <  1:6290297/100‑1 (MQ=255)
                                                  agtCAATGTTAAAGACTCGAAATTTCAAT‑‑‑‑TTTT‑TTTAATCATTTATTAGGAAAGCTTCAAAAACCTCATTACCTATGACTTTCCCTCTATTTGTAAGTGc                   <  1:1741956/100‑1 (MQ=255)
                                                  agtCAATGTTAAAGACTCGAAATTTCAAT‑‑‑‑TTTT‑TTTAATCATTTATTAGGAAAGCATTAAAAACCTCATTACCTATGACTTTCCCTCTATTTGTAAGTGc                   <  1:6164266/100‑1 (MQ=255)
                                                   gtCAATGTTAAAGACTCGAAATTTCAATT‑‑‑TTTT‑TTTAATCATTTATTAGGAAAGCTTCAAAAACCTCATTACCTATGACTTTCCCTCTATTTGTAAGTGc                   <  1:4147288/100‑1 (MQ=255)
                                                   gtCAATGTTAAAGACTCGAAATTTCAAT‑‑‑‑TTTT‑TTTAATCATTTATTAGGAAAGCTTCAAAAACCTCATTACCTATGACTTTCCCTCTATTTGTAAGTGCa                  <  1:1260743/100‑1 (MQ=255)
                                                   gtCAATGTTAAAGACTCGAAATTTCAAT‑‑‑‑TTTT‑TTTAATCATTTATTAGGAAAGCTTCAAAAACCTCATTACCTATGACTTTCCCTCTATTTGTAAGTGCa                  <  1:4487099/100‑1 (MQ=255)
                                                   gtCAATGTTAAAGACTCGAAATTTCAAT‑‑‑‑TTTT‑TTTAATCATTTATTAGGAAAGCTTCAAAAACCTCATTACCTATGACTTTCCCTCTATTTGTAAGTGCa                  <  1:4373000/100‑1 (MQ=255)
                                                   gtCAATGTTAAAGACTCGAAATTTCAAT‑‑‑‑TTTT‑TTTAATCATTTATTAGGAAAGCTTCAAAAACCTCATTACCTATGACTTTCCCTCTATTTGTAAGTGCa                  <  1:2289182/100‑1 (MQ=255)
                                                   gtCAATGTTAAAGACTCGAAATTTCAAT‑‑‑‑TTTT‑TTTAATCATTTATTAGGAAAGCTTCAAAAACCTCATTACCTATGACTTTCCCTCTATTTGTAAGTGCa                  <  1:1715595/100‑1 (MQ=255)
                                                    tCAATGTTAAAGACTCGAAATTTCAAT‑‑‑‑TTTT‑TTTAATCATTTATTAGGAAAGCTTCAAAAACCTCATTCCCTATGACTTTCCCTCTATTTGTAAGTGCaa                 <  1:5792425/100‑1 (MQ=255)
                                                    tCAATGTTAAAGACTCGAAATTTCAAT‑‑‑‑TTTT‑TTTAATCATTTATTAGGAAAGCTTCAAAAACCTCATTACCTATGCCTTTCCCTCTATTTGTAAGTGCaa                 <  1:5598303/100‑1 (MQ=255)
                                                    tCAATGTTAAAGACTCGAAATTTCAAT‑‑‑‑TTTT‑TTTAATCATTTATTAGGAAAGCTTCAAAAACCTCATTACCTATGACTTTCCCTCTATTTGTAAGTGCag                 <  1:4041198/100‑2 (MQ=255)
                                                     cttccgatctgACTCGAAATTTCAATTTT‑TTTT‑TTTAATCATTTATTAGGAAAGCTTCAAAAACCTCATTACCTATGACTTTCCCTCTATTTGTAAGTGc                   <  1:6431343/90‑1 (MQ=255)
                                                     cAATGGTAAAGACTCGAAATTTCAAT‑‑‑‑TTTT‑TTTAATCATTTATTAGGAAAGCTTCAAAAACCTCATTACCTATGACTTTCCCTCGATTTGTAAGTGCagg                <  1:4288336/100‑3 (MQ=255)
                                                      aaTGTTAAAGACTCGAAATTTCAAT‑‑‑‑TTTT‑TTTAATCATTTATTAGGAAAGCTTCAAAAACCTCATTACCTATGACTTTCCCTCTATTTGTAAGTGCaagc               <  1:6216310/100‑3 (MQ=255)
                                                      aaTGTTAAAGACTCGAAATTTCAAT‑‑‑‑TTTT‑TTTAATCATTTATTAGGAAAGCTTCAAAAACCTCATTACCTATGACTTTCCCTCTATTTGTAAGTGCAATc               <  1:5224212/100‑1 (MQ=255)
                                                             aaGACTCGAAATTTCAAT‑‑‑‑TTTT‑TTTAATCATTTATTAGGAAAGCTTCAAAAACCTCATTACCTATGACTTTCCCTCTATTTGTAAGTGCAATCACATCGt        <  1:5863123/100‑1 (MQ=255)
                                                                  tagAAATTTCAAT‑‑‑‑TTTT‑TTTAATCATTTATTAGGAAAGCTTCAAAAACCTCATTACCTATGACTTTCCCTCTATTTGTAAGTGCAATCACATCGTTCttt   <  1:4047827/98‑1 (MQ=255)
                                                                   cGAAATTTCAAT‑‑‑‑TTTT‑TTTAATCATTTATTAGGAAAGCTTCAAAAACCTCATTACCTATGACTTTCCCTCTATTTGTAAGTGCAATCACATCggaagagc  <  1:779802/100‑8 (MQ=255)
                                                                              |                                                                                             
CACTCACTTATCTCAAGTGCTAATTATAATTTATCAAATTGGTCAAAGTAAGTCAATGTTAAAGACTCGAAATTTCAA‑‑‑‑‑TTTT‑TTTAATCATTTATTAGGAAAGCTTCAAAAACCTCATTACCTATGACTTTCCCTCTATTTGTAAGTGCAATCACATCGTTCTTTT  >  NC_007793/1692964‑1693129

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑