Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,590,771 C→G 6.1% intergenic (+63/+36) acrD → / ← ypfM multidrug efflux pump RND permease AcrD/uncharacterized protein YpfM

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,590,7710CG 6.1% 326.9 / 8.7 99intergenic (+63/+36)acrD/ypfMmultidrug efflux pump RND permease AcrD/uncharacterized protein YpfM
Reads supporting (aligned to +/- strand):  ref base C (52/40);  new base G (2/4);  total (55/44)
Fisher's exact test for biased strand distribution p-value = 4.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.75e-01

GGCCACTATTCTGGCTATTTACTTCGTGCCGCTGTTCTTTGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCGTGGTGAA  >  NC_000913/2590624‑2590920
                                                                                                                                                   |                                                                                                                                                     
ggCCACTATTCTGGCTATTTACTTCGTGCCGCTGTTCTTTGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAg                                                                                                                                                     <  1:1395148/150‑1 (MQ=255)
 gCCACTATTCTGGCTATTTACTTCGTGCCGCTGTTCTTTGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGc                                                                                                                                                    <  1:908493/150‑1 (MQ=255)
   cACTATTCTGGCTATTTACTTCGTGCCGCTGTTCTTTGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATCTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCgg                                                                                                                                                  >  1:1194048/1‑150 (MQ=255)
          cTGGCTATTTACTTCGTGCCGCTGTTCTTTGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAATTCTACATGGCAATGTGTCGACTTTTTCAACTTTCCTATTAAGAACCTTCTCAGCGGGTTCTTg                                                                                                                                           <  2:779654/150‑1 (MQ=255)
           tGGCTATTTACTTCGTGCCGCTGTTCTTTGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGc                                                                                                                                          <  1:203868/150‑1 (MQ=255)
            ggCTATTTACTTCGTGCCGCTGTTCTTTGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCt                                                                                                                                         >  1:563361/1‑150 (MQ=255)
             gCTATTTACTTCGTGCCGCTGTTCTTTGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCCGCTGAGCAGGTTCttaatc                                                                                                                                        <  1:604510/150‑5 (MQ=255)
              cTATTTACTTCGTGCCGCTGTTCTTTGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCCGCTGAGCAGGTTCttaatcg                                                                                                                                       <  1:106018/150‑6 (MQ=255)
              cTATTTACTTCGTGCCGCTGTTCTTTGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCCGCTGAGCAGGTTCttaatcg                                                                                                                                       <  1:65041/150‑6 (MQ=255)
               tATTTACTTCGTGCCGCTGTTCTTTGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGtt                                                                                                                                      <  1:766533/150‑1 (MQ=255)
               tATTTACTTCGTGCCGCTGTTCTTTGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCCGCTGAGCAGGTTCttaatcgg                                                                                                                                      <  2:21237/150‑7 (MQ=255)
               tATTTACTTCGTGCCGCTGTTCTTTGTGCTGGTGCGCCGCCGCTTCCCGATGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGtt                                                                                                                                      >  1:340788/1‑150 (MQ=255)
                  ttACTTCGTGCCGCTGTTCTTTGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGt                                                                                                                                   <  1:567810/150‑1 (MQ=255)
                   tACTTCGTGCCGCTGTTCTTTGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTa                                                                                                                                  >  1:1213061/1‑150 (MQ=255)
                          tGCCGCTGTTCTTTGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGtc                                                                                                                           >  1:406052/1‑150 (MQ=255)
                           gCCGCTGTTCTTTGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGtct                                                                                                                          >  1:103335/1‑150 (MQ=255)
                            ccGCTGTTCTTTGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGtctc                                                                                                                         <  1:129502/150‑1 (MQ=255)
                             cGCTGTTCTTTGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCa                                                                                                                        <  1:1280799/150‑1 (MQ=255)
                              gCTGGTCTTTGTGCTGGTGCGGCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAg                                                                                                                       <  2:754205/150‑1 (MQ=255)
                               cTGTTCTTTGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAgg                                                                                                                      >  1:1303915/1‑150 (MQ=255)
                                  ttCTTTGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAAt                                                                                                                   <  1:1339163/150‑1 (MQ=255)
                                   tCTTTGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAAtt                                                                                                                  <  2:290480/150‑1 (MQ=255)
                                    cTTTGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAAtta                                                                                                                 >  2:1000764/1‑150 (MQ=255)
                                    cTTTGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAAtta                                                                                                                 >  2:952347/1‑150 (MQ=255)
                                     tttGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGtc                                                                                                                           <  2:1480308/139‑1 (MQ=255)
                                     tttGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGtc                                                                                                                           >  1:1480308/1‑139 (MQ=255)
                                      ttGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAAttact                                                                                                               <  1:622016/150‑1 (MQ=255)
                                         tgCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAAttacttac                                                                                                            <  2:421164/150‑1 (MQ=255)
                                          gCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACg                                                                                                           <  1:380210/150‑1 (MQ=255)
                                           cTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGa                                                                                                          <  1:481874/150‑1 (MQ=255)
                                            tGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGaa                                                                                                         <  1:1251096/150‑1 (MQ=255)
                                                gcgcCGCCGCTTCCCGCGGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCAt                                                                                                     <  1:931958/150‑1 (MQ=255)
                                                 ggccgccgcTTCCAGCTGAAGCCGCGCCCGGAATAAGCAATAATAAAGGAGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATa                                                                                                    <  2:499108/149‑1 (MQ=255)
                                                 cgccgccgcTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCATTGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATa                                                                                                    <  2:225380/150‑1 (MQ=255)
                                                 cgccgccgcTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATa                                                                                                    >  1:1365212/1‑150 (MQ=255)
                                                 cgccgccgcTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTCGTCTCAGGAATTACTTACGAAGCATa                                                                                                    <  1:1389442/150‑1 (MQ=255)
                                                  gccgccgcTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATaa                                                                                                   >  1:1540855/1‑150 (MQ=255)
                                                  gccgccgcTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATaa                                                                                                   >  2:1180969/1‑150 (MQ=255)
                                                  gccgccgcTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATaa                                                                                                   >  1:1381220/1‑150 (MQ=255)
                                                  gccgccgcTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATaa                                                                                                   >  2:536925/1‑150 (MQ=255)
                                                    cgccgcTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACt                                                                                                 >  1:745000/1‑150 (MQ=255)
                                                      ccgcTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTc                                                                                               <  2:1174355/150‑1 (MQ=255)
                                                         cTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGAt                                                                                            <  2:1213061/150‑1 (MQ=255)
                                                              cGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGt                                                                                       >  2:204669/1‑150 (MQ=255)
                                                                   aaGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTc                                                                                  >  1:1437822/1‑150 (MQ=255)
                                                                    aGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTcc                                                                                 <  1:1302804/150‑1 (MQ=255)
                                                                      ccGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAg                                                                               >  1:1308528/1‑150 (MQ=255)
                                                                       cgcgCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGt                                                                              <  1:952347/150‑1 (MQ=255)
                                                                       cgcgCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGt                                                                              >  1:1318408/1‑150 (MQ=255)
                                                                              ggAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCac                                                                       >  2:164214/1‑149 (MQ=255)
                                                                              ggAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAg                                                                       >  2:286542/1‑150 (MQ=255)
                                                                               gAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGt                                                                      <  2:103335/150‑1 (MQ=255)
                                                                                 aTAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTc                                                                    <  2:634314/150‑1 (MQ=255)
                                                                                  tAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTc                                                                    >  2:211940/1‑149 (MQ=255)
                                                                                  tAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTc                                                                    <  1:211940/149‑1 (MQ=255)
                                                                                      cAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGtt                                                               >  1:1132418/1‑150 (MQ=255)
                                                                                       aaTAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTc                                                              >  2:940005/1‑150 (MQ=255)
                                                                                        aTAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCa                                                             >  1:1557764/1‑150 (MQ=255)
                                                                                        aTAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCa                                                             >  1:29501/1‑150 (MQ=255)
                                                                                           aaaaaGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATc                                                          >  1:1322503/1‑150 (MQ=255)
                                                                                            aaaaGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCa                                                         <  2:948171/150‑1 (MQ=255)
                                                                                              aaGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATa                                                       >  2:128046/1‑150 (MQ=255)
                                                                                                 gCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATaaca                                                    >  2:852112/1‑150 (MQ=255)
                                                                                                       tGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTTATCAGCGGGTTATTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTACAGTTCAAAAGTTCACGTTAAATCATAACAACctct                                              >  2:488157/1‑150 (MQ=255)
                                                                                                       tGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACctct                                              >  1:808242/1‑150 (MQ=255)
                                                                                                        gCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACctctc                                             >  2:339756/1‑150 (MQ=255)
                                                                                                         ccAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACctctct                                            >  1:1627723/1‑150 (MQ=255)
                                                                                                         caaaTGTTTCGCATTTTTAAACTTTCCGATTAATAAACTGATAAGCGGGTTCTTGCTGTTTGTAAATATTCTCACCAAATACTTAAAAAAAATAAGTTGGATAAAGTATTTCCAGTTACCCAGGTAACGTGAAATCATAAACAActctct                                            >  2:715823/3‑150 (MQ=255)
                                                                                                                 tctgcTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGTGTTCTTGCTGTTTGTTCTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTa                                    <  1:328912/146‑1 (MQ=255)
                                                                                                                  cGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTAt                                   >  2:99730/1‑150 (MQ=255)
                                                                                                                  cGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTAt                                   >  1:831889/1‑150 (MQ=255)
                                                                                                                   tccTTTTTCAACTTTCCGATTAAGAACCTTCTCATCGGGTTCTTTGTTTTTGTAATTTTTCTCATGAATTACTTAAGAACCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATg                                  <  1:163027/149‑1 (MQ=255)
                                                                                                                   gCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATg                                  <  2:1437822/150‑1 (MQ=255)
                                                                                                                    ccTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATgg                                 >  1:1039346/1‑150 (MQ=255)
                                                                                                                    ccTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAAATTCGATAAAGAGTTTCCAGTTCCCAAGTTCACTTGAAATAATAACAACATATCTTTAAATTAtaa                                 >  2:163027/1‑148 (MQ=255)
                                                                                                                     cTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATggg                                <  2:1381220/150‑1 (MQ=255)
                                                                                                                     cTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATggg                                <  2:406052/150‑1 (MQ=255)
                                                                                                                       ttttCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTa                              >  1:1602063/1‑150 (MQ=255)
                                                                                                                       ttttCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTa                              <  2:1308528/150‑1 (MQ=255)
                                                                                                                         ttCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGACTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTAtt                            >  2:85623/1‑150 (MQ=255)
                                                                                                                             aCTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTAc                        >  2:798144/1‑150 (MQ=255)
                                                                                                                             aCTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTAc                        <  2:1322503/150‑1 (MQ=255)
                                                                                                                              cTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACg                       <  1:1387982/150‑1 (MQ=255)
                                                                                                                              cTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACg                       >  1:255859/1‑150 (MQ=255)
                                                                                                                                  aacagcaaGAACCCGCTGAGCAGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGTTATTAAACGaaaa                   >  1:441760/7‑150 (MQ=255)
                                                                                                                                   cGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGaaaa                  >  1:1188957/1‑149 (MQ=255)
                                                                                                                                   cGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAAc                  >  1:1194453/1‑150 (MQ=255)
                                                                                                                                     agcaaGAACCCGCTGAGCAGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACaa                >  2:1216598/4‑150 (MQ=255)
                                                                                                                                     aTTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACaa                >  1:859804/1‑150 (MQ=255)
                                                                                                                                       tAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAata              >  2:640562/1‑150 (MQ=255)
                                                                                                                                        aaGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAatat             >  1:1308318/1‑150 (MQ=255)
                                                                                                                                        aaGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAatat             >  2:129747/1‑150 (MQ=255)
                                                                                                                                         aGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAAtata            >  2:311263/1‑150 (MQ=255)
                                                                                                                                            aCCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCg         >  1:70250/1‑150 (MQ=255)
                                                                                                                                             ccTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCgt        <  1:488157/150‑1 (MQ=255)
                                                                                                                                               tGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCgtgg      <  2:354461/150‑1 (MQ=255)
                                                                                                                                                gCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCgtggt     >  2:1162867/1‑150 (MQ=255)
                                                                                                                                                  tCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCGTGGTGa   <  1:852112/150‑1 (MQ=255)
                                                                                                                                                  tCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCGTGGTGa   >  2:1427193/1‑150 (MQ=255)
                                                                                                                                                   cAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCGTGGTGaa  >  1:412673/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GGCCACTATTCTGGCTATTTACTTCGTGCCGCTGTTCTTTGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCGTGGTGAA  >  NC_000913/2590624‑2590920

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