Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,788,671 A→T 6.6% intergenic (+22/‑314) ygaQ → / → csiD putative uncharacterized protein YgaQ/PF08943 family protein CsiD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,788,6710AT 6.6% 559.7 / 15.8 176intergenic (+22/‑314)ygaQ/csiDputative uncharacterized protein YgaQ/PF08943 family protein CsiD
Reads supporting (aligned to +/- strand):  ref base A (83/81);  new base T (4/8);  total (87/89)
Fisher's exact test for biased strand distribution p-value = 3.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

Only 100 of 177 total aligned reads displayed.

CCGAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCATTTCCTGTCTCACTTTGCC  >  NC_000913/2788523‑2788820
                                                                                                                                                    |                                                                                                                                                     
ccGAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGAt                                                                                                                                                      >  2:133657/1‑150 (MQ=255)
ccGAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGAt                                                                                                                                                      >  2:323285/1‑150 (MQ=255)
 cGAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATa                                                                                                                                                     <  2:1233752/150‑1 (MQ=255)
 cGAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATATGCCCCGATa                                                                                                                                                     >  1:1079552/1‑150 (MQ=255)
  gAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATaa                                                                                                                                                    >  1:1481101/1‑150 (MQ=255)
   aaaaaTTATCCCCCCTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATaaa                                                                                                                                                   <  2:2060614/150‑1 (MQ=255)
    aaaaTTATCCCCCCTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATaaaa                                                                                                                                                  >  2:1394886/1‑150 (MQ=255)
     aaaTTATCCCCCCTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAAc                                                                                                                                                 <  1:2245567/150‑1 (MQ=255)
     aaaTTATCCCCCCTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAAc                                                                                                                                                 >  1:489965/1‑150 (MQ=255)
         tATCCCCCCTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACgggg                                                                                                                                             <  2:1481101/150‑1 (MQ=255)
          aTCCCCCCTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGc                                                                                                                                            <  2:1717168/150‑1 (MQ=255)
           tCCCCCCTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGTTTTATCGGGGCa                                                                                                                                           <  2:1152990/150‑1 (MQ=255)
              ccccTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGTTTTATCGGGGCAGat                                                                                                                                        <  1:307752/150‑1 (MQ=255)
              ccccTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGat                                                                                                                                        >  1:1602315/1‑150 (MQ=255)
               cccTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGata                                                                                                                                       >  1:648681/1‑150 (MQ=255)
               cccTGGTAAAGGAGTAACGCCTAGAGGCTGAACGCTAAATAACAATATTCCTGGGCGTAGCGCTCGTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCGCTATTTATAATTAAATGTTATATCTGCCCCGTTTTATCGGGGCAGata                                                                                                                                       <  2:1333931/150‑1 (MQ=255)
                ccTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGataa                                                                                                                                      <  1:810641/150‑1 (MQ=255)
                    gTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATg                                                                                                                                  >  1:1512852/1‑150 (MQ=255)
                          tATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATgtttagt                                                                                                                            <  2:578416/150‑1 (MQ=255)
                              tACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACTGGGCAGATAATATGTTTAGTTTAc                                                                                                                        >  1:1661227/1‑150 (MQ=255)
                              tACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTAc                                                                                                                        >  2:1772462/1‑150 (MQ=255)
                              tACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTAc                                                                                                                        >  1:1950072/1‑150 (MQ=255)
                              tACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTAc                                                                                                                        >  1:1959728/1‑150 (MQ=255)
                              tACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTAc                                                                                                                        >  1:204267/1‑150 (MQ=255)
                                cGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTa                                                                                                                      >  1:2326068/1‑150 (MQ=255)
                                cGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTa                                                                                                                      >  1:2317757/1‑150 (MQ=255)
                                    tAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACgg                                                                                                                  <  1:1091049/150‑1 (MQ=255)
                                    tAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACgg                                                                                                                  <  1:1964382/150‑1 (MQ=255)
                                        ggCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCAt                                                                                                              >  1:754107/1‑150 (MQ=255)
                                        ggCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCAt                                                                                                              <  1:1394886/150‑1 (MQ=255)
                                          cGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCAttt                                                                                                            >  2:2010011/1‑150 (MQ=255)
                                           ggAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCAtttt                                                                                                           >  1:2295567/1‑150 (MQ=255)
                                           ggAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCAtttt                                                                                                           <  1:117896/150‑1 (MQ=255)
                                           ggAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCAtttt                                                                                                           >  2:1917587/1‑150 (MQ=255)
                                                   aaaTATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGaa                                                                                                   >  1:2350238/1‑150 (MQ=255)
                                                       atCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCa                                                                                               <  2:167732/150‑1 (MQ=255)
                                                       atCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCa                                                                                               >  2:1279759/1‑150 (MQ=255)
                                                       atCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCa                                                                                               <  1:2133044/150‑1 (MQ=255)
                                                        tCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCAt                                                                                              >  2:297806/1‑150 (MQ=255)
                                                         cAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCAtt                                                                                             >  1:1759629/1‑150 (MQ=255)
                                                           atatTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTa                                                                                           <  1:2185785/150‑1 (MQ=255)
                                                               tCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGt                                                                                       >  2:1151792/1‑150 (MQ=255)
                                                               tCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGt                                                                                       >  1:1699207/1‑150 (MQ=255)
                                                                ccTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGtt                                                                                      >  2:1539627/1‑150 (MQ=255)
                                                                ccTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGtt                                                                                      >  1:436707/1‑150 (MQ=255)
                                                                ccTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGtt                                                                                      <  1:1248161/150‑1 (MQ=255)
                                                                  tGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGtttt                                                                                    <  2:2174061/150‑1 (MQ=255)
                                                                  tGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGtttt                                                                                    >  2:1607522/1‑150 (MQ=255)
                                                                   gggCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACTGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGttttt                                                                                   <  2:1661227/150‑1 (MQ=255)
                                                                   gggCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGttttt                                                                                   <  1:2445125/150‑1 (MQ=255)
                                                                      cGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTAt                                                                                <  1:2170400/150‑1 (MQ=255)
                                                                      cGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCAATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTAt                                                                                <  1:1430691/150‑1 (MQ=255)
                                                                           cgcTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTTGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGAcc                                                                           >  2:857836/1‑150 (MQ=255)
                                                                              tcttctTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGa                                                                        <  1:1313934/150‑1 (MQ=255)
                                                                              tcttctTTTAGGGAAAACACGACAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGa                                                                        <  1:669008/150‑1 (MQ=255)
                                                                                 tcttTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATaa                                                                     <  1:1462183/150‑1 (MQ=255)
                                                                                 tcttTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATaa                                                                     >  1:1531016/1‑150 (MQ=255)
                                                                                 tcttTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATaa                                                                     >  1:1569443/1‑150 (MQ=255)
                                                                                  cttTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAAt                                                                    >  2:1526074/1‑150 (MQ=255)
                                                                                       aGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTATTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGtt                                                               >  1:1638573/1‑150 (MQ=255)
                                                                                        gggAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTtc                                                              <  2:2051026/150‑1 (MQ=255)
                                                                                          gAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTtctc                                                            >  2:1244612/1‑150 (MQ=255)
                                                                                          gAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATATGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTtctc                                                            >  2:249718/1‑150 (MQ=255)
                                                                                                 ggAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGtt                                                     >  2:1492806/1‑150 (MQ=255)
                                                                                                 ggAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACAAGATAATCTGTTCTCTAATGtt                                                     >  2:140501/1‑150 (MQ=255)
                                                                                                  gagaACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTa                                                    <  1:695516/150‑1 (MQ=255)
                                                                                                  gagaACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTa                                                    <  1:2410153/150‑1 (MQ=255)
                                                                                                  gagaACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTa                                                    <  2:2084496/150‑1 (MQ=255)
                                                                                                  gagaACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTa                                                    <  2:2315206/150‑1 (MQ=255)
                                                                                                       ccgccgAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTcc                                               <  2:1036909/150‑1 (MQ=255)
                                                                                                        cgccgAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTccc                                              >  1:1578422/1‑150 (MQ=255)
                                                                                                        cgccgAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTccc                                              <  2:2295567/150‑1 (MQ=255)
                                                                                                         gccgAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTcccc                                             >  1:2330242/1‑150 (MQ=255)
                                                                                                          ccgAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTccccc                                            <  2:1541392/150‑1 (MQ=255)
                                                                                                              aCTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAAc                                        <  2:2326068/150‑1 (MQ=255)
                                                                                                               cTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAAcc                                       <  2:340984/150‑1 (MQ=255)
                                                                                                                 aTCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTg                                     >  1:747064/1‑150 (MQ=255)
                                                                                                                 aTCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAATTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTg                                     >  2:979532/1‑150 (MQ=255)
                                                                                                                     ctATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGc                                 <  2:1959728/150‑1 (MQ=255)
                                                                                                                     ctATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGc                                 <  2:1950072/150‑1 (MQ=255)
                                                                                                                      tATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGGCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCt                                <  1:269117/150‑1 (MQ=255)
                                                                                                                        tttATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCttt                              <  2:754107/150‑1 (MQ=255)
                                                                                                                         ttATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTa                             <  1:444345/150‑1 (MQ=255)
                                                                                                                           ataATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTTAAGCCATTTACTTTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGt                           <  2:609905/150‑1 (MQ=255)
                                                                                                                            taATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGtt                          <  1:1780039/150‑1 (MQ=255)
                                                                                                                             aaTTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGttt                         >  1:40290/1‑149 (MQ=255)
                                                                                                                             aaTTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTa                         <  2:128332/150‑1 (MQ=255)
                                                                                                                                ttAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTAtta                      >  2:1983293/1‑149 (MQ=255)
                                                                                                                                ttAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTc                      <  2:402157/150‑1 (MQ=255)
                                                                                                                                    tGTTATATCTGCCCAGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCAttt                  >  2:629940/1‑150 (MQ=255)
                                                                                                                                       tataTCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCATTTCCt               >  2:827052/1‑150 (MQ=255)
                                                                                                                                        atatCTGGCCCGATAAAACGGGGCAGATAATATGATTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGGTATTCATTTCCTg              <  1:553596/150‑1 (MQ=255)
                                                                                                                                         tatCTGCCCCGTTTTATCGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCATTTCCTGt             >  1:973392/1‑150 (MQ=255)
                                                                                                                                           tCTGCCCCGTTTTATCGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAg                            >  2:951467/1‑133 (MQ=255)
                                                                                                                                           tCTGCCCCGTTTTATCGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCATTTCCTGtct           >  2:2480586/1‑150 (MQ=255)
                                                                                                                                           tCTGCCCCGTTTTATCGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCATTTCCTGtct           >  1:931991/1‑150 (MQ=255)
                                                                                                                                            cTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCATTTCCTGtctc          <  1:1335903/150‑1 (MQ=255)
                                                                                                                                                cccGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCATTTCCTGTCTCACtt      >  2:252344/1‑150 (MQ=255)
                                                                                                                                                cccGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCATTTCCTGTCTCACtt      >  1:1880902/1‑150 (MQ=255)
                                                                                                                                                    aTAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCATTTCCTGTCTCACTTTGcc  <  2:2203674/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CCGAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCATTTCCTGTCTCACTTTGCC  >  NC_000913/2788523‑2788820

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