Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,866,592:1 +C 100% coding (960/993 nt) rpoS ← RNA polymerase, sigma S (sigma 38) factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,866,5881.C96.4% 582.1 / 12.7 167coding (964/993 nt)rpoSRNA polymerase, sigma S (sigma 38) factor
Reads supporting (aligned to +/- strand):  ref base . (3/3);  new base C (92/69);  total (95/72)
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

Only 100 of 185 total aligned reads displayed.

CATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAG‑CC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACA  >  NC_000913/2866489‑2866682
                                                                                                    |                                                                                               
cATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAg                                                                                                  >  2:1589699/1‑100 (MQ=255)
cATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAg                                                                                                  >  2:2182572/1‑100 (MQ=255)
    ggACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAG‑cc‑cc                                                                                            >  1:1954767/1‑100 (MQ=255)
        cTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAG‑cc‑ccctgc                                                                                        >  1:731368/1‑96 (MQ=255)
        cTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAG‑cc‑ccctgc                                                                                        >  2:3251619/1‑96 (MQ=255)
          tttATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGt                                                                                       <  2:3021450/100‑1 (MQ=255)
           ttATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGtt                                                                                      <  2:1216895/100‑1 (MQ=255)
            tATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGtt                                                                                      <  2:1633015/99‑1 (MQ=255)
                gtgCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAg                                                                                 <  1:1383773/100‑1 (MQ=255)
                gtgCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAg                                                                                 >  1:2183021/1‑100 (MQ=255)
                 tgCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAgg                                                                                >  1:1247263/1‑100 (MQ=255)
                 tgCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAgg                                                                                >  1:2837153/1‑100 (MQ=255)
                  gCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGa                                                                               >  1:2492467/1‑100 (MQ=255)
                  gCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGa                                                                               >  1:3307163/1‑100 (MQ=255)
                      aGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTc                                                                           <  2:3133034/100‑1 (MQ=255)
                       gAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTcg                                                                          >  1:2152048/1‑100 (MQ=255)
                       gAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTcg                                                                          >  1:424493/1‑100 (MQ=255)
                       gAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTcg                                                                          <  2:1044646/100‑1 (MQ=255)
                       gAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTcg                                                                          >  1:1338385/1‑100 (MQ=255)
                          aaGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCgcgc                                                                       <  2:805992/100‑1 (MQ=255)
                          aaGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCgcgc                                                                       <  2:2767386/100‑1 (MQ=255)
                          aaGGCCAGCCTCGCTTGAGACTGGCCTGTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCgcgc                                                                       <  2:2441649/100‑1 (MQ=255)
                            ggCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCCCGCa                                                                      >  2:3392870/1‑99 (MQ=255)
                            ggCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAG‑CCACCTGCGTTTGCAGGCTTTCGCGCa                                                                      >  2:2426296/1‑99 (MQ=255)
                              ccAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAAc                                                                   >  2:1294370/1‑100 (MQ=255)
                              ccAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAAc                                                                   >  1:1816869/1‑100 (MQ=255)
                                   cTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTCCTCGCGGAACAGCGCTTCGTTATTCAGCCC‑CCTGCTTTTGCAGGATTTCGCGCAAACGgcgc                                                              <  1:504508/100‑1 (MQ=255)
                                     cGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCTCGCAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAg                                                            <  2:1165636/100‑1 (MQ=255)
                                     cGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAg                                                            >  2:1195633/1‑100 (MQ=255)
                                      gCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAgg                                                           <  1:3263144/100‑1 (MQ=255)
                                      gCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAgg                                                           >  1:2353886/1‑100 (MQ=255)
                                       cTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGc                                                          >  1:2788255/1‑100 (MQ=255)
                                       cTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGc                                                          >  1:647184/1‑100 (MQ=255)
                                        ttGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGcc                                                         <  1:3815093/100‑1 (MQ=255)
                                         tGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGcc                                                         <  2:740242/99‑1 (MQ=255)
                                           agaCTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTc                                                      >  1:1601978/1‑100 (MQ=255)
                                           agaCTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTc                                                      <  2:1395770/100‑1 (MQ=255)
                                            gaCTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCa                                                     <  2:495587/100‑1 (MQ=255)
                                            gaCTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCa                                                     <  1:1277157/100‑1 (MQ=255)
                                            gaCTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCa                                                     >  1:2375247/1‑100 (MQ=255)
                                             gcTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCaa                                                    <  1:3475552/99‑1 (MQ=255)
                                             aCTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCaa                                                    <  2:1964635/100‑1 (MQ=255)
                                              cTGGCCTTTCTGACAGATGCTTACTTCCTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAAc                                                   <  2:128075/100‑1 (MQ=255)
                                              cTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAAc                                                   >  2:2902942/1‑100 (MQ=255)
                                              cTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTCGCAGGCTTTCGCGCAAACGGCGCAGGCCTTCAAc                                                   >  2:2846563/1‑100 (MQ=255)
                                               tGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAAcc                                                  <  1:3780312/100‑1 (MQ=255)
                                                ngCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAAcc                                                  >  2:2189192/2‑99 (MQ=255)
                                                ggCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCt                                                 >  1:49125/1‑100 (MQ=255)
                                                 gCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTg                                                <  1:830824/100‑1 (MQ=255)
                                                 gCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTg                                                >  2:3293959/1‑100 (MQ=255)
                                                  ccTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGa                                               >  1:206783/1‑100 (MQ=255)
                                                   cTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGaa                                              <  2:2603904/100‑1 (MQ=255)
                                                      tCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATc                                            <  1:1707038/99‑1 (MQ=255)
                                                      tCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCt                                           >  1:1667088/1‑100 (MQ=255)
                                                        tGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTgg                                         >  2:1354226/1‑100 (MQ=255)
                                                          aCAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCg                                       <  2:3327333/100‑1 (MQ=255)
                                                          aCAGATCCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGTGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCg                                       <  2:617657/100‑1 (MQ=255)
                                                           cAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGa                                      <  2:2976079/100‑1 (MQ=255)
                                                             gATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAac                                    >  1:3393899/1‑100 (MQ=255)
                                                             gATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAac                                    >  2:1223899/1‑100 (MQ=255)
                                                               tGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAacac                                  >  2:2136648/1‑100 (MQ=255)
                                                                   tacttactCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTc                              >  1:2756415/1‑100 (MQ=255)
                                                                    ccgatctCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCa                             <  1:2809307/95‑1 (MQ=255)
                                                                       tactCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACgg                          <  1:3680287/100‑1 (MQ=255)
                                                                        actCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTTGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACggg                         <  1:3141069/100‑1 (MQ=255)
                                                                        actCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACggg                         >  1:1180311/1‑100 (MQ=255)
                                                                        actCGCGGAACAGCGCTTCGAGATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACggg                         <  2:3338913/100‑1 (MQ=255)
                                                                         ctCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGt                        >  1:15334/1‑100 (MQ=255)
                                                                          tCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTg                       >  2:866257/1‑100 (MQ=255)
                                                                           cgcgGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCaagatcgg                      >  2:2809287/1‑93 (MQ=255)
                                                                           cgcgGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGa                      >  2:20469/1‑100 (MQ=255)
                                                                           cgcgGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGa                      <  1:3616929/100‑1 (MQ=255)
                                                                             cgGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAgg                    >  1:1548643/1‑100 (MQ=255)
                                                                               gAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGcc                  >  1:2546212/1‑100 (MQ=255)
                                                                               gAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGcc                  >  2:3589778/1‑100 (MQ=255)
                                                                                 aCAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCaa                >  1:3111594/1‑100 (MQ=255)
                                                                                   aGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAAtt              <  1:2040673/100‑1 (MQ=255)
                                                                                   aGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAAtt              >  2:300470/1‑100 (MQ=255)
                                                                                    gcgcTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAAttt             >  2:1840261/1‑100 (MQ=255)
                                                                                     cgcTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCagatcg            >  1:1040223/1‑95 (MQ=255)
                                                                                      gcTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCa           <  2:517399/100‑1 (MQ=255)
                                                                                      gcTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCa           >  1:950572/1‑100 (MQ=255)
                                                                                      gcTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCa           <  2:1848500/100‑1 (MQ=255)
                                                                                      gcTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCa           >  2:2527442/1‑100 (MQ=255)
                                                                                      gcTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCa           <  1:1501651/100‑1 (MQ=255)
                                                                                       cTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCAc          <  1:1096882/100‑1 (MQ=255)
                                                                                       cTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCAc          <  2:1494296/100‑1 (MQ=255)
                                                                                       cTTCGATATGCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCTGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCAc          <  1:1647973/100‑1 (MQ=255)
                                                                                          cGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGAc       >  2:1718339/1‑100 (MQ=255)
                                                                                          cGATATTCAG‑CC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGAcc      >  2:3813557/1‑100 (MQ=255)
                                                                                           gATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAGACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGAcc      <  1:2677359/100‑1 (MQ=255)
                                                                                           gATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGAcc      >  1:2385204/1‑100 (MQ=255)
                                                                                           gATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACATGAATCTGGCGAACACGTTCACGGGTGCGGCCAATTTCACGAcc      >  2:3533834/1‑100 (MQ=255)
                                                                                             tataCAG‑CC‑ACTGCGTTTGCAGGATTTTGCGCAAACGGCGCAGGCCTTCAACCAGAATCTGGCGAACACGTTCACGGGTGAGGCCAATCTCACGACCTa    <  2:792850/99‑1 (MQ=255)
                                                                                             tatTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTa    <  2:3694552/100‑1 (MQ=255)
                                                                                             tatTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTa    <  1:3162928/100‑1 (MQ=255)
                                                                                              atTCAGCCC‑CCTGCCTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTAc   >  2:3572321/1‑100 (MQ=255)
                                                                                               ttCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACa  >  1:544199/1‑100 (MQ=255)
                                                                                               ttCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACa  >  1:3221002/1‑100 (MQ=255)
                                                                                               ttCAGCCC‑CCTGCGTCTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACa  <  1:3764137/100‑1 (MQ=255)
                                                                                                    |                                                                                               
CATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAG‑CC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACA  >  NC_000913/2866489‑2866682

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