Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 568,690 T→A 6.4% intergenic (+43/‑212) 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,6900TA 6.4% 444.3 / 12.1 147intergenic (+43/‑212)emrE/ybcKDLP12 prophage; multidrug resistance protein/DLP12 prophage; putative recombinase
Reads supporting (aligned to +/- strand):  ref base T (65/72);  new base A (6/4);  total (71/76)
Fisher's exact test for biased strand distribution p-value = 5.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

Only 100 of 170 total aligned reads displayed.

GGGGATTTTTCGGCCAACGGCTGGACCTGCCAGCCATTATAGGCATGATGTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTATTGTTGTATATTCAAGATAATGTGACCTTAATTGTCTTTTTAAATAAA  >  NC_000913/568541‑568838
                                                                                                                                                     |                                                                                                                                                    
ggggATTTTTCGGCCAACGGCTGGACCTGCCAGCCATTATAGGCATGATGTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCttatat                                                                                                                                                      >  1:395269/1‑150 (MQ=255)
ggggATTTTTCGGCCAACGGCTGGACCTGCCAGCCATTATAGGCATGATGTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCttatat                                                                                                                                                      >  2:1492328/1‑150 (MQ=255)
  ggATTTTTCGGCCAACGGCTGGACCTGCCAGCCATTATAGGCATGATGTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCttatattt                                                                                                                                                    <  1:434794/150‑1 (MQ=255)
   gATTTTTCGGCCTACGGCTGGACCTGCCATCCATTATAGGCATGATGTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAATCACACCACATTAAATTAATTTGTTTCTAAACGACTAAAATATGGAGGCTCttatattta                                                                                                                                                   <  1:2033710/150‑1 (MQ=255)
   gATTTTTCGGCCAACGGCTGGACCTGCCAGCCATTATATGCATGATGTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCttatattta                                                                                                                                                   <  2:1931904/150‑1 (MQ=255)
      ttttCGGCCAACGGCTGGACCTGCCAGCCATTATAGGCATGATGTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCttatatttatat                                                                                                                                                <  1:913029/150‑1 (MQ=255)
      ttttCGGCCAACGGCTGGACCTGCCAGCCATTATAGGCATGATGTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCttatatttatat                                                                                                                                                <  1:1262178/150‑1 (MQ=255)
         tCGGCCAACGGCTGGACCTGCCAGCCATTATAGGCATGATGTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAg                                                                                                                                             <  1:1239451/150‑1 (MQ=255)
          cGGCCAACGGCTGGACATGCCAGCCATTATAGGCATGATCTTGATTTGTGCCGGTGTGTTGATTAGTAATTTATTGTCACGAAGCACACCACAGTAAAATAATTTGTTTCTAAACGAGTAAAATATTGAGGCTCTTTTATTTATATGAgg                                                                                                                                            >  1:409213/1‑149 (MQ=255)
             ccAACGGCTGGACCTGCCAGCCATTATAGGCATGATGTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCATATAAATATAAGAGCCTc                                                                                                                                         <  1:733395/150‑1 (MQ=255)
              cAACGGCTGGACCTGCCATCCATTATAGTCATGATGTAGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCg                                                                                                                                        <  1:1394342/150‑1 (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)
                   gCTGGACCTGCCAGCCATTATAGGCATGATGGTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACAACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTa                                                                                                                                   <  2:1071561/150‑1 (MQ=255)
                     tGGACCTGCCAGCCATTATAGGCATGATGTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATg                                                                                                                                 <  1:2142275/150‑1 (MQ=255)
                           tGCCAGCCATTATAGGCATGATGTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCATCGTTTTATGCttttt                                                                                                                           >  1:244915/1‑150 (MQ=255)
                             ccagccaTTATAGGCATGATGGTGATTTGTGCCGGGGTGTTGATTATTAATTTATTTTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATCTGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGt                                                                                                                         <  2:338893/150‑1 (MQ=255)
                              cagccaTTATAGGCATGATGTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGATGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGtt                                                                                                                        >  2:210557/1‑150 (MQ=255)
                              cagccaTTATAGGCATGATGTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGtt                                                                                                                        <  2:234310/150‑1 (MQ=255)
                              cagccaTTATAGGCATGATGTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGtt                                                                                                                        >  1:1494845/1‑150 (MQ=255)
                              aaggcaTTATAGGCATGATGTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGtt                                                                                                                        <  2:1377628/146‑1 (MQ=255)
                                   aTTATAGGCATGATGTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGt                                                                                                                   >  1:1395789/1‑150 (MQ=255)
                                     tataGGCATGATGTTGATTTGTGCCTGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCt                                                                                                                 <  1:515515/150‑1 (MQ=255)
                                      ataGGCATGATGTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCATCGTTTTATGCTTTTTGTTAATGTCtt                                                                                                                >  1:1465613/1‑150 (MQ=255)
                                           cATGATGTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTAttt                                                                                                           >  1:1452696/1‑150 (MQ=255)
                                            atgatgTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTAtttt                                                                                                          >  1:1873854/1‑150 (MQ=255)
                                               atgTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTa                                                                                                       <  2:417320/150‑1 (MQ=255)
                                                tgTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTAt                                                                                                      >  2:1712940/1‑150 (MQ=255)
                                                tgTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTAt                                                                                                      <  2:928864/150‑1 (MQ=255)
                                                  ttGATTTTTGCCGGTGTGTTTATTATTAATTAAATGTCACGAAGCACAACACATTAAAATAATTTGTTTCTACACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGt                                                                                                    <  2:1705309/150‑1 (MQ=255)
                                                   tGATTTTTGCCGGTGTGTTGATTATTAATTTATTTTAACGAAGCACACCACATTAAAATAATTTTTTTATAAAACACGAAAATATGGAGGCTCGTATAATTATATGTGCCACCTTTTTTTCTTTTTTTTAATGTAttttatttttatttt                                                                                                   >  2:1633823/1‑147 (MQ=255)
                                                        tgtgCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCtt                                                                                              >  1:34343/1‑150 (MQ=255)
                                                            ccGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTtgt                                                                                          <  2:697623/150‑1 (MQ=255)
                                                             gggTGTGTTGATTATTAATGTATTGTAACGAAGCACACCACATTAATCTAATTTGTTTCTACACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTtgtg                                                                                         <  1:856476/149‑1 (MQ=255)
                                                                 gtgtTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCttt                                                                                     >  1:499875/1‑150 (MQ=255)
                                                                 gtgtTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCttt                                                                                     >  2:141427/1‑150 (MQ=255)
                                                                 gtgtTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCttt                                                                                     >  1:1502220/1‑150 (MQ=255)
                                                                  tgtTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTc                                                                                    <  2:1801406/150‑1 (MQ=255)
                                                                   gtTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCa                                                                                   <  2:1787852/150‑1 (MQ=255)
                                                                     tGATTATTAATTTATTGTCACGAAGCTCACCACATTAATATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGTGCCTAGTTTTATGATTTTTGTTAATGTCATTATTTTTTATGTATTCTTTTGTGCTTTCAAg                                                                                 >  2:317383/1‑150 (MQ=255)
                                                                     tGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCaat                                                                                 >  1:633168/1‑149 (MQ=255)
                                                                     tGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAg                                                                                 >  1:1045692/1‑150 (MQ=255)
                                                                     tGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAg                                                                                 <  2:1200357/150‑1 (MQ=255)
                                                                      gATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGa                                                                                >  1:855517/1‑150 (MQ=255)
                                                                         tattaATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGACTTTATTTTTTATGAATTCTTTTGTGCTTTCAAGATTa                                                                             >  2:218285/1‑150 (MQ=255)
                                                                            taATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATgg                                                                          <  1:1715330/150‑1 (MQ=255)
                                                                            taATTTATTGTCACGAAGCACACCACATGAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATgg                                                                          <  2:1759101/150‑1 (MQ=255)
                                                                             aaTTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGc                                                                         >  1:697287/1‑150 (MQ=255)
                                                                              aTTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCg                                                                        <  2:1883515/150‑1 (MQ=255)
                                                                               tttATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGt                                                                       <  1:2198281/150‑1 (MQ=255)
                                                                                     gTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTTCTTTCAAGATTATGGCGTAAGaaa                                                                 >  2:1561655/1‑150 (MQ=255)
                                                                                     gTCACGAAGCACACCACATCAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGaaa                                                                 >  2:1572205/1‑150 (MQ=255)
                                                                                        aCGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAAtt                                                              >  2:1850440/1‑150 (MQ=255)
                                                                                         ctaaGCACACAACATTAAAATAATTTGTTGCTAAAAGACTAAAATATGGAGGCTCTTATATTTATATGAGCCGCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTg                                                             <  2:1395789/148‑1 (MQ=255)
                                                                                          gAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGc                                                            <  2:1105628/150‑1 (MQ=255)
                                                                                            aGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCaa                                                          <  2:135082/150‑1 (MQ=255)
                                                                                                 accacATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGTTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGa                                                     <  2:898522/150‑1 (MQ=255)
                                                                                                  ccacATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGat                                                    >  2:625011/1‑150 (MQ=255)
                                                                                                     caTTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGattat                                                 >  2:648658/1‑150 (MQ=255)
                                                                                                     aaTTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGattat                                                 <  1:601119/149‑1 (MQ=255)
                                                                                                       ttAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTAttg                                               >  1:1595872/1‑150 (MQ=255)
                                                                                                         aaaaTAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTAttgtt                                             >  2:956181/1‑150 (MQ=255)
                                                                                                         aaaaTAAATTGTTTCTAAAAGTCTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATTCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTAttgtt                                             <  1:2024234/150‑1 (MQ=255)
                                                                                                           aataatTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTATtgttgt                                           >  2:1613491/1‑150 (MQ=255)
                                                                                                           aataatTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTATtgttgt                                           >  1:1137326/1‑150 (MQ=255)
                                                                                                           aataatTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTATtgttgt                                           <  2:1476207/150‑1 (MQ=255)
                                                                                                           aataatTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTATtgttgt                                           >  2:397244/1‑150 (MQ=255)
                                                                                                                 ttGTATCTAAACGCCTAAAATATGGAGGGTATTATATTTATATGAGCCTCGTTTTTTTCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTATTGTTGTATATTc                                     <  2:382238/150‑1 (MQ=255)
                                                                                                                  tGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTATTGTTGTATATTCa                                    <  1:1717756/150‑1 (MQ=255)
                                                                                                                     ttCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTTTTAAAGTCTTTATTTTTTATGTATTATTTTGAGCTTTCAAGATTATTTCTTAAGAAAATTGCAATACTATTATTGTTGTATATTCaata                                 >  1:361292/1‑148 (MQ=255)
                                                                                                                     ttCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTATTGTTGTATATTCAAGa                                 >  1:1399267/1‑150 (MQ=255)
                                                                                                                     ttCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTATTGTTGTATATTCAAGa                                 <  2:733163/150‑1 (MQ=255)
                                                                                                                      tCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTATTGTTGTATATTCAAGAt                                >  1:1572861/1‑150 (MQ=255)
                                                                                                                          aaCGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGGTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTTTGGCGTAAGAAAATTGCAATACGATTATTGTTGTATATTCAAGATAAtt                            >  1:245160/1‑149 (MQ=255)
                                                                                                                          aaCGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTATTGTTGTATATTCAAGATAAtg                            <  2:1494845/150‑1 (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)
                                                                                                                                 aaaaTATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTATTGTTGTATATTCAAGATAATGTGACCtt                     >  1:1986450/1‑150 (MQ=255)
                                                                                                                                  ataaaacgAGGCTCATATAAATATAAGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTATTGTTGTATATTCAAGATAATGTGACCTTa                    >  2:1376590/8‑150 (MQ=255)
                                                                                                                                    atatGGAGGCGCTTATAATTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTATTGTTGTATATTCAAGATAATGTGACCTTAAt                  <  2:346504/150‑1 (MQ=255)
                                                                                                                                     aaacgAGGCTCATATAAATATAAGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCATTTGTGGTTTCAAGATTATGGCGTAAGAAAATTTCAATACGATTATAGTTGTATATTCAAGATAATGAGACATTAAtt                 >  2:128606/5‑150 (MQ=255)
                                                                                                                                      atGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTATTGTTGTATATTCAAGATAATGTGACCTTAATTg                <  2:1030878/150‑1 (MQ=255)
                                                                                                                                      atGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTATTGTTGTATATTCAAGATAATGTGACCTTAATTg                <  2:1021810/150‑1 (MQ=255)
                                                                                                                                        cgAGGCTCATATAAATATAAGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTATTGTTGTATATTCAAGATAATGTGACCTTAATTGTc              >  2:1441413/2‑150 (MQ=255)
                                                                                                                                           ggCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTATTGTTGTATATTCAAGATAATGTGACCTTAATTGTCttt           >  2:1412212/1‑150 (MQ=255)
                                                                                                                                             ctctTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTATTGTTGTATATTCAAGATAATGTGACCTTAATTGTCttttt         <  2:893206/150‑1 (MQ=255)
                                                                                                                                              tctTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTATTGTTGTATATTCAAGATAATGTGACCTTAATTGTCTTTTta        <  2:1502220/150‑1 (MQ=255)
                                                                                                                                                ttatatttatatGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTATTGTTGTATATTCAAGATAATGTGACCTTAATTGTCTTTTtaaa      <  1:2054419/150‑1 (MQ=255)
                                                                                                                                                ttaaatttatatGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATTTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTATTGTTGTATATTCAAGATAATGTGACCTTAATTGTCTTTTtaaa      <  1:487190/150‑1 (MQ=255)
                                                                                                                                                 tatatttatatGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTATTGTTGTATATTCAAGATAATGTGACCTTAATTGTCTTTTtaaat     >  2:356029/1‑150 (MQ=255)
                                                                                                                                                 tatatttatatGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTATTGTTGTATATTCAAGATAATGTGACCTTAATTGTCTTTTtaaat     <  1:1065605/150‑1 (MQ=255)
                                                                                                                                                   tatttatatGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTATTGTTGTATATTCAAGATAATGTGACCTTAATTGTCTTTTtaaataa   <  2:653647/150‑1 (MQ=255)
                                                                                                                                                   tatttatatGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTATTGTTGTATATTCAAGATAATGTGACCTTAATTGTCTTTTtaaataa   >  2:699244/1‑150 (MQ=255)
                                                                                                                                                    atttatatGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTATTGTTGTATATTCAAGATAATGTGACCTTAATTGTCTTTTtaaataaa  >  2:31267/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GGGGATTTTTCGGCCAACGGCTGGACCTGCCAGCCATTATAGGCATGATGTTGATTTGTGCCGGTGTGTTGATTATTAATTTATTGTCACGAAGCACACCACATTAAAATAATTTGTTTCTAAACGACTAAAATATGGAGGCTCTTATATTTATATGAGCCTCGTTTTATGCTTTTTGTTAATGTCTTTATTTTTTATGTATTCTTTTGTGCTTTCAAGATTATGGCGTAAGAAAATTGCAATACGATTATTGTTGTATATTCAAGATAATGTGACCTTAATTGTCTTTTTAAATAAA  >  NC_000913/568541‑568838

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: