Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,590,767 T→C 6.8% intergenic (+59/+40) 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,7670TC 6.8% 443.6 / 14.2 146intergenic (+59/+40)acrD/ypfMmultidrug efflux pump RND permease AcrD/uncharacterized protein YpfM
Reads supporting (aligned to +/- strand):  ref base T (65/71);  new base C (2/8);  total (67/79)
Fisher's exact test for biased strand distribution p-value = 1.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 149 total aligned reads displayed.

TTTCGGCCACTATTCTGGCTATTTACTTCGTGCCGCTGTTCTTTGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCGTGG  >  NC_000913/2590620‑2590916
                                                                                                                                                   |                                                                                                                                                     
tttCGGCCACTATTCTGGCTATTTACTTCGTGCCGCTGTTCTTTGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGc                                                                                                                                                     <  1:1112109/150‑1 (MQ=255)
 ttCGGCCACTATTCTGGCTATTTACTTCGTGCCGCTGTTCTTTGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCt                                                                                                                                                    <  1:1694488/150‑1 (MQ=255)
  tCGGCCACTATTCTGGCTATTTACTTCGTGCCGCTGTTCTTTGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTc                                                                                                                                                   <  1:868535/150‑1 (MQ=255)
  tCGGCCACTATTCTGGCTATTTACTTCGTGCCGCTGTTCTTTGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCCGCtg                                                                                                                                                   >  1:1756154/1‑149 (MQ=255)
      ccACTATTCTGGCTATTTACTTCGTGCCGCTGTTCTTTGTGCTGGTGCGCCGCCGCTTCCCGCTGGAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCg                                                                                                                                               >  2:115762/1‑150 (MQ=255)
       cACTATTCTGGCTATTTACTTCGTGCCGCTGTTCTTTGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCgg                                                                                                                                              <  2:269950/150‑1 (MQ=255)
        aCTATTCTGGCTATTTACTTCGTGCCGCTGTTCTTTGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCggg                                                                                                                                             <  1:910919/150‑1 (MQ=255)
        aCTATTCTGGCTATTTACTTCGTGCCGCTGTTCTTTGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCggg                                                                                                                                             <  2:1181672/150‑1 (MQ=255)
          tATTCTGGCTATTTACTTCGTGCCGCTGTTCTTTGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGtt                                                                                                                                           <  2:686325/150‑1 (MQ=255)
          tATTCTGGCTATTTACTTCGTGCCGCTGTTCTTTGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGtt                                                                                                                                           <  1:691973/150‑1 (MQ=255)
            ttCTGGCTATTTACTTCGTGCCGCTGTTCTTTGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCt                                                                                                                                         <  2:558223/150‑1 (MQ=255)
            ttCTGGCTATTTACTTCGTGCCGCTGTTCTTTGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCt                                                                                                                                         <  1:1877424/150‑1 (MQ=255)
             tCTGGCTATTTACTTCGTGCCGCTGTTCTTTGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCtt                                                                                                                                        <  1:1825165/150‑1 (MQ=255)
             tCTGGCTATTTACTTCGTGCCGCTGTTCTTTGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCCGCTGAGCAGGTTCtt                                                                                                                                        <  1:2241975/150‑1 (MQ=255)
             tCTGGCTATTTACTTCGTGCCGCTGTTCTTTGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCCGCTGAGCAGGTTCtt                                                                                                                                        <  1:294743/150‑1 (MQ=255)
              cTGGCTATTTACTTCGTGCCGCTGTTCTTTGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCCGCTGAGCAGGTTCtta                                                                                                                                       <  1:1808106/150‑2 (MQ=255)
               tGGCTATTTACTTCGTGCCGCTGTTCTTTGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGc                                                                                                                                      <  2:939728/150‑1 (MQ=255)
               tGGCTATTTACTTCGTGCCGCTGTTCTTTGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGc                                                                                                                                      >  1:1441567/1‑150 (MQ=255)
                 gCTATTTACTTCGTGCCGCTGTTCTTTGTTCTGGTTCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCATTAAAAAAGTCTACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCCGCTGAGCAGGTTCttaatc                                                                                                                                    <  2:423052/150‑5 (MQ=255)
                 gCTATTTACTTCGTGCCGCTGTTCTTTGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCCGCTGAGCAGGTTCttaatc                                                                                                                                    <  1:1518541/150‑5 (MQ=255)
                     tttACTTCGTGCCGCTGTTCTTTGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTg                                                                                                                                <  1:1173862/150‑1 (MQ=255)
                          ttCGTGCCGCTGTTCTTTGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACtt                                                                                                                           <  1:1047160/150‑1 (MQ=255)
                           tCGTGCCGCTGTTCTTTGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACttt                                                                                                                          >  2:810089/1‑150 (MQ=255)
                            cGTGCCGCTGTTCTTTGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTg                                                                                                                         >  2:75952/1‑150 (MQ=255)
                               gCCGCTGTTCTTTGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGtct                                                                                                                      <  1:150618/150‑1 (MQ=255)
                                ccGCTGTTCTTTGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGtctc                                                                                                                     >  1:1057520/1‑150 (MQ=255)
                                ccGCTGTTCTTTGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGtctc                                                                                                                     <  1:895816/150‑1 (MQ=255)
                                 cGCTGTTCTTTGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCggg                                                                                                                                             <  2:1793005/125‑1 (MQ=255)
                                  gCTGTTCTTTGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAg                                                                                                                   <  2:1356077/150‑1 (MQ=255)
                                   cTGTTCTTTGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAgg                                                                                                                  <  2:2009395/150‑1 (MQ=255)
                                   cTGTTCTTTGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAgg                                                                                                                  <  1:1865356/150‑1 (MQ=255)
                                    tGTTCTTTGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGa                                                                                                                 >  2:1671858/1‑150 (MQ=255)
                                     gTTCTTTGTGCTGGTGCGGCGCGGCTTCCCGCTGATGCCGCGGCCGGTATTAGCAATTAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTGAGCGGGTTCTTGCTGTTTGTATTTTGTCTCAGGaa                                                                                                                <  2:1013446/150‑1 (MQ=255)
                                     gTTCTTTGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGaa                                                                                                                >  1:1839740/1‑150 (MQ=255)
                                     gTTCTTTGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGaa                                                                                                                <  2:1908201/150‑1 (MQ=255)
                                      ttCTTTGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAAt                                                                                                               <  1:1036011/150‑1 (MQ=255)
                                        cTTTGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAAtta                                                                                                             >  2:1655357/1‑150 (MQ=255)
                                        cTTTGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAAtta                                                                                                             >  2:981778/1‑150 (MQ=255)
                                         tttGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAAttac                                                                                                            <  2:1149149/150‑1 (MQ=255)
                                          ttGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAAttact                                                                                                           <  1:1636110/150‑1 (MQ=255)
                                          ttGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAAttact                                                                                                           <  2:488443/150‑1 (MQ=255)
                                           tgtgCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAAttactt                                                                                                          >  1:1473066/1‑150 (MQ=255)
                                              gCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACg                                                                                                       >  1:1037431/1‑150 (MQ=255)
                                                 ggTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAg                                                                                                    >  1:225714/1‑150 (MQ=255)
                                                  gTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGc                                                                                                   >  2:2325565/1‑150 (MQ=255)
                                                     cgccgccgcTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATa                                                                                                >  1:5455/1‑150 (MQ=255)
                                                         gccgcTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACtt                                                                                            >  1:1211876/1‑150 (MQ=255)
                                                           cgcTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCg                                                                                          <  2:1562602/150‑1 (MQ=255)
                                                            gcTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGa                                                                                         >  1:321231/1‑150 (MQ=255)
                                                            gcTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGa                                                                                         <  2:635500/150‑1 (MQ=255)
                                                              ttCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATa                                                                                       >  1:180909/1‑150 (MQ=255)
                                                               tCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATaa                                                                                      >  2:1584782/1‑150 (MQ=255)
                                                                cccGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATaaa                                                                                     >  2:1026622/1‑150 (MQ=255)
                                                                  cGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGt                                                                                   <  1:1123179/150‑1 (MQ=255)
                                                                  cGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGt                                                                                   >  1:500261/1‑150 (MQ=255)
                                                                     tGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGCCATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCtt                                                                                >  2:1983341/1‑150 (MQ=255)
                                                                      gAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCATTGTGGCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCttt                                                                               <  2:1275068/150‑1 (MQ=255)
                                                                        aGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTcc                                                                             <  2:1131712/150‑1 (MQ=255)
                                                                        aGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTcc                                                                             >  2:2360117/1‑150 (MQ=255)
                                                                        aGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTcc                                                                             <  2:1890928/150‑1 (MQ=255)
                                                                         gCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCa                                                                            >  2:30381/1‑150 (MQ=255)
                                                                         gCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCa                                                                            >  1:2455932/1‑150 (MQ=255)
                                                                          ccGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAg                                                                           <  1:2072916/150‑1 (MQ=255)
                                                                              gcCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGCAGCATACCTTCGATAAAGTCTTTCCAGTTcc                                                                       >  2:1208748/1‑150 (MQ=255)
                                                                              gcCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTcc                                                                       <  2:2130081/150‑1 (MQ=255)
                                                                                 cGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCa                                                                    >  1:796280/1‑150 (MQ=255)
                                                                                        aGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACg                                                             >  1:2130593/1‑150 (MQ=255)
                                                                                         gCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAACGTCTTGCCAGTTCCCCAGTTCACGt                                                            >  2:1134462/1‑150 (MQ=255)
                                                                                           aaTAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTc                                                          <  1:30381/150‑1 (MQ=255)
                                                                                                aacaGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCATTAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTTCCCAGTTCTAGTTCAATCa                                                     <  2:1641613/147‑1 (MQ=255)
                                                                                                 aaaGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTATCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATctt                                                    >  1:1569279/1‑148 (MQ=255)
                                                                                                  aaGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGTGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATa                                                   <  2:2388480/150‑1 (MQ=255)
                                                                                                  aaGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATa                                                   <  2:1319507/150‑1 (MQ=255)
                                                                                                  aaGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTAAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATa                                                   >  2:2263445/1‑150 (MQ=255)
                                                                                                    ggCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATaac                                                 >  1:1552513/1‑150 (MQ=255)
                                                                                                     gCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATaaca                                                <  2:1473066/150‑1 (MQ=255)
                                                                                                               aaTGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTta                                      <  1:1697153/150‑1 (MQ=255)
                                                                                                                 tgtgTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTtata                                    <  2:2092362/150‑1 (MQ=255)
                                                                                                                   tgtCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAAt                                  <  2:180909/150‑1 (MQ=255)
                                                                                                                   tgtCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAAt                                  <  1:814650/150‑1 (MQ=255)
                                                                                                                    gtCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAAtt                                 >  1:590338/1‑150 (MQ=255)
                                                                                                                      cGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTAt                               >  2:5641/1‑150 (MQ=255)
                                                                                                                       gCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATg                              <  2:1738929/150‑1 (MQ=255)
                                                                                                                            tttCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTAt                         <  2:225714/150‑1 (MQ=255)
                                                                                                                            tttCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTAt                         >  2:336153/1‑150 (MQ=255)
                                                                                                                             ttCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTAtt                        >  1:617534/1‑150 (MQ=255)
                                                                                                                              tCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTc                       <  2:321231/150‑1 (MQ=255)
                                                                                                                               cAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCt                      >  2:1407044/1‑150 (MQ=255)
                                                                                                                                 aCTTTCCGATTAAGAACCTGCTCATCTGGTTCTTGCTGTTTGTACTTTGTCTCTGGAATTACTTACGTAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAAGTACGTCTCTTATAATTATGGGTATTCTAc                    <  2:1338795/150‑1 (MQ=255)
                                                                                                                                  cTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACg                   >  1:1548132/1‑150 (MQ=255)
                                                                                                                                  cTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACg                   >  1:1486963/1‑150 (MQ=255)
                                                                                                                                  cTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACg                   >  1:868316/1‑150 (MQ=255)
                                                                                                                                       cGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAAc              >  1:1758848/1‑150 (MQ=255)
                                                                                                                                       cGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAAc              <  1:1627320/150‑1 (MQ=255)
                                                                                                                                            aaGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTACCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTAAGGAAAAAatat         >  1:100734/1‑150 (MQ=255)
                                                                                                                                             aGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAAtata        >  2:1029661/1‑150 (MQ=255)
                                                                                                                                             aGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAAtata        >  1:2328204/1‑150 (MQ=255)
                                                                                                                                               aaCCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATAcc      >  2:2395921/1‑150 (MQ=255)
                                                                                                                                               aaCCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATAcc      >  1:628536/1‑150 (MQ=255)
                                                                                                                                                   tGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCgtgg  <  1:2325565/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TTTCGGCCACTATTCTGGCTATTTACTTCGTGCCGCTGTTCTTTGTGCTGGTGCGCCGCCGCTTCCCGCTGAAGCCGCGCCCGGAATAAGCAATAAAAAAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCGTGG  >  NC_000913/2590620‑2590916

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: