Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 520,993 T→G 100% 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.7% 513.8 / ‑2.8 148V193G (GTC→GGC) ybbPputative ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base G (71/75);  total (73/75)
Fisher's exact test for biased strand distribution p-value = 2.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.65e-01

Only 100 of 149 total aligned reads displayed.

CGGATTCCGGTTTTAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGTCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGAACAACGCTGGTACGG  >  NC_000913/520897‑521089
                                                                                                |                                                                                                
cGGATTCCGGTTTTAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCAc                                                                                               >  1:2098413/1‑100 (MQ=255)
     tCCGGTTTTAACCCCTTCCAGATGTCTCCGCTTCTGATGATTAATCTTGTGGATTTCGATACAACCGGAGCCGTGCAGCCGGGGCGGCGGGGCACCTgg                                                                                           >  2:2614848/1‑99 (MQ=255)
     tCCGGTTTTAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGCGTCGGGGCACCTGGc                                                                                          >  2:1074449/1‑100 (MQ=255)
     tCCGGTTTTAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGc                                                                                          <  1:634290/100‑1 (MQ=255)
     tCCGGGTTTAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGc                                                                                          >  2:2366662/1‑100 (MQ=255)
      ccGGTTTTAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCg                                                                                         >  2:970452/1‑100 (MQ=255)
      ccGGTTTTAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCg                                                                                         >  1:3910790/1‑100 (MQ=255)
      ccGGTTTTAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCg                                                                                         >  1:3767357/1‑100 (MQ=255)
      ccGGTTTTAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCg                                                                                         >  1:3742580/1‑100 (MQ=255)
      ccGGTTTTAACCCCTTCCAGATGGCCCCGCGTCTGACGATGAAACTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGGCGCGGGGCACCTGGCg                                                                                         >  2:1611927/1‑100 (MQ=255)
        ggTTTTAACCCCTTCCAGATGGCTCCGCGTCTGATGGTGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGtt                                                                                       <  1:1570040/100‑1 (MQ=255)
        ggTTTTAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGtt                                                                                       <  1:1144786/100‑1 (MQ=255)
         gTTTTAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTta                                                                                      <  1:3430562/100‑1 (MQ=255)
         gTTTTAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTta                                                                                      <  2:654966/100‑1 (MQ=255)
           tttAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTtata                                                                                    <  2:2817607/100‑1 (MQ=255)
           tttAACCCCTTCCAGATGGCTCCGCGTCTGATGAGGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTgata                                                                                    <  2:2785705/100‑1 (MQ=255)
            ttAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATaa                                                                                   <  1:2756885/100‑1 (MQ=255)
             tAACCCCTTGCAGATGGCTCCGCGTCTGATGTTGAACCCGGCCACTGTCGTTCAAGCCGGAACCGGGCAGCCGGGGCGAGGGGGCACCtcacg                                                                                         >  2:3483097/1‑89 (MQ=255)
             tAACCCCTTCCAGATGGCTGCGCGTCTGATGATGAATCTGGAGGATGTCGATAAAAGCGGAGCAGTGAAGCCGGGGAGACGGGGCACCCGGAGTta                                                                                      >  2:3955591/1‑96 (MQ=255)
             tAACCCCTTCCAGATGGCTCCGCGTCTGATGGTGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATaaa                                                                                  >  1:1773859/1‑100 (MQ=255)
             tAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATaaa                                                                                  <  1:4166594/100‑1 (MQ=255)
               aCCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCCGACCGCCGTCGATAAACCCGGAGCCGCGCACCCCGGGAGACGGGGCACCTGGcc                                                                                         >  2:2559669/1‑90 (MQ=255)
                ccccTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTc                                                                               >  2:605334/1‑100 (MQ=255)
                ccccTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGAAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTc                                                                               >  2:928296/1‑100 (MQ=255)
                 cccTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTcg                                                                              <  1:734297/100‑1 (MQ=255)
                   cTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTcggc                                                                            >  1:2721074/1‑100 (MQ=255)
                   cTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTcggc                                                                            >  1:3121575/1‑100 (MQ=255)
                   cTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTcggc                                                                            <  2:2188956/100‑1 (MQ=255)
                    ttCCAGATGTCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGACGGGGAGTCGGGGCACCTGGCGTTATAAATTcggcg                                                                           <  2:776119/100‑1 (MQ=255)
                     tCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTcggcgg                                                                          <  2:839862/100‑1 (MQ=255)
                     tCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGAATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTcggcgg                                                                          <  2:646038/100‑1 (MQ=255)
                       cAGATGCCTCCGACTCTGATGAACACTCTGGAGAATGTCGATAAAACCGGAGCCGTGCAGCAGCGGAACTGGGGCACCTGTCGTGATAAGTTcggc                                                                            >  2:2423983/1‑96 (MQ=255)
                         gATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAAc                                                                      >  2:1540078/1‑100 (MQ=255)
                          aTGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACg                                                                     <  2:4118460/100‑1 (MQ=255)
                             gCTCCGCGTCTGATGCTGAATCTGTCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCTCCTGGCGTTATAAATTCGGCGGCAACgaga                                                                  <  2:3163093/100‑1 (MQ=255)
                             gCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACgaga                                                                  <  1:2251522/100‑1 (MQ=255)
                             gCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGCGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACgagn                                                                  <  2:2269096/100‑2 (MQ=255)
                              cTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGaa                                                                 >  2:3343553/1‑100 (MQ=255)
                              cTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGaa                                                                 <  1:970456/100‑1 (MQ=255)
                               tCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGGGGCAACGAGAAc                                                                <  2:171980/100‑1 (MQ=255)
                                ccGCTTCTGATGATGAATCGGGCGGATGTCTATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAAcc                                                               <  1:2853670/100‑1 (MQ=255)
                                 cgcgTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAAcc                                                               <  2:1475747/99‑1 (MQ=255)
                                 cgcgTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTACAGCCGGGGAGTCGGGGCACCCGGCGCTATACATTCGGCGGCGACGAGAACCa                                                              >  2:1702862/1‑100 (MQ=255)
                                   cgTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGc                                                            <  1:1525353/100‑1 (MQ=255)
                                     tCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGTCACCTGGCGTTATAAATTCGGCGGCAACGGGAACCAGCTc                                                          >  2:335378/1‑100 (MQ=255)
                                        gatgatgaATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGAc                                                       >  1:3249384/1‑100 (MQ=255)
                                         atgatgaATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACg                                                      <  2:1885187/100‑1 (MQ=255)
                                         atgatgaATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACg                                                      >  1:4116561/1‑100 (MQ=255)
                                          tgatgaATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACgg                                                     >  1:2565038/1‑100 (MQ=255)
                                           gatgaATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGc                                                    >  2:1018266/1‑100 (MQ=255)
                                           gatgaATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTCATAAATTCGGCGGCAACGAGAACCAGCTCGACGGc                                                    <  1:934159/100‑1 (MQ=255)
                                             tgaATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTa                                                  <  1:236378/100‑1 (MQ=255)
                                             tgaATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTa                                                  <  1:2989633/100‑1 (MQ=255)
                                               aaTCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATg                                                <  2:909366/100‑1 (MQ=255)
                                               aaTCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATg                                                >  1:1924202/1‑100 (MQ=255)
                                               aaTCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATg                                                >  1:2501358/1‑100 (MQ=255)
                                                 tCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATga                                               <  2:3375277/99‑1 (MQ=255)
                                                    ggCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCGCGACGGCTATGAGgaa                                           >  2:2754266/1‑100 (MQ=255)
                                                     gCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGaaan                                          <  2:3053688/100‑2 (MQ=255)
                                                     gCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAAt                                          <  2:550353/100‑1 (MQ=255)
                                                       ggATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATgg                                        <  2:2640053/100‑1 (MQ=255)
                                                       ggAGGTCGATAAAACCGGAGCCGGGTATCCGGGGAGTCGGGGCACCTCGCGGTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATg                                         <  2:3666312/99‑1 (MQ=255)
                                                        gATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGgt                                       <  2:254667/100‑1 (MQ=255)
                                                         aTGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGgtt                                      >  2:1768003/1‑100 (MQ=255)
                                                         aTGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGgtt                                      <  1:117386/100‑1 (MQ=255)
                                                          tGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGgttg                                     >  1:2479510/1‑100 (MQ=255)
                                                          tGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGgttg                                     <  2:1432050/100‑1 (MQ=255)
                                                           gTCGATAAAACCGGAGCCGTGCAGCCGGGGCGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATgg                                        >  2:1605551/1‑96 (MQ=255)
                                                           gTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTGGACGGCTATGAGAAATGgttg                                     >  1:1706234/1‑99 (MQ=255)
                                                           gTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGgttgt                                    <  1:605338/100‑1 (MQ=255)
                                                             cgctAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTAGCGTTATAAATTCGGCGGCAACGAGAAGCAGCTCGACGGCTATGAGACATGgtactta                                  >  2:2682416/4‑95 (MQ=255)
                                                              gATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTAc                                 >  1:2201358/1‑100 (MQ=255)
                                                              gATAAAACCGGAGCCGTACAGCCGGGGAGTCGGGGCACCGGGCGATATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGGGAAGTGggtgt                                    >  2:269706/1‑97 (MQ=255)
                                                               aTAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTAcc                                <  1:515767/100‑1 (MQ=255)
                                                                 aaaaCCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTc                              <  2:3656852/100‑1 (MQ=255)
                                                                  aaaCCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCa                             <  2:3131283/100‑1 (MQ=255)
                                                                  aaaCCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCa                             <  2:3422091/100‑1 (MQ=255)
                                                                  aaaCCGGAGCCGTGCAGCCGGGGAGTCGGGGCAACTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCa                             >  2:3812580/1‑100 (MQ=255)
                                                                   aaCCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCCGCTCGACGGCTATGAGAAATGGTTGTTACCTCAg                            >  1:3833088/1‑100 (MQ=255)
                                                                   aaCCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAg                            <  1:528691/100‑1 (MQ=255)
                                                                    aCCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTCGGTACCTCAGc                           >  2:1290989/1‑100 (MQ=255)
                                                                     ccGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCt                          >  1:2135536/1‑100 (MQ=255)
                                                                      cGGAGCCTTGCAGCCGGGGAGTCGAGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCtt                         <  2:3800521/100‑1 (MQ=255)
                                                                      cGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGGCGGCTATGAGAAAGGGTGGTTACCTCAGCtt                         >  1:3323313/1‑100 (MQ=255)
                                                                       ggAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTa                        <  2:125186/100‑1 (MQ=255)
                                                                          gCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAAc                     <  1:1158992/100‑1 (MQ=255)
                                                                           ccGTGCAGCCGGGGAGTCGGGGCACCTGGCATTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAAcc                    <  2:3487944/100‑1 (MQ=255)
                                                                            cGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAAccc                   >  1:2908538/1‑100 (MQ=255)
                                                                            cGTGCAGCCGGGGAGTCGGGGCACCTGGCGTCATAAATTCGGGGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAAccc                   >  1:2819379/1‑100 (MQ=255)
                                                                               gCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGaa                >  1:2699228/1‑100 (MQ=255)
                                                                               gCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCCCAGCTTAAACCCGaa                >  2:3271335/1‑100 (MQ=255)
                                                                                  gCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGaacaa             >  1:1925954/1‑100 (MQ=255)
                                                                                    cGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGAACAACg           <  1:3118384/100‑1 (MQ=255)
                                                                                     ggggAGTCGGTGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAACTGCTTGTTACCTCAGCTTAAACCCGAACAACGc          <  2:2746858/100‑1 (MQ=255)
                                                                                     ggggAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGAACAACGc          <  2:3216206/100‑1 (MQ=255)
                                                                                       ggAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGAACAACGCTg        >  1:2390131/1‑100 (MQ=255)
                                                                                       ggAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGAACAACGCTg        >  1:804457/1‑100 (MQ=255)
                                                                                        gAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGAACAACGCTgg       <  2:2889730/100‑1 (MQ=255)
                                                                                          gTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGAACAACGCTGGTa     <  1:3591961/100‑1 (MQ=255)
                                                                                          gTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCCATGAGAGATGGGTGTTACCTCAGCTTAAACCCGAACAACGCAGGTa     >  2:731688/1‑100 (MQ=255)
                                                                                                |                                                                                                
CGGATTCCGGTTTTAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGTCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGAACAACGCTGGTACGG  >  NC_000913/520897‑521089

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