Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 520,993 T→G V193G (GTC→GGC)  ybbP → putative ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913520,9930TG98.9% 335.5 / NA 92V193G (GTC→GGC) ybbPputative ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base G (47/44);  total (48/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGATTCCGGTTTTAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGTCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGAACAACGCTGGTACGG  >  NC_000913/520897‑521089
                                                                                                |                                                                                                
cGGATTCCGGTTTTAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGacac                                                                                               >  2:674574/1‑100 (MQ=255)
cGGATTCCGGTTTTAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCAc                                                                                               >  1:2190342/1‑100 (MQ=255)
    ttCCGGTTTTAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTgg                                                                                           >  2:1425149/1‑100 (MQ=255)
      ccGGTTTTAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCg                                                                                         >  2:1940370/1‑100 (MQ=255)
      ccGGTTTTAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCg                                                                                         >  1:2233709/1‑100 (MQ=255)
       cGGTTTTAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGt                                                                                        <  1:1797324/100‑1 (MQ=255)
       cGGTTTTAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGt                                                                                        <  1:2219671/100‑1 (MQ=255)
       cGGTTTTAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGt                                                                                        <  1:319793/100‑1 (MQ=255)
        ggTTTTAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGtt                                                                                       <  2:1131835/100‑1 (MQ=255)
          ttttAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTtat                                                                                     <  2:1415577/100‑1 (MQ=255)
          ttttAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTtat                                                                                     <  2:1006581/100‑1 (MQ=255)
          ttttAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTtat                                                                                     <  2:1614256/100‑1 (MQ=255)
           tttAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGTGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTtata                                                                                    <  2:1251678/100‑1 (MQ=255)
           tttAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTtata                                                                                    <  2:545532/100‑1 (MQ=255)
           tttAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTtata                                                                                    >  1:831618/1‑100 (MQ=255)
             tAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATaaa                                                                                  >  2:1097219/1‑100 (MQ=255)
             tAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATaaa                                                                                  <  2:1260166/100‑1 (MQ=255)
              aaCTCCTTCCAGATGGGGCCGCGTCTGATGATTAATCTGCTGAATGTCGTGAAAACGGGAGCCGTGCAACGGGGGAGGCGGGGCACATGGCGGtata                                                                                    >  2:1938943/1‑97 (MQ=255)
              aaCCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAAt                                                                                 >  1:90438/1‑100 (MQ=255)
                ccccTTCCAGATGGCTCCGCGTCTGATGATGAGTCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTc                                                                               >  1:2213748/1‑100 (MQ=255)
                ccccTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTc                                                                               >  1:796405/1‑100 (MQ=255)
                ccccTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTc                                                                               >  2:1502892/1‑100 (MQ=255)
                ccccTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTc                                                                               >  1:258903/1‑100 (MQ=255)
                ccccTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTc                                                                               >  1:2239950/1‑100 (MQ=255)
                    ttCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTcggcg                                                                           <  2:2325434/100‑1 (MQ=255)
                     tCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGCTGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTcggcgg                                                                          <  1:1643343/100‑1 (MQ=255)
                     tCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTcggcgg                                                                          >  1:487206/1‑100 (MQ=255)
                      ccAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCggcggc                                                                         <  1:1801027/100‑1 (MQ=255)
                       cAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCa                                                                        >  1:1072189/1‑100 (MQ=255)
                         gATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAAc                                                                      >  1:1856839/1‑100 (MQ=255)
                          aTGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACg                                                                     <  1:1097221/100‑1 (MQ=255)
                            ggCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACgag                                                                   <  1:1425153/100‑1 (MQ=255)
                            ggCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACgag                                                                   >  2:235516/1‑100 (MQ=255)
                             ggtCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACgagn                                                                  <  2:2274518/98‑2 (MQ=255)
                             gCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACgaga                                                                  <  1:1317038/100‑1 (MQ=255)
                               tCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAAc                                                                >  1:1467336/1‑100 (MQ=255)
                                   cgTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGc                                                            <  1:297425/100‑1 (MQ=255)
                                   cgTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGc                                                            >  1:31822/1‑100 (MQ=255)
                                      cTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCg                                                         <  1:147699/100‑1 (MQ=255)
                                      cTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCg                                                         >  2:1768977/1‑100 (MQ=255)
                                        gatgatgaATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGAc                                                       >  1:490323/1‑100 (MQ=255)
                                         ctgatgaATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACg                                                      <  2:1869623/99‑1 (MQ=255)
                                         atgatgaATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACg                                                      <  2:669853/100‑1 (MQ=255)
                                          tgatgaATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACgg                                                     <  1:338540/100‑1 (MQ=255)
                                           gatgaATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACgag                                                    >  1:1869629/1‑98 (MQ=255)
                                           gatgaATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGc                                                    <  1:1940377/100‑1 (MQ=255)
                                             tgaATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCt                                                   >  1:573541/1‑99 (MQ=255)
                                             tgaATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTa                                                  <  2:1072187/100‑1 (MQ=255)
                                             tgaATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTa                                                  >  1:1138251/1‑100 (MQ=255)
                                             tgaATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTa                                                  >  2:240239/1‑100 (MQ=255)
                                             tgaATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTa                                                  >  1:1260483/1‑100 (MQ=255)
                                             tgaATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTa                                                  >  1:1945881/1‑100 (MQ=255)
                                             tgaATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGACGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGGCCCAGCTCGAGGGc                                                    >  2:437546/1‑98 (MQ=255)
                                              gaATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCCGGGCACCTGGCGTGACACATTCGGCGGCAACGAGAACCAGCTCGGCGGCGAt                                                 >  2:1624120/1‑100 (MQ=255)
                                                  cTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATgaga                                             <  2:2233701/100‑1 (MQ=255)
                                                  cTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATgaga                                             <  1:1768982/100‑1 (MQ=255)
                                                   tGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGaa                                            <  1:1847917/100‑1 (MQ=255)
                                                      cGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCATCTCGACGGCTATGAGAAATg                                         <  2:573540/100‑1 (MQ=255)
                                                      cGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATg                                         <  2:90436/100‑1 (MQ=255)
                                                      cGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATg                                         >  1:142536/1‑100 (MQ=255)
                                                          tGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCTCCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGgttg                                     <  2:2277990/100‑1 (MQ=255)
                                                          tGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGgttg                                     <  2:744533/100‑1 (MQ=255)
                                                           gTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGgttgt                                    <  2:107584/100‑1 (MQ=255)
                                                            tCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGgttgtt                                   >  1:1909369/1‑100 (MQ=255)
                                                            tCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGgttgtt                                   >  1:114556/1‑100 (MQ=255)
                                                             cGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGgttgtca                                  >  2:1422844/1‑98 (MQ=255)
                                                             cGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTa                                  <  1:571817/100‑1 (MQ=255)
                                                             cGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTa                                  >  1:1401850/1‑100 (MQ=255)
                                                              gATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTAc                                 >  2:72206/1‑100 (MQ=255)
                                                              gATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTAc                                 <  1:1835406/100‑1 (MQ=255)
                                                                 aaaaCCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTc                              <  2:1197724/100‑1 (MQ=255)
                                                                 aaaaCCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTc                              <  2:1404971/100‑1 (MQ=255)
                                                                 aaaaCCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTc                              <  2:2213740/100‑1 (MQ=255)
                                                                  aaaCCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCa                             >  1:811324/1‑100 (MQ=255)
                                                                  aaaCCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCa                             >  2:1704016/1‑100 (MQ=255)
                                                                   aaCCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAg                            >  1:187925/1‑100 (MQ=255)
                                                                     ccGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCt                          <  1:1704021/100‑1 (MQ=255)
                                                                      cGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCACCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCtt                         <  1:1645529/100‑1 (MQ=255)
                                                                       ggAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTa                        >  1:514813/1‑100 (MQ=255)
                                                                          gCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAAc                     >  1:789062/1‑100 (MQ=255)
                                                                            cGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAAccc                   <  2:796404/100‑1 (MQ=255)
                                                                             gTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCg                  >  2:370691/1‑100 (MQ=255)
                                                                               gCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGaa                <  2:31821/100‑1 (MQ=255)
                                                                                 aGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGaaca              <  2:2209406/100‑1 (MQ=255)
                                                                                  gCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGaacaa             >  1:26160/1‑100 (MQ=255)
                                                                                  gCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGaacaa             >  1:369348/1‑100 (MQ=255)
                                                                                    cGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGAACAACg           >  1:1903214/1‑100 (MQ=255)
                                                                                     ngggAGACGGGTCACGAGGCGTTATACATTCGGCGGCCACGAGAACCCACTCGACGGCTATGAGAAATGGTTGTTACCTCAGCATAAACCCGAACAACg           >  2:888917/2‑99 (MQ=255)
                                                                                     ggggAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGAACAACGc          <  1:235518/100‑1 (MQ=255)
                                                                                      gggAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGAACAACGCt         <  2:142534/100‑1 (MQ=255)
                                                                                         aGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGAACAACGCTGGt      >  2:1078453/1‑100 (MQ=255)
                                                                                          gTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGAACAACGCTGGTa     >  2:802915/1‑100 (MQ=255)
                                                                                          gTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGAACAACGCTGGTa     >  2:1476839/1‑100 (MQ=255)
                                                                                             ggggCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGAACAACGCTGGTACgg  >  1:1559307/1‑100 (MQ=255)
                                                                                                |                                                                                                
CGGATTCCGGTTTTAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGTCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGAACAACGCTGGTACGG  >  NC_000913/520897‑521089

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