Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 955,154 G→A 34.5% 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,1540GA34.5% 135.9 / 122.2 133intergenic (+105/‑208)cyoE/alaCprotoheme IX farnesyltransferase/aminotransferase
Reads supporting (aligned to +/- strand):  ref base G (45/40);  new base A (28/17);  total (73/60)
Fisher's exact test for biased strand distribution p-value = 3.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.65e-02

Only 100 of 136 total aligned reads displayed.

GACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCAGGCC  >  NC_002947/955005‑955303
                                                                                                                                                      |                                                                                                                                                      
gACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAgg                                                                                                                                                        <  1:360339/150‑1 (MQ=255)
 aCAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTACTTTTCCTGTGCCGGCCTCTTCGCGGGTAAAACCGCTCCAACAAGGt                                                                                                                                                       >  2:5027549/1‑150 (MQ=255)
 aCAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTGTTTGGGGGTATTGTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCGCCTACCAGGt                                                                                                                                                       >  2:3931561/1‑150 (MQ=255)
   aGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGTGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTc                                                                                                                                                     <  2:5305193/150‑1 (MQ=255)
   aGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTc                                                                                                                                                     <  2:5152556/150‑1 (MQ=255)
    gCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTcc                                                                                                                                                    <  1:113328/150‑1 (MQ=255)
    gCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTcc                                                                                                                                                    <  1:1700784/150‑1 (MQ=255)
     ccAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCt                                                                                                                                                   <  2:2352210/150‑1 (MQ=255)
     ccAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCt                                                                                                                                                   <  2:2352124/150‑1 (MQ=255)
     ccAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAAACCGCGCCAACAAGGttacg                                                                                                                                                   >  1:1782552/1‑147 (MQ=255)
     ccAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGGAAACCCGCTCCAACAAGGGTccg                                                                                                                                                   >  1:5425632/1‑149 (MQ=255)
        gACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTgga                                                                                                                                                >  2:3364682/1‑148 (MQ=255)
         aCCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAg                                                                                                                                               <  2:2140003/150‑1 (MQ=255)
           cGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGa                                                                                                                                             >  2:3764361/1‑150 (MQ=255)
            gCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAg                                                                                                                                            <  1:3764361/150‑1 (MQ=255)
             cTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGc                                                                                                                                           <  1:3220338/150‑1 (MQ=38)
               gCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCgg                                                                                                                                         <  1:4379376/150‑1 (MQ=38)
               gCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCgg                                                                                                                                         <  1:4415439/150‑1 (MQ=38)
               gCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCgg                                                                                                                                         <  1:1183493/150‑1 (MQ=38)
               gCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCgg                                                                                                                                         <  1:4494051/150‑1 (MQ=255)
                cGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCggg                                                                                                                                        <  2:1046327/150‑1 (MQ=38)
                cGGAGGGGCTGATGACCTCCGCACGCTGAAACTGCTGCCTGTGTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCggg                                                                                                                                        <  2:3786949/150‑1 (MQ=25)
                    tgtgCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCCCGCCAACCAGGGTCCTGTGGGGGCGGGTTTa                                                                                                                                    >  1:828041/1‑150 (MQ=255)
                          gatgaCCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCgg                                                                                                                                         >  2:3163673/1‑139 (MQ=37)
                          gatgaCCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCgg                                                                                                                                         <  2:1891244/139‑1 (MQ=37)
                          gatgaCCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCgg                                                                                                                                         >  1:1891244/1‑139 (MQ=37)
                          gatgaCCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCgg                                                                                                                                         <  1:3163673/139‑1 (MQ=37)
                           atgaCCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTGTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCGACCAGGGTCCGGTGGGAGCGGGTTTACCCGCGa                                                                                                                             >  2:1032555/1‑150 (MQ=38)
                            tgaCCTCCGCACGCTGAAACTGCGGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTCCAGGAACCTTGTAGGAGCGGGTTTACCCGCGaa                                                                                                                            <  1:4955519/150‑1 (MQ=25)
                              aCCTACGCCCGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGTTTTTTTTTTGGGGTTATTTTTCCTGTGCGGGCCCCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGTGCGGGTTTACCCGCGAaga                                                                                                                          <  2:1830479/150‑1 (MQ=38)
                              aCCTACGCACGCTGAAACGGCTGCCTGTTTCACGAAAACCCCGGCCACGGT‑GCGGGGGTTTTTTTTTGGGGGTTTTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGATCCAACAAGGTGCGTGTAGTAGGGGGTTTACCCGCGAaga                                                                                                                          <  2:2993281/150‑1 (MQ=37)
                                    gcacgcTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCGCCGACAAGGTTGCGGTTGGAGCGGGGTTTCACCCGGAGAGGCCAg                                                                                                                    >  2:1007907/1‑150 (MQ=255)
                                    gcacgcTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCCCTCCGAAAAGGTTGCAGTAGGAGCGGGTTTTCCCGCGACGAGGCCAg                                                                                                                    >  1:5426086/1‑150 (MQ=38)
                                          tGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTc                                                                                                                                                     >  2:3437973/1‑111 (MQ=255)
                                          tGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTc                                                                                                                                                     <  1:3437973/111‑1 (MQ=255)
                                              aCTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCggg                                                                                                                                        >  1:3330643/1‑120 (MQ=25)
                                              aCTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCggg                                                                                                                                        <  2:3330643/120‑1 (MQ=25)
                                                 gcggcCGGTTACACCAACACCCCGGCCAGGGG‑CCCGGGGTTTTTTTTTGGGGGTATTTTCCCTGGCCCGCCCTCTTCCCGGGAAAACCCGCCCCCCCCAGGTCCCTGTCGGCGCGGGTTTACCCGCGAAGAGGCCAGCACA‑Gaaaaataa                                                                                                      <  2:4592341/147‑1 (MQ=255)
                                                   tgcCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCGCCTACAAGGTTCCGGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACA‑Gaaaaataaaa                                                                                                    >  2:1662098/1‑150 (MQ=255)
                                                   tgcCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCCCTCCTACAAGGTTACTGTATGAGCGGGTTTAACCGCGAAGAGGCCAGCACA‑Gaaaaaaaaaa                                                                                                    >  2:3333794/1‑150 (MQ=38)
                                                             aCGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCggg                                                                                                                                        >  1:2176946/1‑105 (MQ=21)
                                                             aCGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGAtt                                                                                          >  2:3993422/1‑150 (MQ=255)
                                                             aCGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTT‑CTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCgg                                                                                                                                         >  2:1964702/1‑103 (MQ=255)
                                                             aCGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTT‑CTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCgg                                                                                                                                         <  1:1964702/103‑1 (MQ=255)
                                                                 aaaCCCCGG‑CACAGTGGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCgg                                                                                                                                         >  2:482084/1‑100 (MQ=37)
                                                                 aaaCCCCGG‑CACAGTGGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCgg                                                                                                                                         <  1:482084/100‑1 (MQ=37)
                                                                   aCCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGc                                                                                                                                           <  1:5177240/96‑1 (MQ=21)
                                                                   aCCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGc                                                                                                                                           >  2:5177240/1‑96 (MQ=21)
                                                                      ccGGCCCCTGG‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCGGGCCACTTCGCGGGAAACCCCGCTCCTACCAGGTGCCTGTCGGAGCGGGTTTCCCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCaaaaaa                                                                                 <  2:4173464/150‑1 (MQ=255)
                                                                       cGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCATAGCGGGTAAAACCGCGCCTACAAGGGTGCGGTAGGAGCGGGTTTAACCGCGAAGAGGCCACCCCA‑GAAAAATAAAAAAACCATATTTTAaaaaaaa                                                                                >  2:2188443/1‑149 (MQ=255)
                                                                                    cGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCAGATTTTCAAAAAATAGATATGAAAAg                                                                    >  2:4951937/1‑150 (MQ=255)
                                                                                       tttttttttttGGGGTTATTTTACCTGGGCCGGCCTCTTCGCGGGAAAACCCGCTCCACCACGTTTCCTGTAGGCGCGGGTTAACCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTc                                                                 <  1:3336432/148‑1 (MQ=37)
                                                                                          ttttttctGGGGGTATTTTGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCAACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCtaa                                                              <  1:5434776/150‑1 (MQ=255)
                                                                                           tttttttGGGGGTATTTTTCCTGTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCtaat                                                             <  1:877636/150‑1 (MQ=255)
                                                                                                 tGGGGGTATTTGGCCTGTGCCGGCCTCTTCGCGGGTAACCCCGCTCCACCACGGGTCCTTTAGGTGCGGTTTTACCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACa                                                       <  1:3507718/150‑1 (MQ=37)
                                                                                                    gtttttttttCGTGTGCCGGCCTCTTCGGGGGTAACCCGGCTCCAACACGTTTCCTTTTGTCGGGGTTTTACCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGa                                                    <  1:2502349/145‑1 (MQ=255)
                                                                                                                       gCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCGTGTAGGAGCGGGTTTACCCCCGAAGAGGCCCGCACA‑CAAAAAAAAAAATACCAGATTTTAAAAAAAAAGAACTGAAACGTTCGAATAAAACAGGAAAATATCCTTTACAGGCaa                                 >  2:2138341/1‑149 (MQ=18)
                                                                                                                       gCCTCTTCGCGGGTAAACCCCCTCCGACAGGAACCGGGTAGGAGCGGGTTTACCCGCGCAGAGGCCAGCCCA‑GAAAAATAAAAAAACCAGATTGTCAAAAAAAAGATCTGAAAAGTTCGAATAAAACAGGAAAGTATCCTTTACAGGccg                                 >  2:2138154/1‑148 (MQ=21)
                                                                                                                        ccTCTTCGCGGTAAAACCCGCTCCTACAAGGTTCGTGTTGGAGCGGTTTTACCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATc                                <  1:4340252/150‑1 (MQ=38)
                                                                                                                            ttCGCGGGTAAACCCGCTCCAACAGGAACCGTGTAGGAGCGGGTTTAACCGCGAAGAGGCCAGCACA‑GAAAAAAAAAAATACCAGAATTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGt                            >  2:3402035/1‑150 (MQ=25)
                                                                                                                               gcgGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATc                                             <  2:4415141/130‑1 (MQ=255)
                                                                                                                               gcgGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTc                         >  1:1360847/1‑150 (MQ=255)
                                                                                                                                cgGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCg                        <  1:3710269/150‑1 (MQ=255)
                                                                                                                                   gTAAACCCGCTCCTACATGGTTCCTATAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCa                     <  2:2066803/150‑1 (MQ=255)
                                                                                                                                    tAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAAGTATCCTTTACAGGCATCCTGTTTCGGCAc                    >  1:1437512/1‑150 (MQ=255)
                                                                                                                                       aCCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAt                 >  2:4297839/1‑150 (MQ=255)
                                                                                                                                        cccGCTCCTACAGGAACCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtc                >  2:884061/1‑150 (MQ=38)
                                                                                                                                        cccGCTCCTACAGGAACCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtc                >  1:115429/1‑150 (MQ=38)
                                                                                                                                        cccGCTCCTACAGGAACCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtc                >  2:4086511/1‑150 (MQ=38)
                                                                                                                                        cccGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTc                         >  2:3071622/1‑142 (MQ=255)
                                                                                                                                        cccGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTc                         <  1:3071622/142‑1 (MQ=255)
                                                                                                                                        cccGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACA‑GAAAAATAACAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATGAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtc                >  1:2128316/1‑150 (MQ=255)
                                                                                                                                        cccGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtc                >  1:2676038/1‑150 (MQ=255)
                                                                                                                                         ccGCTCCTACAGGAACCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACA‑GAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtct               >  1:1761576/1‑150 (MQ=37)
                                                                                                                                         ccGCTCCTACAGGAACCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTa                                               <  2:2295738/118‑1 (MQ=25)
                                                                                                                                         ccGCTCCTACAGGAACCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTa                                               >  1:2295738/1‑118 (MQ=25)
                                                                                                                                         ccGCTCCTACAGGAACCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACa                                      >  1:2932731/1‑127 (MQ=37)
                                                                                                                                         ccGCTCCTACAGGAACCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCActctct               >  1:4405234/1‑150 (MQ=37)
                                                                                                                                         ccGCTCCTACAGGAACCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtct               >  2:1241656/1‑150 (MQ=38)
                                                                                                                                         ccGCTCCTACAGGAACCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtct               >  2:1241273/1‑150 (MQ=38)
                                                                                                                                         ccGCTCCTACAGGAACCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtct               >  2:4041913/1‑150 (MQ=38)
                                                                                                                                         ccGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCaaa                                                                                    >  1:3966032/1‑81 (MQ=255)
                                                                                                                                         ccGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCaaa                                                                                    <  2:3966032/81‑1 (MQ=255)
                                                                                                                                         ccGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtct               >  2:4704243/1‑150 (MQ=255)
                                                                                                                                         ccGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtct               >  1:3365591/1‑150 (MQ=255)
                                                                                                                                         ccGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtct               >  1:4735091/1‑150 (MQ=255)
                                                                                                                                           gCTCCTACAGGAACCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtctt              >  1:770506/1‑150 (MQ=25)
                                                                                                                                           gCTCCTACAGGAACCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACA‑GAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtcttc             >  2:976561/1‑150 (MQ=37)
                                                                                                                                           gCTCCTACAGGAACCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtcttc             >  2:4210138/1‑150 (MQ=38)
                                                                                                                                           gCTCCTACAGGAACCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtcttc             >  2:2193331/1‑150 (MQ=38)
                                                                                                                                              ccTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGa          <  1:3993422/150‑1 (MQ=255)
                                                                                                                                              ccTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGa          <  2:1782552/150‑1 (MQ=255)
                                                                                                                                                tACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGttt                          <  2:3755471/132‑1 (MQ=255)
                                                                                                                                                tACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGAtt        >  2:4794324/1‑150 (MQ=255)
                                                                                                                                                tACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGAtt        >  2:4794411/1‑150 (MQ=255)
                                                                                                                                                tACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGAtt        >  1:1316695/1‑150 (MQ=255)
                                                                                                                                                  cAAGGTCCCTGTCGGAGCGGCTTTCCCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCa      <  1:3333794/150‑1 (MQ=255)
                                                                                                                                                   ctggTTCCGGTTGGTGCGGGTTGACCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCAg     <  1:5122558/148‑1 (MQ=255)
                                                                                                                                                   aaGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtcttc             <  2:4592684/142‑1 (MQ=255)
                                                                                                                                                   aaGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtcttc             >  1:4592684/1‑142 (MQ=255)
                                                                                                                                                      |                                                                                                                                                      
GACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACA‑GAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCAGGCC  >  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: ‑