Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,511,584 T→A 9.1% intergenic (+46/‑153) PP_4839 → / → cycA iron‑regulated membrane protein/D‑alanine/beta‑alanine/D‑serine/glycine permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,511,5840TA 9.1% 379.6 / 2.2 152intergenic (+46/‑153)PP_4839/cycAiron‑regulated membrane protein/D‑alanine/beta‑alanine/D‑serine/glycine permease
Reads supporting (aligned to +/- strand):  ref base T (72/57);  new base A (13/3);  total (90/62)
Fisher's exact test for biased strand distribution p-value = 6.18e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 158 total aligned reads displayed.

TGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCGTGTTAGCCTAGCCGGCTACCTCTGACAAGAAGACTGCAG  >  NC_002947/5511435‑5511733
                                                                                                                                                     |                                                                                                                                                      
tggtgctggGGGTGATGGCCCCGCTGGGGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCt                                                                                                                                                        <  2:7974426/150‑1 (MQ=255)
cggcgctggGGGTGATGGTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCt                                                                                                                                                        <  2:7974672/146‑1 (MQ=255)
  gggcggggggTGATGGGCCCGCCGGCGGGGGTGTCGATGGTGGTGATGGGGGGGGTGGATAGTTTGGTGGGGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTAc                                                                                                                                                      <  2:8607191/145‑1 (MQ=38)
   tgctggGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACg                                                                                                                                                     >  1:7891048/1‑150 (MQ=255)
     cggggggCGATGGTCCCGCCGGCGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAg                                                                                                                                                   <  1:3221660/148‑1 (MQ=255)
           gTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGGGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCACCAAACAGGGACcg                                                                                                                                             >  2:4408473/1‑150 (MQ=39)
               tGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGtt                                                                                                                                         <  2:121734/150‑1 (MQ=255)
                gTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCCCTCCAACGACGGAcccccgtg                                                                                                                                        >  1:2665397/1‑144 (MQ=39)
                 ttCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGt                                                                                                                                       >  1:4323772/1‑150 (MQ=255)
                 gtCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGt                                                                                                                                       <  2:4937691/149‑1 (MQ=255)
                 gtCCCGCTGGTGGGGGGGGCGAGGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGt                                                                                                                                       <  2:2279699/149‑1 (MQ=255)
                 gtCCCGCTGGCGGGGGTGTCGATGGTGGTGATGGGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGttct                                                                                                                                       <  2:6451339/149‑3 (MQ=255)
                 ggcccGCTGGGGGGGGGGTCGATGGGGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGt                                                                                                                                       <  2:4937678/148‑1 (MQ=255)
                   cccGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAAACCCCACCAAAGACGGACCCCGTTGTAg                                                                                                                                     >  1:7056900/1‑150 (MQ=255)
                   cccGCGGGGGGGGGGGTCGATGGTGGGGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAg                                                                                                                                     <  2:5516881/150‑1 (MQ=255)
                     cGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAGGa                                                                                                                                   >  2:1565291/1‑150 (MQ=255)
                     cGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTACTACGAGGGACCGCGTTGTAGGa                                                                                                                                   >  2:3185969/1‑150 (MQ=255)
                     cGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAAACCCCACCAAAGAGGGACCCCGTGGTAGGa                                                                                                                                   >  1:6956837/1‑150 (MQ=255)
                       cTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGc                                                                                                                                 <  2:142642/150‑1 (MQ=255)
                         ggtggGGGTGTCGAGGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGCgg                                                                                                                               <  1:7713376/150‑1 (MQ=255)
                          gggggggCGTCGAGGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGCggg                                                                                                                              <  2:7486294/148‑1 (MQ=255)
                             ggggTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGCGGGttt                                                                                                                           >  2:1985988/1‑150 (MQ=255)
                               ggTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCCCGCCAAAGAGGGACCGCGGTGGAGGAGCGGGGTTAc                                                                                                                        >  1:508806/1‑150 (MQ=255)
                               ggTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAAACCCCTCCAACGAGGGACCGCGTTGTAGGAGCGGGTTTAc                                                                                                                        >  1:3345138/1‑150 (MQ=255)
                                      ctggtggGGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTGGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTCCGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGaa                                                                                                                 <  2:6529231/149‑1 (MQ=255)
                                         gtggtgATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGGTGTAGGAGCGGGTTTACCCGCGAagag                                                                                                              >  2:1027553/1‑150 (MQ=255)
                                         gtggtgATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTGGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGCGGGGATACCCGCGAacag                                                                                                              >  2:5560511/1‑150 (MQ=255)
                                          tggtgATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAAACCGCTACAACGAGGGACCGCGTTGTATGAGCGGGTTTACCCGCGAAGAgg                                                                                                             >  2:287145/1‑150 (MQ=255)
                                            gggagggggggggTGGATAGTTTGGGGGTGCGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAACCCCGCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGcc                                                                                                           <  2:7140573/143‑1 (MQ=38)
                                             tgctggggggggTGTTTGTTTTGGGGGTGCGGCGCCGCGGGTTCTTGGCCAGGGTTTGCGCGTTTTTTTGCGTGTGCTGCCCTCTTGGGGGGTAAACCGGCACCACCGCGGGACGGGTTTTTTGGTGGGTTTTTACCCGCGAAGAGGCCg                                                                                                          <  2:8088455/144‑1 (MQ=255)
                                                   ggggggTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAg                                                                                                    >  2:6098469/1‑150 (MQ=255)
                                                    gggggTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAgg                                                                                                   >  2:4765640/1‑150 (MQ=255)
                                                    gggggTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCCCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAgg                                                                                                   >  1:5761152/1‑150 (MQ=255)
                                                      gggTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGCgg                                                                                                                               >  1:1614231/1‑121 (MQ=255)
                                                         tGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGaaa                                                                                              >  2:7095814/1‑150 (MQ=255)
                                                             tAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTAAGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCt                                                                                          >  2:3265446/1‑150 (MQ=255)
                                                               gtgtgggggTGCGGCGCCGCCGGGCCTGGGCCAGGGCTTGAGCCTTGTGTGGCCTGTGCTGGCCTATGCGCGGGAAAACCCGCTCCTACGTGGGACCGCGTTGTAGCAGGGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGa                                                                                        <  2:1026748/147‑1 (MQ=37)
                                                               gTTGGGGGGTGCGGCGCCGCCGGGTCTGGGCCCGCGCGGGCGCCTGGTGGGGCCTGGGCGGGCCTCTGCGCGGGTAAACCCGCTCCTCCGAGGGCCCGCGTTGTCGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGa                                                                                        <  2:2696547/150‑1 (MQ=21)
                                                                tttGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAAAACGCACCTAAGAGGGACCCCCGTGTTGGTGCGGGTTTACCCCCGAAGAGGCCCGTAAAGGCGAAAAGCTGAg                                                                                       >  1:2275954/1‑150 (MQ=37)
                                                                 tgggtggtGAGGCGCCGCCGGGTCTGGGCCAGCGCGTGCGCCTTGTGGTGCCTGGGCTGGCCTGTGCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGc                                                                                      <  2:8624465/148‑1 (MQ=37)
                                                                    gtggtgAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTGCTAAGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATTCTGAGCCGa                                                                                   >  1:5723340/1‑150 (MQ=255)
                                                                      ggtgAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAAACCCCACCTACCAGGGACCGCGGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAAAGCGGAGCCGaaa                                                                                 >  1:7724977/1‑150 (MQ=38)
                                                                            gcgcCGCGGGGTCTTGGCCAGGGGTTGCGCCTTGTGTTGCCTGTGGTGGCCTCTGCGCGGTTAACCCCGCTCCTACGAGGGACCGCGTTGTTGGTGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCag                                                                           <  2:4061741/150‑1 (MQ=37)
                                                                            gcgcCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGCGGTTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCag                                                                           <  2:2756551/150‑1 (MQ=255)
                                                                            gcgcCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCCCTCCGACGAGGGCCCGCGGTGTTGGGGCGGGTTTACCCCCGCAGCGGCCGGTACAGGCGGAATGCGGGGGCGAAAATGccg                                                                           >  1:6008172/1‑148 (MQ=37)
                                                                             ggcggcggggTCTTGGCCAGGGTTTGAGCTTTGTGTTGCGGGTGGTGGCCTCTTCGGGGGTAAACCCGCTCCTACTTGGGACCGGTTTGTTGGAGGGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCaga                                                                          <  2:2304778/143‑1 (MQ=25)
                                                                               ccgccgGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGt                                                                        >  1:6547481/1‑150 (MQ=255)
                                                                                cgccgGGTGTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTc                                                                       <  2:5737355/150‑1 (MQ=255)
                                                                                  ccgGGTCTGGGCCCGCGCTTGAGCCTGGGGGTGCGGGGGCTGGCCTGTGCGCGGTTAAACCCGCGCCTACGAGGGACCGCGTTGTCGGAGGGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGa                                                                     <  2:1620264/150‑1 (MQ=25)
                                                                                      gTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCGCCTAAGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAAAGCTGAGCCGAAAATGCAGAGTCGATATg                                                                 >  1:2822671/1‑150 (MQ=255)
                                                                                         ttGGGCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGTGGGACCGCGTTGTAGGAGCGGGTTTAGCCGCGAAGTGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTcg                                                              >  2:8796888/1‑150 (MQ=255)
                                                                                             cccgggcTTGAGCTTTGTGTGGCCTGGGCGGGCCTCTTCGCGGTTAAACCCGCCCCTCGGTGGGACCGCGTTGTTGGTGCGGTTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGcc                                                          <  1:5638872/145‑1 (MQ=37)
                                                                                                  gcTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCCCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAAATCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGcccccct                                                     >  2:3850711/1‑147 (MQ=255)
                                                                                                   cTTGAGCGTTGTGTGGCCGGGGCTGGCCTCTTCGCGGGTAACCCCGCTCCGACGAGGGACCGCGTTGTCGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATc                                                    <  1:5705546/150‑1 (MQ=38)
                                                                                                     tgcgCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTTTAGGAGGGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTg                                                  <  2:946919/147‑1 (MQ=255)
                                                                                                     tGAGCCTTGTGTTGCCTGTGCGGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTg                                                  <  2:5379331/150‑1 (MQ=255)
                                                                                                       aGCGTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTCGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGa                                                <  2:6201579/150‑1 (MQ=255)
                                                                                                           ttGTGTTGCCTGTGCTGGCCTCTTCGGGGGTAAACCCGCTCCTACGTGGGACCGCGTTGTCGGAGGGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAAcgc                                            <  1:8602682/150‑1 (MQ=255)
                                                                                                           ttGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGtt                                                                                                                                         <  2:301441/58‑1 (MQ=255)
                                                                                                           ttGTGTGGCCTGGGCTGGCCTCTGCGGGGGTAAACCCGCCCCAACGCGGGACCGGGTTGTTGGAGCGTGTTAACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAAcgc                                            <  1:2769337/150‑1 (MQ=255)
                                                                                                             gtgtTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAAcgcgc                                          <  2:2819970/150‑1 (MQ=255)
                                                                                                                  gCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGCGGGACCGCGGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGAGCGGGCACGCGCGTGtt                                     >  2:8264459/1‑150 (MQ=38)
                                                                                                                  gCCTGTGCTGGCCTCTTCGCGGGGTAACCCGCTCCGACGAGGGGCCGCGGTGTAGGAGCGGGGTTACCCGCGAAGAGGCCGGTACAGGCGAAATTCTGAGCCGAAAATGCAGAGTCGATCAGGCGCGCCCCCATCTGGGACGCCCGGGtt                                     >  1:4991244/1‑150 (MQ=25)
                                                                                                                     tgtaccggCCTCTTCGCGGGTAAACCCGCTCCTACAACGCGGTCCCTCGTAGGAGCGGGTTTTCCCGCGGAGAGGCCCGTACAGGCGAAAATCTGAGCCGAAAATGCAGATTAGATCTGTCGCGCCCCCATCTGGGACGCGCGTGTtagc                                  >  2:1243818/7‑150 (MQ=9)
                                                                                                                     ggcgggggCCTCTCCGGGGGTAACCCCGCCCCACCGGGGCCCGGGGTTTTGGGCGGGGTTTACCCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCGTGTtagc                                  <  1:1825486/144‑1 (MQ=255)
                                                                                                                       tgCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGCgg                                                                                                                               <  2:3119484/56‑1 (MQ=255)
                                                                                                                       tgCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGCgg                                                                                                                               <  2:3119827/56‑1 (MQ=255)
                                                                                                                       tgCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGCgg                                                                                                                               >  1:3119484/1‑56 (MQ=255)
                                                                                                                       tgCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGCgg                                                                                                                               >  1:3119827/1‑56 (MQ=255)
                                                                                                                        gCTGGCCTCTTCGCGGGTAAACCCGCTCCTCCGAGGGACCGCGTTGTCGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCGTGTtagccta                               <  2:587870/150‑1 (MQ=255)
                                                                                                                          tGGCCTCTTCGCGGGTAAAACCGCTCCAACCAGGGACCGCGTTGTAGGAGCGGGTTTTCCCGCGAAGAGGCCCGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTAGCGCCCCGCGCTGGAACGCGCGTGTtagccaagc                             >  2:2384821/1‑150 (MQ=38)
                                                                                                                          tGGCCTCTTCGCGGGTAAAACCGCAACTAACAGGGACCGCGTTGTTGGAGCCGGTTTAACCCCCAAGTGGCCCGTAAACGCGAAATTCTGAGCCGAAAATTCAGAGTCGATCTGTCGCGCCCCGCTCTGGAAAGCGCGGGTtagccgagc                             >  2:2547485/1‑150 (MQ=255)
                                                                                                                              ctctTCGCGGGTAAAACCGCTCCTACGAGGGACCGCGGTGTAGGAGCGGGGTTACCCGCGAAGAGGCCCGTACAGGCGAAATGCTTAGCCGAAAAATCAGAGTAGATCTGTCGCGCCCCCATCAGGAACGCGCGTGTTATCCTAGCCGGc                         >  2:3700160/1‑150 (MQ=255)
                                                                                                                               tctTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCGTGTTAGCCGAGCCGGCt                        >  2:6371822/1‑150 (MQ=255)
                                                                                                                               tctTCGCGGGTAAACCCGCTCCAACGAGGGACCGCGGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAAAGCTGAGCCCAAAATGCAGAGTCGATCTGTCGCCCCCCGCTCTGGAACCCGCGTGTTAGCCGAGCCGGca                        >  2:3204004/1‑149 (MQ=37)
                                                                                                                                      gggTAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCGTGTTAGCCTAGCCGGCTACCTCTg                 >  2:2731754/1‑150 (MQ=255)
                                                                                                                                        gTAAACCCGCTCCTACAACGCGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCgt                                        <  2:3266315/125‑1 (MQ=18)
                                                                                                                                         tAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCGTGTTAGCCTAGCCGGCTACCTCTGACa              >  1:8048421/1‑150 (MQ=255)
                                                                                                                                         tAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCGTGTTAGCCTAGCCGGCTACCTCTGACa              >  2:4127916/1‑150 (MQ=255)
                                                                                                                                         tAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCGTGTTAGCCTAGCCGGCTACCTCTGACa              >  1:5236158/1‑150 (MQ=255)
                                                                                                                                          aaaCCCGCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCGTGTTAGCCTAGCCGGCTACCTCTGACaa             >  2:1894939/1‑150 (MQ=255)
                                                                                                                                             cccGCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCGTGTTAGCCTAGCCGGCTACCTCTGACaagaa          >  2:7751889/1‑150 (MQ=255)
                                                                                                                                             cccGCCCCTCCGCGGGCCCGGTTTGTCGGCGCGGGTTTCCCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCGTGTTAGCCTAGCCGGCTACCTCTGACaagaa          <  1:2393494/150‑1 (MQ=255)
                                                                                                                                              ccGCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCGTGTTAGCCTAGCCGGCTACCTCTGACaagaag         >  1:3466804/1‑150 (MQ=255)
                                                                                                                                              ccGCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCGTGTTAGCCTAGCCGGCTACCTCTGACaagaag         >  2:8816441/1‑150 (MQ=255)
                                                                                                                                              ccGCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCGTGTTAGCCTAGCCGGCTACCTCTGACaagaag         >  2:8142341/1‑150 (MQ=255)
                                                                                                                                              ccGCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCGTGTTAGCCTAGCCGGCTACCTCTGACaagaag         >  1:3466632/1‑150 (MQ=255)
                                                                                                                                              ccGCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCGTGTTAGCCTAGCCGGCTACCTCTGACaagaag         >  2:6581608/1‑150 (MQ=255)
                                                                                                                                              ccGCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCGTGTTAGCCTAGCCGGCTACCTCTGACaagaag         >  2:6569021/1‑150 (MQ=255)
                                                                                                                                              ccGCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCGTGTTAGCCTAGCCGGCTACCTCTGACaagaag         >  1:5002062/1‑150 (MQ=255)
                                                                                                                                              ccGCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCGTGTTAGCCTAGCCGGCTACCTCTGACaagaag         >  1:6588197/1‑150 (MQ=255)
                                                                                                                                               cGCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCGTGTTAGCCTAGCCGGCTACCTCTGACAagaaga        >  2:7709268/1‑150 (MQ=255)
                                                                                                                                                gCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCGTGTTAGCCTAGCCGGCTACCTCTGACAAGAAGAc       >  1:1896822/1‑150 (MQ=255)
                                                                                                                                                  tCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCGTGTTAGCCTAGCCGGCTACCTCTGACAAGAAGACTg     >  1:4803640/1‑150 (MQ=255)
                                                                                                                                                   ccTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCGTGTTAGCCTAGCCGGCTACCTCTGACAAGAAGACTGc    >  2:3846762/1‑150 (MQ=255)
                                                                                                                                                    cTACGCGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCGTGTTAGCCTAGCCGGCTACCTCTGACAAGAAGACTGCa   <  1:7130177/150‑1 (MQ=255)
                                                                                                                                                    cTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCGTGTTAGCCTAGCCGGCTACCTCTGACAATAAGACTGCa   >  1:3163924/1‑150 (MQ=255)
                                                                                                                                                    cTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCGTGTTAGCCTAGCCGGCTACCTCTGACAAGAAGACTGCa   >  2:4613264/1‑150 (MQ=255)
                                                                                                                                                    cTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCGTGTTAGCCTAGCCGGCTACCTCTGACAAGAAGACTGCa   >  1:1498265/1‑150 (MQ=255)
                                                                                                                                                     tACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCGTGTTAGCCTAGCCGGCTACCTCTGACAAGAAGACTGCAg  >  1:446308/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
TGGTGCTGGGGGTGATGTTCCCGCTGGTGGGGGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCGTGTTAGCCTAGCCGGCTACCTCTGACAAGAAGACTGCAG  >  NC_002947/5511435‑5511733

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