Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,705,984 T→G 33.7% intergenic (‑194/+114) dppB ← / ← dppA dipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,705,9840TG33.7% 492.6 / inf 291intergenic (‑194/+114)dppB/dppAdipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base T (82/111);  new base G (92/6);  total (174/117)
Fisher's exact test for biased strand distribution p-value = 2.63e-19
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.34e-02

Only 100 of 295 total aligned reads displayed.

TTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCT  >  NC_000913/3705885‑3706083
                                                                                                   |                                                                                                           
  gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                                           <  1:1774437/100‑1 (MQ=255)
   ctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTAt                                                                                                          >  1:2797547/1‑100 (MQ=255)
    tgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATc                                                                                                         >  2:1239567/1‑100 (MQ=255)
     gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                                        >  2:1896680/1‑100 (MQ=255)
      ccGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCg                                                                                                       >  2:634613/1‑100 (MQ=255)
      caggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCg                                                                                                       <  1:2549909/98‑1 (MQ=255)
      caggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCg                                                                                                       <  1:3440299/98‑1 (MQ=255)
        ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGc                                                                                                     <  1:1981352/100‑1 (MQ=37)
        ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGc                                                                                                     <  2:1005494/100‑1 (MQ=37)
        ggATGGCGCTGCACGCAGACGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTGGTTTATGCCGTTTGCGGCTTTAACGCCtttttcg                                                                                                       >  2:1585974/1‑93 (MQ=25)
         gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                                    <  1:552597/100‑1 (MQ=255)
         gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTGTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                                    <  2:295458/100‑1 (MQ=255)
          aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                                    <  2:3391859/99‑1 (MQ=255)
          aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCt                                                                                                   <  2:1868115/100‑1 (MQ=255)
           tGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑a                                                                                                >  1:2932017/1‑100 (MQ=255)
              cgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑aa                                                                                            >  1:3822560/1‑100 (MQ=255)
                cTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATACCGGATCCGGCGTTAACGCCTTATCCGCCCT‑‑AC‑aaa                                                                                           >  2:2724628/1‑99 (MQ=255)
                 tGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAAt                                                                                         >  1:1325645/1‑100 (MQ=255)
                 tGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAAt                                                                                         >  1:3256633/1‑100 (MQ=255)
                  gCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATc                                                                                        >  1:1916829/1‑100 (MQ=255)
                  gCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATc                                                                                        >  1:3011338/1‑100 (MQ=255)
                   gaCGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACACAGTACGCTATGTTTATGCCGAATGCGTCTTTAACGCCTTATCCGTCCT‑‑Acgca                                                                                            >  2:1485222/2‑93 (MQ=25)
                   cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCg                                                                                       <  1:3510278/100‑1 (MQ=255)
                   cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCg                                                                                       <  2:3814777/100‑1 (MQ=255)
                    aCGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGt                                                                                      <  2:834424/100‑1 (MQ=255)
                     cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGt                                                                                      >  2:1686935/1‑99 (MQ=255)
                     cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTg                                                                                     <  2:814803/100‑1 (MQ=255)
                     cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTg                                                                                     <  1:3787153/100‑1 (MQ=255)
                          agGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAAt                                                                                <  1:3620351/100‑1 (MQ=255)
                            gCACAGCGCTACCGGGTATTCAGATGGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTc                                                                              <  2:1924443/100‑1 (MQ=255)
                            gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTc                                                                              <  1:1653585/100‑1 (MQ=255)
                             cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCa                                                                             >  1:1399735/1‑100 (MQ=255)
                              acaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCaa                                                                            <  1:1145316/100‑1 (MQ=255)
                               caGCGCCACCGGGTATACAGATTGTTGACTAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGTCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAAt                                                                           <  2:2364817/100‑1 (MQ=255)
                                aGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATa                                                                          <  2:128994/100‑1 (MQ=255)
                                     cACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTg                                                                     <  2:1773465/100‑1 (MQ=255)
                                     cACCGGGTATACAAATTGTTGACAAAGTGCGATTTGTTTATGGCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTg                                                                     >  1:2662340/1‑100 (MQ=255)
                                         gggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑Ga                                                               >  1:819259/1‑100 (MQ=255)
                                             ataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAt‑tt                                                          <  2:433215/100‑1 (MQ=255)
                                              taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAt‑ttt                                                         >  1:1104764/1‑100 (MQ=255)
                                              taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAt‑ttt                                                         >  1:163557/1‑100 (MQ=255)
                                                cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑g                                                     <  2:2290613/100‑1 (MQ=255)
                                                cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑g                                                     <  2:3715525/100‑1 (MQ=255)
                                                    ttgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAgg                                                 <  1:698438/100‑1 (MQ=255)
                                                     cataaaCAAAGCGCACTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGc                                                >  2:3227649/6‑100 (MQ=18)
                                                      gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGcc                                               <  2:2211694/100‑1 (MQ=255)
                                                       taaaCAAAGCGCACTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCt                                              >  2:1983579/4‑100 (MQ=18)
                                                         gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGa                                            >  2:2556248/1‑100 (MQ=255)
                                                             aaGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATAAg                                        <  1:1216033/100‑1 (MQ=255)
                                                              aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGc                                       <  1:723624/100‑1 (MQ=255)
                                                                 gcgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGa                                    >  2:1394693/1‑100 (MQ=255)
                                                                      ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCa                               <  2:1123426/100‑1 (MQ=255)
                                                                        gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATc                             <  2:3528445/100‑1 (MQ=255)
                                                                        gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATc                             >  2:2122604/1‑100 (MQ=255)
                                                                          gtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAACATATTG‑CA‑ATTC‑TCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                           >  1:2104335/5‑100 (MQ=18)
                                                                          gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAAATCTTGCCAATTCAATATATTG‑CA‑GGAC‑CAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                           <  1:2324286/96‑1 (MQ=18)
                                                                           taggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAAATCTTGCCAATTCAATATATTG‑CA‑GGAC‑CAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                          >  1:890024/4‑100 (MQ=18)
                                                                           taggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAAGTCGTATAAATTCAATGTATTG‑CA‑GTTT‑ATTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAgg                          >  1:450411/4‑100 (MQ=25)
                                                                           tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAgg                          <  2:3018724/100‑1 (MQ=255)
                                                                            aggCCGGATAAGGAGATCACGCCGCATCCGGCCT‑‑AC‑AATACTGTGCAAATTCAATATATTG‑CA‑GGAA‑ACAC‑‑GTAGGCCTGATAAGCCAAGCGCATCAGGc                         >  2:2434693/3‑100 (MQ=17)
                                                                             ggCCGGATAAGGGGCTTACGCCGCATCCGGCAT‑‑AC‑AAGAACGTGCAAATTCAACAAATTGTGATGCAT‑CCT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                         >  2:746507/2‑100 (MQ=18)
                                                                             ggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑AC‑AAATTC‑TGCAAATTCAACAAATTG‑C‑‑GAATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAt                       >  1:2442379/2‑100 (MQ=18)
                                                                             ggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑AC‑AAAAGCATGCAAATTCAATATATTG‑CA‑GAGT‑TCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  2:1631156/2‑100 (MQ=17)
                                                                             ggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑AC‑AAAAGCATGCAAATTCAATATATTG‑CA‑GAGT‑TCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  1:840017/2‑100 (MQ=17)
                                                                             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAAATCTTGCCAATTCAATATATTG‑CA‑GGAC‑CAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  1:599739/2‑100 (MQ=18)
                                                                             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAGATCGTGCAAATTCAATATATTG‑CA‑TGAG‑CCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  2:135189/2‑100 (MQ=17)
                                                                             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAAACTTTGCAAATTCAATATATTG‑C‑‑‑AAT‑CTCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  1:2194927/2‑100 (MQ=18)
                                                                             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAAACTTTGCAAATTCAATATATTG‑C‑‑‑AAT‑CTCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  1:1823785/2‑100 (MQ=18)
                                                                             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAAACCTTGCTAATTCAATATATTG‑CAGGGAC‑TAT‑‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCa                        >  1:2888249/2‑100 (MQ=18)
                                                                              gCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑AC‑AAAATCGTGCTAATTCAATATATTG‑CA‑GAAA‑CCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                       >  1:2328186/1‑99 (MQ=18)
                                                                              gCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAAATCTTGCCAATTCAATATATTG‑CA‑GGAC‑CAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                       >  2:3093156/1‑99 (MQ=18)
                                                                              gCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAAAGCATGAAAATTCAATATATTG‑CA‑GGAG‑CTGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                       >  1:3892116/1‑99 (MQ=17)
                                                                              gCCGGATAAGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt                       >  2:202301/1‑100 (MQ=255)
                                                                               ctggataaggcgTTCACGCCGCATCCAGCCT‑‑AC‑AAAATTGTGCAAAATCAATGGATTG‑CA‑CAGC‑TAGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaga                      >  1:2728497/9‑98 (MQ=25)
                                                                               ccggataaggcgTTTACGCCGCATCCGGCCT‑‑AC‑AAAAGCATGCAAATTCAATATATTG‑CA‑GAGT‑TCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagt                      >  1:1746472/9‑98 (MQ=17)
                                                                               ccggataaggcgTTCACGCCGCATCCGGCCT‑‑AC‑AAAATCGCTTAAATTCAATATATTG‑CA‑CCAT‑CCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaat                      >  1:915809/9‑98 (MQ=21)
                                                                               ccggataaggcgTTCACGCCGCATCCGGCCT‑‑AC‑AAAAGCATGAAAATTCAATATATTG‑CA‑GGAG‑CTGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagt                      >  1:3074817/9‑98 (MQ=18)
                                                                               ccggataaggcgTTCACGCCGCATCCGGCCT‑‑AC‑AAAAGCATGAAAATTCAATATATTG‑CA‑GGAG‑CTGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagt                      >  1:302215/9‑98 (MQ=17)
                                                                               ccGGATAAGGCGTTCACGCCGCATCCGGCATGAAC‑AAAAACATGCAAATTCAATAGATTG‑CA‑GAGA‑TTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  1:201161/1‑100 (MQ=17)
                                                                                 ggataaggcgTTTACGCCGCATCTGGCCT‑‑AC‑AAATTCTTCCAAATTCGATATATTG‑CT‑GAAA‑TCAT‑‑GTAGGCCTGATAAGCGTATCGCATCAGGCAattt                    >  1:148022/7‑100 (MQ=17)
                                                                                 ggataaggcgTTCACGCCGCATCCGGCCT‑‑ACAAAAATCTTGCCAATTCAATATATTG‑CA‑GGAC‑CAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                    >  2:1516479/7‑100 (MQ=18)
                                                                                 ggataaggcgTTCACGCCGCATCCGGCCT‑‑AC‑AAAGTCGTATAAATTCAATGTATTG‑CA‑GTTT‑ATTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaatct                    >  1:3203297/7‑96 (MQ=21)
                                                                                 ggataaggcgTTCACGCCGCATCCGGCCT‑‑AC‑AAAAGCATGAAAATTCAATATATTG‑CA‑GGAG‑CTGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                    >  1:212303/7‑100 (MQ=18)
                                                                                 ggATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtttt                    <  2:3883383/100‑1 (MQ=255)
                                                                                  gataaggcgTTTACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCTAATTCAATATATTG‑CA‑GAAA‑CCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgtttt                   >  2:1845336/6‑100 (MQ=18)
                                                                                   aTGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                  <  1:1586966/100‑1 (MQ=255)
                                                                                    taaggcgTTTACGCCGCATCCGGCCT‑‑AC‑AAAAGCGTGCAAATTCAATACATTG‑CA‑TGGG‑CCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAattttac                 >  2:1429123/4‑98 (MQ=18)
                                                                                    tGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAACTCGTGTAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgc                 <  2:2668498/100‑1 (MQ=255)
                                                                                      cggcgTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct               >  1:36715/1‑100 (MQ=255)
                                                                                      cggcgTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct               >  1:2247401/1‑100 (MQ=255)
                                                                                         cgTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg            <  2:3595180/100‑1 (MQ=255)
                                                                                            tAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa         <  1:1116322/100‑1 (MQ=255)
                                                                                             aaCGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAt        >  2:1267549/1‑100 (MQ=255)
                                                                                              aCGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc       <  2:1180630/100‑1 (MQ=255)
                                                                                              aCGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc       >  2:1502175/1‑100 (MQ=255)
                                                                                               cGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCa      <  1:3749156/100‑1 (MQ=255)
                                                                                                gCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg     >  1:3155332/1‑100 (MQ=255)
                                                                                                gCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg     <  1:3345244/100‑1 (MQ=255)
                                                                                                gCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTG‑AA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg     <  2:257960/100‑1 (MQ=255)
                                                                                                  cTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTc   <  2:2867253/100‑1 (MQ=255)
                                                                                                   |                                                                                                           
TTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTG‑CA‑GAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCT  >  NC_000913/3705885‑3706083

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