Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 568,691 T→A 5.4% intergenic (+44/‑211) emrE → / → ybcK DLP12 prophage; multidrug resistance protein/DLP12 prophage; putative recombinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913568,6910TA 5.4% 435.9 / 11.1 147intergenic (+44/‑211)emrE/ybcKDLP12 prophage; multidrug resistance protein/DLP12 prophage; putative recombinase
Reads supporting (aligned to +/- strand):  ref base T (66/73);  new base A (5/3);  total (71/76)
Fisher's exact test for biased strand distribution p-value = 4.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.20e-01

Only 100 of 169 total aligned reads displayed.

GGGATTTTTCGGCCAACGGCTGGACCTGCCAGCCATTATAGGCATGATGTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTATTGTTGTATATTCAAGATAATGTGACCTTAATTGTCTTTTTAAATAAAAA  >  NC_000913/568542‑568840
                                                                                                                                                     |                                                                                                                                                     
gggATTTTTCGGCCAACGGCTGGACCTGCCAGCCATTATAGGCATGATGTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCttatatt                                                                                                                                                       >  1:667680/1‑150 (MQ=255)
  gATTTTTCGGCCTACGGCTGGACCTGCCATCCATTATAGGCATGATGTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAATCACACCACATTAAATTAATTTGTTTCTAAACGACTAAAATATGGAGGCTCttatattta                                                                                                                                                     <  1:2033710/150‑1 (MQ=255)
  gATTTTTCGGCCAACGGCTGGACCTGCCAGCCATTATATGCATGATGTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCttatattta                                                                                                                                                     <  2:1931904/150‑1 (MQ=255)
  gATTTTTCGGCCAACGGCTGGACCTGCCAGCCATTATAGGCATGATGTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCttatattta                                                                                                                                                     >  1:959717/1‑150 (MQ=255)
  gATTTTTCGGCCAACGGCTGGACCTGCCAGCCATTATAGGCATGATGTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCttatattta                                                                                                                                                     <  1:43269/150‑1 (MQ=255)
     ttttCGGCCAACGGCTGGACCTGCCAGCCATTATAGGCATGATGTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCttatatttatat                                                                                                                                                  <  1:1262178/150‑1 (MQ=255)
      tttCGGCCAACGGCTGGACCTGCCAGCCATTATAGGCATGATGTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATg                                                                                                                                                 <  2:1681021/150‑1 (MQ=255)
        tCGGCCAACGGCTGGACCTGCCAGCCATTATAGGCATGATGTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAg                                                                                                                                               <  1:1239451/150‑1 (MQ=255)
         cGGCCAACGGCTGGACCTGCCAGCCATTATAGGCATGATGTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGc                                                                                                                                              >  2:1784949/1‑150 (MQ=255)
             cAACGGCTGGACCTGCCATCCATTATAGTCATGATGTAGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCg                                                                                                                                          <  1:1394342/150‑1 (MQ=255)
             cAACGGCTGGACCTGCCAGCCATTATAGGCATGATGTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCg                                                                                                                                          >  1:776937/1‑150 (MQ=255)
             cAACGGCTGGACCTGCCAGCCATTATAGGCATGATGTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCg                                                                                                                                          >  1:1489350/1‑150 (MQ=255)
              aaCGGCTGGACCTGCCAGCCATTATAGGCATGATGTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGt                                                                                                                                         >  2:604441/1‑150 (MQ=255)
               aCGGCTGGACCTGCCAGCCATTATAGGCATGATGTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCATATAAATATAAGAGCCTccat                                                                                                                                        <  2:542160/150‑4 (MQ=255)
                cGGCTGGACCTGCCAGCCATTATAGGCATGATGTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGttt                                                                                                                                       <  1:1908038/150‑1 (MQ=255)
                 ggCTGGACCTGCCAGCCATTATAGGCATGATGTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGtttt                                                                                                                                      >  1:782804/1‑150 (MQ=255)
                    tGGACCTGCCAGCCATTATAGGCATGATGTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATg                                                                                                                                   <  1:2142275/150‑1 (MQ=255)
                      gACCTGCCAGCCATTATAGGCATGATGTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGTCTCTTATATTTATATGATCATCTTTTTATgat                                                                                                                                 >  1:382238/1‑148 (MQ=255)
                         cTGCCAGCCATTATAGGCATGATGTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCtttt                                                                                                                              >  2:149995/1‑150 (MQ=255)
                           gccagccaTTATAGGCATGATGTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGTTTTTTg                                                                                                                            >  1:898522/1‑150 (MQ=255)
                            ccagccaTTATAGGCATGATGGTGATTTGTGCCGGGGTGTTGATTATTAATTTATTTTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATCTGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGt                                                                                                                           <  2:338893/150‑1 (MQ=255)
                             cagccaTTATAGGCATGATGTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGtt                                                                                                                          >  1:1494845/1‑150 (MQ=255)
                             cagccaTTATAGGCATGATGTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGtt                                                                                                                          <  2:234310/150‑1 (MQ=255)
                             aaggcaTTATAGGCATGATGTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGtt                                                                                                                          <  2:1377628/146‑1 (MQ=255)
                                  aTTATAGGCATGATGTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGt                                                                                                                     >  1:1395789/1‑150 (MQ=255)
                                          cATGATGTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTAttt                                                                                                             >  1:1452696/1‑150 (MQ=255)
                                           atgatgTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTAtttt                                                                                                            >  1:1873854/1‑150 (MQ=255)
                                               tgTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTAt                                                                                                        >  2:1712940/1‑150 (MQ=255)
                                               tgTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTAt                                                                                                        <  2:928864/150‑1 (MQ=255)
                                                gTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATg                                                                                                       >  2:1374100/1‑150 (MQ=255)
                                                 ttGATTTTTGCCGGTGTGTTTATTATTAATTAAATGTCACGAAGCACAACACATTAAAATAATTTGTTTCTACACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGt                                                                                                      <  2:1705309/150‑1 (MQ=255)
                                                       tgtgCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCtt                                                                                                >  1:34343/1‑150 (MQ=255)
                                                       tgtgCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCtt                                                                                                >  2:131392/1‑150 (MQ=255)
                                                           ccGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTtgt                                                                                            <  2:697623/150‑1 (MQ=255)
                                                           ccGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTtgt                                                                                            >  1:1951260/1‑150 (MQ=255)
                                                                gtgtTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCttt                                                                                       >  2:141427/1‑150 (MQ=255)
                                                                 tgtTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTc                                                                                      <  2:1801406/150‑1 (MQ=255)
                                                                 ggtTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATATGTTTCTAAACGACTAAAATCTGGAGGCGCTTATATTTATATGAGCCGCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTc                                                                                      <  1:1453834/149‑1 (MQ=255)
                                                                  gtTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCa                                                                                     <  2:1787852/150‑1 (MQ=255)
                                                                    tGATTATTAATTTATTGTCACGAAGCTCACCACATTAATATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGTGCCTAGTTTTATGATTTTTGTTAATGTCATTATTTTTTATGTATTCTTTTGTGCTTTCAAg                                                                                   >  2:317383/1‑150 (MQ=255)
                                                                    tGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAg                                                                                   >  1:1045692/1‑150 (MQ=255)
                                                                    tGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAg                                                                                   <  2:1200357/150‑1 (MQ=255)
                                                                     gATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGa                                                                                  >  1:855517/1‑150 (MQ=255)
                                                                     gATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGa                                                                                  <  2:1401413/150‑1 (MQ=255)
                                                                        tattaATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGACTTTATTTTTTATGAATTCTTTTGTGCTTTCAAGATTa                                                                               >  2:218285/1‑150 (MQ=255)
                                                                           taATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATgg                                                                            <  1:1715330/150‑1 (MQ=255)
                                                                            aaTTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGc                                                                           >  1:697287/1‑150 (MQ=255)
                                                                              tttATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGt                                                                         <  1:2198281/150‑1 (MQ=255)
                                                                                    gTCACGAAGCACACCACATCAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGaaa                                                                   >  2:1572205/1‑150 (MQ=255)
                                                                                    gTCACGAAGCACACCACATCAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGaaa                                                                   >  2:1564444/1‑150 (MQ=255)
                                                                                       aCGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAAtt                                                                >  2:1850440/1‑150 (MQ=255)
                                                                                        cGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGACTCGTTTTATGCTTTTTTTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTg                                                               <  2:1388481/150‑1 (MQ=255)
                                                                                         gAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGc                                                              <  2:1105628/150‑1 (MQ=255)
                                                                                              acacCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGAAATAc                                                         >  2:1924538/1‑150 (MQ=255)
                                                                                                accacATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGTTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGa                                                       <  2:898522/150‑1 (MQ=255)
                                                                                                 ccacATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGat                                                      >  2:625011/1‑150 (MQ=255)
                                                                                                    caTTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGattat                                                   >  2:648658/1‑150 (MQ=255)
                                                                                                      ttAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTAttg                                                 >  1:1595872/1‑150 (MQ=255)
                                                                                                       tAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTAttgt                                                <  2:2137204/150‑1 (MQ=255)
                                                                                                       tAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTAttgt                                                <  2:1489350/150‑1 (MQ=255)
                                                                                                        aaaaTAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTAttgtt                                               >  2:956181/1‑150 (MQ=255)
                                                                                                        aaaaTAAATTGTTTCTAAAAGTCTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATTCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTAttgtt                                               <  1:2024234/150‑1 (MQ=255)
                                                                                                         aaaTATTTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTAttgttg                                              <  1:459777/150‑1 (MQ=255)
                                                                                                          aataatTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTATtgttgt                                             >  2:1613491/1‑150 (MQ=255)
                                                                                                          aataatTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTATtgttgt                                             >  2:397244/1‑150 (MQ=255)
                                                                                                          aataatTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTATtgttgt                                             >  1:1137326/1‑150 (MQ=255)
                                                                                                           ataatTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTATTGTTGta                                            <  1:2199951/150‑1 (MQ=255)
                                                                                                           ataatTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTAAATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTATTGTTGta                                            <  2:448152/150‑1 (MQ=255)
                                                                                                             aatTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTATTGTTGtata                                          <  1:149995/150‑1 (MQ=255)
                                                                                                                ttGTATCTAAACGCCTAAAATATGGAGGGTATTATATTTATATGAGCCTCGTTTTTTTCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTATTGTTGTATATTc                                       <  2:382238/150‑1 (MQ=255)
                                                                                                                 tGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTATTGTTGTATATTCa                                      <  1:1717756/150‑1 (MQ=255)
                                                                                                                    ttCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTTTTAAAGTCTTTATTTTTTATGTATTATTTTGAGCTTTCAAGATTATTTCTTAAGAAAATTGCAATACTATTATTGTTGTATATTCaata                                   >  1:361292/1‑148 (MQ=255)
                                                                                                                    ttCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTATTGTTGTATATTCAAGa                                   <  2:733163/150‑1 (MQ=255)
                                                                                                                    ttCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTATTGTTGTATATTCAAGa                                   >  1:1399267/1‑150 (MQ=255)
                                                                                                                    ttCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTATTGTTGTATATTCAAGa                                   >  1:2043504/1‑150 (MQ=255)
                                                                                                                     tCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTATTGTTGTATATTCAAGAt                                  >  1:1572861/1‑150 (MQ=255)
                                                                                                                         aaCGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGGTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTTTGGCGTAAGAAAATTGCAATACGATTATTGTTGTATATTCAAGATAAtt                              >  1:245160/1‑149 (MQ=255)
                                                                                                                          aCGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTATTGTTGTATATTCAAGATAAtgt                             <  1:583796/150‑1 (MQ=255)
                                                                                                                          aCGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTATTGTTGTATATTCAAGATAAtgt                             <  1:70713/150‑1 (MQ=255)
                                                                                                                             aCTAAAATATGGAGGCTCTTATATTTATATGAGCGTCGTTTTATGGTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACTATTTTTTTTGTATATTCAAGATAAtgtctc                          >  2:1523184/1‑147 (MQ=255)
                                                                                                                             aCTAAAATATGGAGGCTCTTATATTTAAATGAGCATCTTTTTAAGCTTTTTGTTACTGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACTATTATTGTTGTATATTCAAGATAATGTGAc                          <  1:1434579/150‑1 (MQ=255)
                                                                                                                              cTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTAGTTAATCTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTATTGTTGTATATTCAAGATAATGTGAca                         >  2:2018960/1‑149 (MQ=255)
                                                                                                                                   atatGGAGGCGCTTATAATTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTATTGTTGTATATTCAAGATAATGTGACCTTAAt                    <  2:346504/150‑1 (MQ=255)
                                                                                                                                    aaacgAGGCTCATATAAATATAAGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCATTTGTGGTTTCAAGATTATGGCGTAAGAAAATTTCAATACGATTATAGTTGTATATTCAAGATAATGAGACATTAAtt                   >  2:128606/5‑150 (MQ=255)
                                                                                                                                     atGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTATTGTTGTATATTCAAGATAATGTGACCTTAATTg                  <  2:1021810/150‑1 (MQ=255)
                                                                                                                                       cgAGGCTCATATAAATATAAGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTATTGTTGTATATTCAAGATAATGTGACCTTAATTGTc                >  2:1441414/2‑150 (MQ=255)
                                                                                                                                       cgAGGCTCATATAAATATAAGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTATTGTTGTATATTCAAGATAATGTGACCTTAATTGTc                >  2:1441413/2‑150 (MQ=255)
                                                                                                                                            ctctTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTATTGTTGTATATTCAAGATAATGTGACCTTAATTGTCttttt           <  2:893206/150‑1 (MQ=255)
                                                                                                                                             tctTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGTCGTAAGAAAATTGCAATACGATTATTGTTGTATATTCAAGATAATGTGACCTTAATTGTCTTTTta          <  1:562477/150‑1 (MQ=255)
                                                                                                                                             tctTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTATTGTTGTATATTCAAGATAATGTGACCTTAATTGTCTTTTta          <  1:604441/150‑1 (MQ=255)
                                                                                                                                             tctTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTATTGTTGTATATTCAAGATAATGTGACCTTAATTGTCTTTTta          <  2:1502220/150‑1 (MQ=255)
                                                                                                                                              ctTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTATTGTTGTATATTCAAGATAATGTGACCTTAATTGTCTTTTtaa         <  1:509674/150‑1 (MQ=255)
                                                                                                                                               ttatatttatatGAGCCTCGTTGTATGCTTTTTGTTAATGTCTATATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTATTGTTGTATATTCAAGATAATGTGACCTTAATTGTCTTTTtaaa        <  2:1049186/150‑1 (MQ=255)
                                                                                                                                               ttaaatttatatGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATTTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTATTGTTGTATATTCAAGATAATGTGACCTTAATTGTCTTTTtaaa        <  1:487190/150‑1 (MQ=255)
                                                                                                                                                tatatttatatGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTATTGTTGTATATTCAAGATAATGTGACCTTAATTGTCTTTTtaaat       <  1:1065605/150‑1 (MQ=255)
                                                                                                                                                tatatttatatGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTATTGTTGTATATTCAAGATAATGTGACCTTAATTGTCTTTTtaaat       >  2:356029/1‑150 (MQ=255)
                                                                                                                                                  tatttatatGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTATTGTTGTATATTCAAGATAATGTGACCTTAATTGTCTTTTtaaataa     >  2:699244/1‑150 (MQ=255)
                                                                                                                                                  tatttatatGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTATTGTTGTATATTCAAGATAATGTGACCTTAATTGTCTTTTtaaataa     <  2:653647/150‑1 (MQ=255)
                                                                                                                                                  tatttatatGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTATTGTTGTATATTCAAGATAATGTGACCTTAATTGTCTTTTtaaataa     >  2:1079315/1‑150 (MQ=255)
                                                                                                                                                   atttatatGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTATTGTTGTATATTCAAGATAATGTGACCTTAATTGTCTTTTtaaataaa    >  2:31267/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGGATTTTTCGGCCAACGGCTGGACCTGCCAGCCATTATAGGCATGATGTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTATTGTTGTATATTCAAGATAATGTGACCTTAATTGTCTTTTTAAATAAAAA  >  NC_000913/568542‑568840

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