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.7% 324.7 / 3.8 92coding (964/993 nt)rpoSRNA polymerase, sigma S (sigma 38) factor
Reads supporting (aligned to +/- strand):  ref base . (2/1);  new base C (56/33);  total (58/34)
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

TCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAG‑C‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACA  >  NC_000913/2866491‑2866682
                                                                                                  |                                                                                                
tCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAG‑c‑c                                                                                               >  2:1791328/1‑100 (MQ=255)
tCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAG‑c‑c                                                                                               >  1:1847416/1‑100 (MQ=255)
tCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAG‑c                                                                                                 >  2:2593729/1‑99 (MQ=255)
 cGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAG‑c‑c‑c                                                                                             <  1:719033/100‑1 (MQ=255)
     ccTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAG‑c‑c‑ccctg                                                                                         <  2:17888/100‑4 (MQ=255)
      cTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAG‑c‑c‑ccctgc                                                                                        >  2:1645111/1‑96 (MQ=255)
        tttATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAG‑C‑C‑CCTGCGtt                                                                                      <  1:726627/100‑1 (MQ=255)
            ttGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGc                                                                                   <  1:1078315/100‑1 (MQ=255)
            ttGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGc                                                                                   <  2:1280563/100‑1 (MQ=255)
             tgtgCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCa                                                                                  <  2:983071/100‑1 (MQ=255)
             tgtgCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCa                                                                                  >  1:694922/1‑100 (MQ=255)
             tgtgCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCa                                                                                  <  1:2487211/100‑1 (MQ=255)
             tgtgCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCa                                                                                  <  1:1648783/100‑1 (MQ=255)
             tgtgCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCa                                                                                  <  1:1237930/100‑1 (MQ=255)
              gtgCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAg                                                                                 <  2:1764984/100‑1 (MQ=255)
              gtgCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAg                                                                                 >  1:2147151/1‑100 (MQ=255)
               tgCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAgg                                                                                <  2:1046115/100‑1 (MQ=255)
                gCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAG‑CAC‑CCTGCGTTTGCAGGa                                                                               >  1:670721/1‑100 (MQ=255)
                 cacaGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGAt                                                                              >  2:1358423/1‑100 (MQ=255)
                 cacaGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGAt                                                                              <  1:828506/100‑1 (MQ=255)
                  acaGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTCCTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGAtt                                                                             >  1:1682501/1‑100 (MQ=255)
                   caGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGAttt                                                                            >  1:1386728/1‑100 (MQ=255)
                   caGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGAttt                                                                            >  1:1386715/1‑100 (MQ=255)
                   caGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGAttt                                                                            >  2:1218230/1‑100 (MQ=255)
                      aaaaGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTcgc                                                                         <  1:412069/100‑1 (MQ=255)
                      aaaaGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTcgc                                                                         <  1:1953006/100‑1 (MQ=255)
                       aaaGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTcgcg                                                                        <  1:134417/100‑1 (MQ=255)
                        aaGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCgcgc                                                                       <  2:812348/100‑1 (MQ=255)
                          ggCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCaa                                                                     >  1:2245676/1‑100 (MQ=255)
                          ggCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCaa                                                                     >  1:995770/1‑100 (MQ=255)
                             tccgatctcTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACg                                                                  <  1:675391/92‑1 (MQ=255)
                             cAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACg                                                                  >  1:2453361/1‑100 (MQ=255)
                               gCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGgc                                                                >  2:1653947/1‑100 (MQ=255)
                                 cTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGgcgc                                                              <  2:971109/100‑1 (MQ=255)
                                 cTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGgcgc                                                              >  2:802896/1‑100 (MQ=255)
                                   cGCTTGAGACTGGCCTTTCTGACGGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAg                                                            >  2:786248/1‑100 (MQ=255)
                                   cGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTATTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAg                                                            >  1:527324/1‑100 (MQ=255)
                                   cGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAg                                                            <  1:2551658/100‑1 (MQ=255)
                                    gCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAgg                                                           >  1:684233/1‑100 (MQ=255)
                                    gCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAgg                                                           >  1:572622/1‑100 (MQ=255)
                                     cTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACgagatcgga                                                          >  2:675387/1‑92 (MQ=255)
                                      ttGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGcc                                                         <  1:802149/100‑1 (MQ=255)
                                       tGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCt                                                        >  1:490208/1‑100 (MQ=255)
                                         agaCTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTc                                                      <  2:963560/100‑1 (MQ=255)
                                         agaCTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTc                                                      >  1:228212/1‑100 (MQ=255)
                                          gaCTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCa                                                     >  1:1278269/1‑100 (MQ=255)
                                          gaCTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCa                                                     <  1:193427/100‑1 (MQ=255)
                                          gaCTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCa                                                     <  1:73533/100‑1 (MQ=255)
                                          gaCTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCa                                                     >  2:919934/1‑100 (MQ=255)
                                           aCTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCaa                                                    >  1:2516250/1‑100 (MQ=255)
                                            ntGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCa                                                     >  2:2294313/2‑98 (MQ=255)
                                            cTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAAc                                                   >  2:427572/1‑100 (MQ=255)
                                             tGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAAcc                                                  >  2:228184/1‑100 (MQ=255)
                                             tGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAAcc                                                  <  1:1248784/100‑1 (MQ=255)
                                               gCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTg                                                >  1:264582/1‑100 (MQ=255)
                                                ccTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGa                                               >  2:2517896/1‑100 (MQ=255)
                                                 cTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGaa                                              >  1:488197/1‑100 (MQ=255)
                                                     cTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAG‑C‑CACCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTg                                          >  1:555132/1‑100 (MQ=255)
                                                       gACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGc                                        >  1:2532422/1‑100 (MQ=255)
                                                       gACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGc                                        <  1:2253459/100‑1 (MQ=255)
                                                         cAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGa                                      >  1:245112/1‑100 (MQ=255)
                                                           gATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAac                                    >  1:832620/1‑100 (MQ=255)
                                                           gATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAac                                    <  2:394081/100‑1 (MQ=255)
                                                             tGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGcctcacc                                  >  2:2410490/1‑94 (MQ=255)
                                                               cttacttaCTCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGt                                >  1:15773/1‑100 (MQ=255)
                                                                  acttactCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTc                              >  1:553995/1‑99 (MQ=255)
                                                                   cttactCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCAc                            >  1:2323835/1‑100 (MQ=255)
                                                                    ttactCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACg                           <  1:238720/100‑1 (MQ=255)
                                                                      actCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACggg                         <  1:1284996/100‑1 (MQ=255)
                                                                       ctCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGt                        >  2:1506308/1‑100 (MQ=255)
                                                                        tCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTg                       <  2:53020/100‑1 (MQ=255)
                                                                        tCGCGGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTg                       <  2:1228330/100‑1 (MQ=255)
                                                                         cgcgGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGa                      >  2:1487927/1‑100 (MQ=255)
                                                                         cgcgGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGa                      >  1:2583066/1‑100 (MQ=255)
                                                                         cgcgGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGACGAACACGTTCACGGGTGa                      >  1:1480184/1‑100 (MQ=255)
                                                                          gcgGAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTggg                     >  1:2145973/1‑98 (MQ=255)
                                                                            ggAACAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGc                   <  2:363302/100‑1 (MQ=255)
                                                                               aCAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCaa                >  1:2324797/1‑100 (MQ=255)
                                                                                cAGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGACTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGCGTGCGGCCAAt               >  2:264032/1‑100 (MQ=255)
                                                                                 ngcgcTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAAtt              >  1:894656/2‑100 (MQ=255)
                                                                                 aGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAAtt              >  2:523565/1‑100 (MQ=255)
                                                                                 aGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAAt               >  2:1478154/1‑99 (MQ=255)
                                                                                 aGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAAACTGAATCTGGAGAATACGTTCACGGGTGAGGCCAttt              >  2:1975419/1‑100 (MQ=255)
                                                                                 aGCGCTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCACAAACGGCGCAGGCCTTCACCCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAAtt              >  2:1422974/1‑100 (MQ=255)
                                                                                  gcgcTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAAttt             >  1:2420792/1‑100 (MQ=255)
                                                                                    gcTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCa           <  1:1764162/100‑1 (MQ=255)
                                                                                     cTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCAc          <  2:2245663/100‑1 (MQ=255)
                                                                                     cTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCAc          >  1:1547518/1‑100 (MQ=255)
                                                                                     cTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCAc          <  1:523569/100‑1 (MQ=255)
                                                                                     cTTCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCAc          >  1:695608/1‑100 (MQ=255)
                                                                                       tCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGa        >  2:1228633/1‑100 (MQ=255)
                                                                                       tCGATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGa        >  1:890538/1‑100 (MQ=255)
                                                                                         gATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGAcc      >  2:1426987/1‑100 (MQ=255)
                                                                                         gATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGAcc      >  2:137299/1‑100 (MQ=255)
                                                                                         gATATTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGAcc      >  1:1389848/1‑100 (MQ=255)
                                                                                          atatTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCt     <  1:2391333/100‑1 (MQ=255)
                                                                                           tatTCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTa    >  2:1639156/1‑100 (MQ=255)
                                                                                             ttCAGCC‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACa  <  2:1480176/100‑1 (MQ=255)
                                                                                                  |                                                                                                
TCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAG‑C‑C‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACA  >  NC_000913/2866491‑2866682

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: ‑