Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,590,775 G→A 6.0% intergenic (+67/+32) 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,7750GA 6.0% 307.9 / 6.3 99intergenic (+67/+32)acrD/ypfMmultidrug efflux pump RND permease AcrD/uncharacterized protein YpfM
Reads supporting (aligned to +/- strand):  ref base G (52/40);  new base A (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 = 6.18e-01

Only 100 of 102 total aligned reads displayed.

CCACTATTCTGGCTATTTACTTCGTGCCGCTGTTCTTTGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCGTGGTGAAGC  >  NC_000913/2590626‑2590922
                                                                                                                                                     |                                                                                                                                                   
ccACTATTCTGGCTATTTACTTCGTGCCGCTGTTCTTTGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCg                                                                                                                                                     >  2:231817/1‑150 (MQ=255)
  aCTATTCTGGCTATTTACTTCGTGCCGCTGTTCTTTGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCggg                                                                                                                                                   <  1:1027733/150‑1 (MQ=255)
  aCTATTCTGGCTATTTACTTCGTGCCGCTGTTCTTTGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCggg                                                                                                                                                   <  2:491179/150‑1 (MQ=255)
      ttattgCTATTTACTTAGTTCATCTGTTCTTTGTGCTGGTTAGCCTCCGCTTCCCGCTTAAGCCGCGACATGAATAATCAATAAAAAATTCGACATGCAAATTTGTAGCCTTTTTCAACTTTCCGATTAAGAACATTATCAGCGGTTTCt                                                                                                                                               <  2:56498/145‑1 (MQ=255)
       tCTGGCTATTTACTTCGTGCCGCTGTTCTTTGTGCTGGTGCGCCGCCGCTTCCCTCTGAAGCCGCGCCCGTAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCTCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCtt                                                                                                                                              <  1:998965/150‑1 (MQ=255)
       tCTGGCTATTTACTTCGTGCCGCTGTTCTTTGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCCGCTGAGCAGGTTCtt                                                                                                                                              <  1:229863/150‑1 (MQ=255)
        cTGGCTATTTACTTCGTGCCGCTGTTCTTTGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTg                                                                                                                                             >  2:800649/1‑150 (MQ=255)
          ggCTATTTACTTCGTGCCGCTGTTCTTTGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCt                                                                                                                                           <  1:231817/150‑1 (MQ=255)
          ggCTATTTACTTCGTGCCGCTGTTCTTTGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCt                                                                                                                                           >  2:308423/1‑150 (MQ=255)
           gCTATTTACTTCGTGCCGCTGTTCTTTGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCCGCTGAGCAGGTTCttaatc                                                                                                                                          <  1:1155745/150‑5 (MQ=255)
            cTATTTACTTCGTGCCGCGGTTCTTTGTGCTGGTGCGCCGCCGCTGCCCGCGGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCCGCTGAGCAGGTTCttaatcg                                                                                                                                         <  2:715476/150‑6 (MQ=255)
             tATTTACTTCGTGCCGCTGTTCTTTGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGtt                                                                                                                                        >  1:519718/1‑150 (MQ=255)
             tATTTACTTCGTGCCGCTGTTCTTTGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGtt                                                                                                                                        >  1:293628/1‑150 (MQ=255)
             tATTTACTTCGTGCCGCTGTTCTTTGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGtt                                                                                                                                        >  2:1530869/1‑150 (MQ=255)
                ttACTTCGTGCCGCTGTTCTTTGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGt                                                                                                                                     <  1:1136107/150‑1 (MQ=255)
                 tACTTCGTGCCGCTGTTCTTTGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTa                                                                                                                                    <  1:1111851/150‑1 (MQ=255)
                    ttCGTGCCGCTGTTCTTTGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACtt                                                                                                                                 <  1:964800/150‑1 (MQ=255)
                        tGCCGCTGTTCTTTGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGtc                                                                                                                             >  1:122299/1‑150 (MQ=255)
                           cGCTGTTCTTTGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCa                                                                                                                          >  1:1505657/1‑150 (MQ=255)
                            gCTGTTCTTTGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAg                                                                                                                         >  1:1010700/1‑150 (MQ=255)
                            gCTGTTCTTTGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAg                                                                                                                         <  1:1511281/150‑1 (MQ=255)
                            gCTGTTCTTTGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCCGCTGAGCAGGTTCttaatcg                                                                                                                                         >  1:1228686/1‑129 (MQ=255)
                            gCTGTTCTTTGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCCGCTGAGCAGGTTCttaatcg                                                                                                                                         <  2:1228686/134‑6 (MQ=255)
                                  cTTTGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGttt                                                                                                                                       >  2:815195/1‑130 (MQ=255)
                                  cTTTGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGttt                                                                                                                                       <  1:815195/130‑1 (MQ=255)
                                  cTTTGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAAtta                                                                                                                   <  2:293628/150‑1 (MQ=255)
                                       tgCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAAttacttac                                                                                                              >  1:403680/1‑150 (MQ=255)
                                        gCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACg                                                                                                             <  2:1188747/150‑1 (MQ=255)
                                        gCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACg                                                                                                             >  1:1016895/1‑150 (MQ=255)
                                        gCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACg                                                                                                             >  2:685488/1‑150 (MQ=255)
                                         cTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTGTTGCTGTTTGTACTTTGTCGCAGGAATTACTTACGa                                                                                                            >  2:1449899/1‑150 (MQ=255)
                                           ggTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAg                                                                                                          <  1:1046905/150‑1 (MQ=255)
                                               cgccgccgcTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATa                                                                                                      <  2:1173243/150‑1 (MQ=255)
                                                gccgacgcTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAAAGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGCTTCTTGCTGTTTGAACTTTGTCTCACGAATTAATTACGAAGCaaaa                                                                                                     >  2:199866/1‑150 (MQ=255)
                                                   gccgcTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGtt                                                                                                                                        <  2:536307/112‑1 (MQ=255)
                                                   gccgcTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGtt                                                                                                                                        >  1:536307/1‑112 (MQ=255)
                                                    ccgcTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTc                                                                                                 <  1:374628/150‑1 (MQ=255)
                                                      gcTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGa                                                                                               >  1:1048872/1‑150 (MQ=255)
                                                       cTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCACCTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGAt                                                                                              <  2:880803/150‑1 (MQ=255)
                                                       cTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGAt                                                                                              <  2:1505657/150‑1 (MQ=255)
                                                         tCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATaa                                                                                            >  1:465308/1‑150 (MQ=255)
                                                          cccGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATaaa                                                                                           >  1:796559/1‑150 (MQ=255)
                                                           ccGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAg                                                                                          >  2:1224857/1‑150 (MQ=255)
                                                            cGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGt                                                                                         >  1:1125134/1‑150 (MQ=255)
                                                                 aaGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTc                                                                                    <  2:1048699/150‑1 (MQ=255)
                                                                  aGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTcc                                                                                   <  2:465308/150‑1 (MQ=255)
                                                                  aGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTcc                                                                                   <  2:30911/150‑1 (MQ=255)
                                                                  aGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTcc                                                                                   <  1:308423/150‑1 (MQ=255)
                                                                  aGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTcc                                                                                   <  1:800649/150‑1 (MQ=255)
                                                                   gCCGCGCCCGTAATAAGCAATAAAAAAGGCGACATGCCATTGTGGCGCCTTTTTCAACTTTCCGATTAAGAACCTTCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCa                                                                                  <  1:1469240/150‑1 (MQ=255)
                                                                   gCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCa                                                                                  >  1:1236298/1‑150 (MQ=255)
                                                                     cgcgCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGt                                                                                <  1:650983/150‑1 (MQ=255)
                                                                      ccgcCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGtt                                                                               >  1:431262/2‑150 (MQ=255)
                                                                        gcCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTcc                                                                             >  1:406786/1‑150 (MQ=255)
                                                                         cccGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTccc                                                                            >  1:255761/1‑150 (MQ=255)
                                                                          ccGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTcccc                                                                           <  1:685488/150‑1 (MQ=255)
                                                                           cGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCa                                                                          >  2:567394/1‑150 (MQ=255)
                                                                            ggAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAATGTCTTTCCAGTTCCCCAg                                                                         <  1:910593/150‑1 (MQ=255)
                                                                              aaTAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGtt                                                                       >  2:1474253/1‑150 (MQ=255)
                                                                               aTAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTc                                                                      <  1:1530869/150‑1 (MQ=255)
                                                                               aTAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTc                                                                      >  2:612645/1‑150 (MQ=255)
                                                                                tAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTccagttcaccagttca                                                                     >  1:946551/1‑150 (MQ=255)
                                                                                  aGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACg                                                                   >  1:470922/1‑150 (MQ=255)
                                                                                  aGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACg                                                                   >  1:259964/1‑150 (MQ=255)
                                                                                    cAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGtt                                                                 >  2:411990/1‑150 (MQ=255)
                                                                                    cAATAAAAAAGGCGAAATGCCAATGTGTCGACTTTTTAAACTTTAAGATTAAGAACCTGCTCAGAGGGTTATTGCTGTTTTTACTTTGTCTCAGGAATTACTTAAGAAGAATAACTTATATAAAGTCTTTCAAGTTACCAAGttaaattt                                                                 >  2:1195224/1‑144 (MQ=255)
                                                                                      aTAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCa                                                               <  2:1010700/150‑1 (MQ=255)
                                                                                        aaaaaaGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAAt                                                             >  1:541781/1‑150 (MQ=255)
                                                                                         aaaaaGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATc                                                            <  2:403680/150‑1 (MQ=255)
                                                                                          aaaaGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCa                                                           >  1:583310/1‑150 (MQ=255)
                                                                                             aGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATaa                                                        >  1:712958/1‑150 (MQ=255)
                                                                                               gCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATaaca                                                      <  1:567394/150‑1 (MQ=255)
                                                                                                 gACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATaacaac                                                    >  1:1124823/1‑150 (MQ=255)
                                                                                                  aCATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAAcc                                                   <  1:612645/150‑1 (MQ=255)
                                                                                                    aTGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACctc                                                 >  2:45096/1‑150 (MQ=255)
                                                                                                     tGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTCAAGTTTCTACTTTGTCTCAGCAATTACTTACGAACCATAACTTCGATAAAGTCTTTCAAGTTCCCCAGTTAACCTTCAATCATAACAAccaaa                                                >  2:552283/1‑147 (MQ=255)
                                                                                                       ccAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACctctct                                              <  2:29488/150‑1 (MQ=255)
                                                                                                        cAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCtt                                             <  1:78235/150‑1 (MQ=255)
                                                                                                              gtCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAAtt                                       >  1:701284/1‑150 (MQ=255)
                                                                                                                cGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTTCTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTAt                                     <  2:1498321/150‑1 (MQ=255)
                                                                                                                  ccTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATgg                                   >  1:722316/1‑150 (MQ=255)
                                                                                                                   cTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATggg                                  <  1:1504311/150‑1 (MQ=255)
                                                                                                                     ttttCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTa                                <  1:1474253/150‑1 (MQ=255)
                                                                                                                       ttCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTAtt                              >  1:1209669/1‑150 (MQ=255)
                                                                                                                        tCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTc                             >  1:1418228/1‑150 (MQ=255)
                                                                                                                         cAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCt                            >  2:507501/1‑150 (MQ=255)
                                                                                                                           aCTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTAc                          >  1:614503/1‑150 (MQ=255)
                                                                                                                              ttCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGa                       <  2:155090/150‑1 (MQ=255)
                                                                                                                                 cGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACACCGTCTCTTATAATTATGGGTATTCTACGGCAAc                    >  2:766816/1‑150 (MQ=255)
                                                                                                                                  cagcaaGAACCCGCTGAGCAGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACa                   >  1:891185/5‑150 (MQ=255)
                                                                                                                                      aaGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAatat               >  2:1336161/1‑150 (MQ=255)
                                                                                                                                       aGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAAtata              >  1:1226143/1‑150 (MQ=255)
                                                                                                                                       aGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAAtata              <  1:1224857/150‑1 (MQ=255)
                                                                                                                                        taaCCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATAc             <  2:662888/149‑1 (MQ=255)
                                                                                                                                           ccTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCgt          >  1:1374620/1‑150 (MQ=255)
                                                                                                                                             tGCTCAGCTGGTTCTTGCTGTTTGTAATTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCAAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACCATAAACCgttg        >  2:1298413/1‑148 (MQ=255)
                                                                                                                                               cTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCgtggtg      >  2:1372238/1‑150 (MQ=255)
                                                                                                                                               cTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCgtggtg      <  2:1124823/150‑1 (MQ=255)
                                                                                                                                                  aGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCGTGGTGAAg   >  1:318763/1‑150 (MQ=255)
                                                                                                                                                   gCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTACCCAGTTCACGTTCAATCATAACAACATCTCTTATAATTATGGGTATTCTACGGAAACAATATACCGTGGTGAAGc  >  1:605221/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CCACTATTCTGGCTATTTACTTCGTGCCGCTGTTCTTTGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCGTGGTGAAGC  >  NC_000913/2590626‑2590922

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: