Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 955,154 G→A 20.9% intergenic (+105/‑208) cyoE → / → alaC protoheme IX farnesyltransferase/aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947955,1540GA20.9% 822.8 / 216.6 411intergenic (+105/‑208)cyoE/alaCprotoheme IX farnesyltransferase/aminotransferase
Reads supporting (aligned to +/- strand):  ref base G (157/164);  new base A (37/48);  total (194/217)
Fisher's exact test for biased strand distribution p-value = 3.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.82e-04

Only 100 of 421 total aligned reads displayed.

GACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCAGGCC  >  NC_002947/955005‑955303
                                                                                                                                                        |                                                                                                                                                      
gACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAgg                                                                                                                                                        <  2:1708323/150‑1 (MQ=255)
 aCAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGt                                                                                                                                                       <  2:5020116/150‑1 (MQ=255)
 aCAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGCAAACCCGCTCCTACAAGGt                                                                                                                                                       >  2:2738402/1‑150 (MQ=255)
   aGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTc                                                                                                                                                     <  1:8238716/150‑1 (MQ=255)
     ccAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCt                                                                                                                                                   <  2:8602627/150‑1 (MQ=255)
     ccAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCGACAAGGTTCCt                                                                                                                                                   >  2:5268511/1‑150 (MQ=255)
        gACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTa                                                                                                                                                <  2:8326774/150‑1 (MQ=255)
        gACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTa                                                                                                                                                <  2:8146751/150‑1 (MQ=255)
        gACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTa                                                                                                                                                <  1:6809255/150‑1 (MQ=255)
        gACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTa                                                                                                                                                <  1:2790856/150‑1 (MQ=255)
        gACCGCCGCGGATGTGCTGAGGACCTCCGCCCGCTGACACTGCGGCCTGTGTCACGAAAACCCCGGCCACTGTGCCGGGGATTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTa                                                                                                                                                <  2:5691940/150‑1 (MQ=38)
         aCCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAg                                                                                                                                               <  1:320048/150‑1 (MQ=255)
          ccGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAgg                                                                                                                                              <  2:6625990/150‑1 (MQ=255)
             cTGCGGATGTGCTGATGACCTACGCACGCTGTAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGc                                                                                                                                           <  2:6502237/150‑1 (MQ=38)
             cTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGc                                                                                                                                           <  1:7523382/150‑1 (MQ=38)
              tGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGc                                                                                                                                           <  1:7251919/150‑1 (MQ=37)
              tGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCgg                                                                                                                                         <  2:2791005/150‑1 (MQ=25)
               gCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGTCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTTCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCgg                                                                                                                                         <  1:6380271/150‑1 (MQ=255)
               gCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCgg                                                                                                                                         <  1:4018328/150‑1 (MQ=38)
               gCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCgg                                                                                                                                         <  2:5992223/150‑1 (MQ=38)
               gCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCgg                                                                                                                                         <  1:1815863/150‑1 (MQ=255)
               gCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCgg                                                                                                                                         <  1:6721747/150‑1 (MQ=255)
               gCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGACGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCgg                                                                                                                                         <  2:2407181/150‑1 (MQ=255)
                cGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCggg                                                                                                                                        <  2:3824404/150‑1 (MQ=38)
                cGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCggg                                                                                                                                        <  2:4356105/150‑1 (MQ=38)
                cGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCggg                                                                                                                                        <  2:6381034/150‑1 (MQ=38)
                cGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCggg                                                                                                                                        <  2:3034314/150‑1 (MQ=255)
                cGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCggg                                                                                                                                        <  1:5496723/150‑1 (MQ=255)
                cGGATGTGCTGATGACCGACGCACGCTGAAACTGCGGCCTGTTTCACGAAAACCCCGGCCACGGTGCCGGGGTTTTTTTTTGGGGGTCTTTTGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCGGAACCTTGTAGGAGCggg                                                                                                                                        <  2:8787656/150‑1 (MQ=25)
                cGGAGGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTTGGGGGTATTTTTCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCgg                                                                                                                                         <  2:4030943/150‑1 (MQ=255)
                 ggATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCGGGt                                                                                                                                       <  2:751066/150‑1 (MQ=37)
                 ggATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGt                                                                                                                                       <  1:8484089/150‑1 (MQ=255)
                       gCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTAcc                                                                                                                                  <  1:4953000/150‑1 (MQ=255)
                               ccTACGCACGCTTAAACTGCTGCCTTTTTAACGAAAACCCCGGCCACTGTGCCGGGTTTTTTTTTTGGGGTTTTTTTTCCTGTGCGGGCCTCTTCGCGGGTAAACCCGCTCCTACATGGTTCCTGTAGGTGGGGTTTTACCCGCGAagag                                                                                                                         <  2:3702881/150‑1 (MQ=37)
                                cTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTAAAAGGTTGCTGTAGGAGCGGGTTTAACCCCGAAGAgg                                                                                                                        >  2:4738864/1‑150 (MQ=255)
                                cTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGGAAACCCCCGCCGACAAGGTGGCGGGAGGAGCGGGGTTACCCGCGCAGAgg                                                                                                                        >  2:4565795/1‑150 (MQ=255)
                                 tACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCGCCGACAAGGGTCCGGGAGGAGCGGGTTTACCCGCGAAGAGGc                                                                                                                       >  1:8747581/1‑150 (MQ=38)
                                   cgcacgcTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGGAAAACCGCTCCGACAAGGGTCGGGGAGGGGCGGGGTTAACCGCGGAGAgggcc                                                                                                                     >  1:1306187/1‑147 (MQ=255)
                                        gcTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTTAACCCGCTCCTACAGGAACCTTGTTGGAGCggg                                                                                                                                        <  1:1729125/126‑1 (MQ=25)
                                              cctgcggCCTGTGACACGAAAACCCCGGCCCCGGGGCGGGGTTTTTTTTTTGGGGGTTTTTTACCTGGGCCGGCAACTCCGCGGGTAAACCCGCTCCTACACGGTTCCTGTCGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGaaaaa                                                                                                         <  1:3633301/149‑1 (MQ=21)
                                              aCTGCTGCCTGTGACACGAACACCCCGGCCACTGTGCCGGGGTATTTTTTGGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAACCCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGaaaaa                                                                                                         <  2:5390547/150‑1 (MQ=37)
                                                  ctgcCTGTTTCACGAAACCCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTCCAGGAACATTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGaaaaatacc                                                                                                     <  2:4302060/150‑3 (MQ=14)
                                                     ccTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCGCCAACAAGGTTGCGGTTGGTGCGGGGTTTCCCGCGTAGAGGCCCGCCCAGgcaaataaaaat                                                                                                  >  2:1798857/1‑150 (MQ=255)
                                                       tGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCgg                                                                                                                                         <  1:5073784/110‑1 (MQ=25)
                                                         tttCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTc                                                                                                                                                     <  2:5630411/96‑1 (MQ=255)
                                                                  aaCCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCgg                                                                                                                                         <  1:7609639/99‑1 (MQ=18)
                                                                  aaCCCCGGCACAGTGGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTGAAACCCCGCCAACAAGGTTGCGGGAGGTGCGGGTTTTACCCCCCAGAAGCCAGCACACAAAAAAAAAAATAACAGAATTTaaa                                                                                     >  1:2063236/1‑150 (MQ=255)
                                                                    ccccGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAAACCGCTCCTACAAGGTTGCGGTAGGAGCGGGTTTACCCGCGCAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTGaaaa                                                                                   >  2:5579092/1‑150 (MQ=38)
                                                                      ccGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGGTAAACCCCGACGACAACGGTGCGGTGGGAGCCGGGTTTCCCCCGCAGAGGCCCCCCCAGAAAAAAAAAAATACCAGATTTTAaaaaaa                                                                                 >  1:7132738/1‑150 (MQ=255)
                                                                      ccGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCaaaaaa                                                                                 >  2:5098729/1‑150 (MQ=255)
                                                                        ggCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATa                                                                               >  2:8082273/1‑150 (MQ=255)
                                                                           cACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCCGGTAAAACCGCACCAACAAAGTTAAGGGAGGAGCGGGTTTACCCCCCAACAGGCCAGCAAAGCAAAAAAAAAATACCAGAATTTAaaaaaaaaaaa                                                                            >  1:7928002/1‑145 (MQ=255)
                                                                            aCTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGTTCCGGGAGGAGCGGGGTTACCCGCGAAGAGGCCCGCCCAGAAAAATAAAAATACCTGATTTTGAAAAAATAGATc                                                                           >  2:7270316/1‑150 (MQ=37)
                                                                                       ggTTTTTTTTTGGGGGTATTTTTCCTGTGCGGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGTTACCTGTCGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTc                                                                 <  1:8583000/150‑1 (MQ=255)
                                                                                       ggTTTTTTTTGGGGGGTATTTTTCCGGTGCCGGCCTCGGCGCGGGAAACCCCGCACCTCCACGGTTCCTGTAGGCGCGGTTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTc                                                                 <  1:577542/150‑1 (MQ=255)
                                                                                           ttttttttGGGGTTATTTTTCCTGTGCCGGCCTCTTCGCGGTTAAACCCGCTCCAACAAGGTTCCTGTCGGAGCGGGTTAACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCtaa                                                              <  2:2205830/150‑1 (MQ=38)
                                                                                           ttttttcgGGGGGTATTTGGCCTGCGCCGGCCTCTTCGCGGGTAACCCCGCCCCTACCAGGTTCCTGTCGGAGCGGGTTTCCCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCtaa                                                              <  2:4802189/150‑1 (MQ=37)
                                                                                                    gggtTATTTTTCCTGTGCGGGCCTCTTGGCGGGAAAACCCGCTCCTACATGGTTCGTGTTGGTGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAgg                                                     <  1:4257158/150‑1 (MQ=37)
                                                                                                      ggTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTATAATAAAACAGGaa                                                   >  2:5846366/1‑150 (MQ=255)
                                                                                                       gTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGaaa                                                  <  2:8565984/150‑1 (MQ=255)
                                                                                                              gccTGTGCCGGCCTCTTCGCGGTTAAACCCGCTCCTCCATGGTGCCTGTAGGAGCGGTTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCt                                           <  1:603710/149‑1 (MQ=38)
                                                                                                               ccTGGGCCGGCCTCTTCGCGGGTAAACCGGCTCCTACACGGTTCCTGTAGGTGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCtt                                          <  2:2559593/150‑1 (MQ=255)
                                                                                                                 ggtgCCGGCCCCTCCGCGGGTAAACCCGCCCCAACACGGTGCGTTTAGGAGGGGTTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTa                                        <  1:6184785/149‑1 (MQ=255)
                                                                                                                        gCCTCTTCGCGGGTAAACCCGCTCCTACCCGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCAt                                 <  2:6981406/150‑1 (MQ=255)
                                                                                                                               cgcgGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGttt                          <  2:7936580/150‑1 (MQ=255)
                                                                                                                                 cgGGAAAACCCGCTCCTACAAGGTTCCGGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCg                        <  2:7595259/150‑1 (MQ=255)
                                                                                                                                   gggAAAACCCGCCCCTACAAGGTTCCTTTAGTAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCgg                       <  1:5572759/150‑1 (MQ=255)
                                                                                                                                     gTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAACTTATCCTTTACAGGCATCCTGTTTCGGCa                     >  2:6779707/1‑150 (MQ=255)
                                                                                                                                      tAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCAc                    >  1:5622385/1‑150 (MQ=255)
                                                                                                                                          cccGCTCCTCCCCGGTCCCGTTCGGCGCGGTTTAACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtc                <  1:4565795/150‑1 (MQ=255)
                                                                                                                                          cccGCTCCTACAGGAACCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtc                >  2:1997432/1‑150 (MQ=38)
                                                                                                                                          cccGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtc                >  1:3709299/1‑150 (MQ=255)
                                                                                                                                          cccGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtc                >  1:2629324/1‑150 (MQ=255)
                                                                                                                                           ccGCTCCTACAGGAACCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtct               >  2:3583686/1‑150 (MQ=38)
                                                                                                                                           ccGCTCCTACAGGAACCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtct               >  1:971012/1‑150 (MQ=38)
                                                                                                                                           ccGCTCCTACAGGAACCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtct               >  2:122773/1‑150 (MQ=38)
                                                                                                                                           ccGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtct               >  1:1001604/1‑150 (MQ=255)
                                                                                                                                           ccGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtct               >  1:6608016/1‑150 (MQ=255)
                                                                                                                                           ccGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtct               >  1:7028271/1‑150 (MQ=255)
                                                                                                                                           ccGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtct               >  2:3123788/1‑150 (MQ=255)
                                                                                                                                           ccGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTGCTAATAAAACAGGAAATTATCCTTAACAGTCATCCTGTTTCGGCACTAtct               >  1:881742/1‑150 (MQ=255)
                                                                                                                                             gCTCCTACAGGAACCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtcttc             >  2:1647905/1‑150 (MQ=37)
                                                                                                                                             gCTCCTACAGGAACCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtcttc             >  1:2915543/1‑150 (MQ=37)
                                                                                                                                             gCTCCTACAGGAACCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtcttc             >  2:4105862/1‑150 (MQ=38)
                                                                                                                                             gCTCCTACAGGAACCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtcttc             >  2:7026870/1‑150 (MQ=38)
                                                                                                                                             gCTCCTACACGGTTCCTGTTGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAGAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtcttc             <  2:1099974/150‑1 (MQ=255)
                                                                                                                                             gCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAACCAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtcttc             <  2:8747581/150‑1 (MQ=255)
                                                                                                                                             gCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtcttc             >  1:7760951/1‑150 (MQ=255)
                                                                                                                                               tCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTg           >  1:5329623/1‑150 (MQ=255)
                                                                                                                                               tCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTg           >  1:994714/1‑150 (MQ=255)
                                                                                                                                                 ctccccggTTCCTGTCGGCGCGGGTTTCCCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGAt         <  1:4138247/144‑1 (MQ=255)
                                                                                                                                                  tACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGAtt        >  2:1413931/1‑150 (MQ=255)
                                                                                                                                                  tACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGAtt        >  2:6327242/1‑150 (MQ=255)
                                                                                                                                                  tACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTACTGAtt        >  1:7372062/1‑150 (MQ=255)
                                                                                                                                                   aCAAGGTTCCTGGCGGAGCGGTTTCCCCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTc       <  1:1225808/150‑1 (MQ=255)
                                                                                                                                                    attggTTCCTGTTGTTGCGGTTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCa      <  1:385033/147‑1 (MQ=255)
                                                                                                                                                     aaGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCAg     >  2:7821104/1‑150 (MQ=255)
                                                                                                                                                     aaGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCAg     >  1:242069/1‑150 (MQ=255)
                                                                                                                                                     aaGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCAg     >  2:3730959/1‑150 (MQ=255)
                                                                                                                                                        gTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCAGgcc  >  2:3298086/1‑150 (MQ=255)
                                                                                                                                                        |                                                                                                                                                      
GACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCAGGCC  >  NC_002947/955005‑955303

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