Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,866,592 +C 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.C98.7% 570.6 / ‑0.6 158coding (964/993 nt)rpoSRNA polymerase, sigma S (sigma 38) factor
Reads supporting (aligned to +/- strand):  ref base . (2/0);  new base C (79/77);  total (81/77)
Fisher's exact test for biased strand distribution p-value = 4.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 168 total aligned reads displayed.

TCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAG‑CC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACA  >  NC_000913/2866491‑2866682
                                                                                                  |                                                                                               
 cGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAG‑cc‑c                                                                                             <  1:1877897/100‑1 (MQ=255)
   gACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAG‑cc‑ccc                                                                                           >  2:253738/1‑99 (MQ=255)
   gACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAG‑cc‑ccc                                                                                           >  1:1739014/1‑99 (MQ=255)
     ccTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAG‑cc‑ccctg                                                                                         >  1:1679780/1‑97 (MQ=255)
      cTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAG‑cc‑ccctgc                                                                                        >  1:2208730/1‑96 (MQ=255)
      cTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAG‑cc‑ccctgc                                                                                        >  1:1580293/1‑96 (MQ=255)
       ttttATTGTGCACAGAAAAGGCCAGACTCGCGTGAGACTGGCCTTGCTGACAGATGCTTACGTAGTCACGGAAGAACGACTTGATATTGAG‑c                                                                                                >  2:2963248/1‑92 (MQ=255)
           aTTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTg                                                                                    <  1:196242/100‑1 (MQ=255)
              gtgCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAg                                                                                 <  2:1297219/100‑1 (MQ=255)
               tgCACAGAAAAGGCCAGCCTCGCTTGCGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAg                                                                                 <  2:150047/99‑1 (MQ=255)
               tgCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAgg                                                                                >  1:2472301/1‑100 (MQ=255)
               tgCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAgg                                                                                >  2:2468881/1‑100 (MQ=255)
                gCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGa                                                                               >  2:2889067/1‑100 (MQ=255)
                 cacaGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGAt                                                                              >  1:15433/1‑100 (MQ=255)
                  acaGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGAtt                                                                             <  1:43589/100‑1 (MQ=255)
                   caGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGAttt                                                                            <  2:221466/100‑1 (MQ=255)
                     gAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTcg                                                                          >  2:1855173/1‑100 (MQ=255)
                      aaaaGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTcgc                                                                         <  2:1117236/100‑1 (MQ=255)
                           gCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCaaa                                                                    >  1:1862532/1‑100 (MQ=255)
                             cAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACg                                                                  <  2:3261975/100‑1 (MQ=255)
                                ccTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTCCTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGgc                                                                >  2:723207/1‑99 (MQ=255)
                                ccTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGgcg                                                               >  2:723015/1‑100 (MQ=255)
                                   cGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAg                                                            <  2:1749459/100‑1 (MQ=255)
                                   cGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAg                                                            <  1:1360471/100‑1 (MQ=255)
                                    gCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGCTATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAgg                                                           >  1:84191/1‑100 (MQ=255)
                                    gCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAgg                                                           <  1:285132/100‑1 (MQ=255)
                                     cTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGc                                                          >  2:278480/1‑100 (MQ=255)
                                     cTTGAGAATGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGc                                                          >  1:1771314/1‑100 (MQ=255)
                                     cTTAAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGc                                                          >  2:2509219/1‑100 (MQ=255)
                                       tGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCt                                                        >  2:1313589/1‑100 (MQ=255)
                                       tGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCt                                                        <  1:1658172/100‑1 (MQ=255)
                                       tGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCt                                                        <  2:1031055/100‑1 (MQ=255)
                                        gaggcTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCtt                                                       <  2:1143590/96‑1 (MQ=255)
                                         agaCTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAATGGCGCAGGCCTTc                                                      <  1:576340/100‑1 (MQ=255)
                                         agaCTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTc                                                      <  1:644613/100‑1 (MQ=255)
                                          gaCTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCa                                                     >  1:2423158/1‑100 (MQ=255)
                                          gaCTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCa                                                     <  2:2373670/100‑1 (MQ=255)
                                          gaCTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCa                                                     >  1:2749004/1‑100 (MQ=255)
                                           aCTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCaa                                                    >  1:1109637/1‑100 (MQ=255)
                                           aCTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCaa                                                    <  1:2134646/100‑1 (MQ=255)
                                            ntGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAAc                                                   >  1:152850/2‑100 (MQ=255)
                                            cTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAAc                                                   >  2:3094266/1‑100 (MQ=255)
                                             tGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAAcc                                                  <  2:1580289/100‑1 (MQ=255)
                                             tGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAAcc                                                  <  2:1403849/100‑1 (MQ=255)
                                              ggCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCt                                                 >  1:1675493/1‑100 (MQ=255)
                                              ggCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCt                                                 <  2:1918602/100‑1 (MQ=255)
                                              ggCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCt                                                 >  1:3009175/1‑100 (MQ=255)
                                               gCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTg                                                <  1:407308/100‑1 (MQ=255)
                                                 cTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGaa                                              >  1:825702/1‑100 (MQ=255)
                                                 cTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGaa                                              <  1:42638/100‑1 (MQ=255)
                                                      tGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTgg                                         >  2:2913610/1‑100 (MQ=255)
                                                      tGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTgg                                         <  2:3087227/100‑1 (MQ=255)
                                                      tGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTgg                                         >  2:1394729/1‑100 (MQ=255)
                                                      tGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTgg                                         <  2:1605611/100‑1 (MQ=255)
                                                        aCAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCg                                       <  2:942567/100‑1 (MQ=255)
                                                         cAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGa                                      >  1:2680926/1‑100 (MQ=255)
                                                          aGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGa                                      <  2:2059170/99‑1 (MQ=255)
                                                           gATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGaa                                     >  2:21651/1‑99 (MQ=255)
                                                           gATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAac                                    >  2:489738/1‑100 (MQ=255)
                                                           gATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAac                                    <  1:525089/100‑1 (MQ=255)
                                                            aTGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAaca                                   >  1:1172161/1‑100 (MQ=255)
                                                             tGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAacac                                  >  1:2331058/1‑100 (MQ=255)
                                                             tGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAacac                                  >  1:2383948/1‑100 (MQ=255)
                                                               cttacttaCTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGt                                >  2:2858255/1‑100 (MQ=255)
                                                               cttacttaCTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGt                                >  1:786656/1‑100 (MQ=255)
                                                               cttacttaCTCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGt                                >  1:3257983/1‑100 (MQ=255)
                                                                    ttactCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACg                           <  2:1739009/100‑1 (MQ=255)
                                                                     tactCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACgg                          >  2:592343/1‑100 (MQ=255)
                                                                     tactCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACgg                          <  2:525394/100‑1 (MQ=255)
                                                                     tactCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACgg                          <  1:2440372/100‑1 (MQ=255)
                                                                       ctCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGt                        >  2:1515377/1‑100 (MQ=255)
                                                                       ctCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGt                        <  1:253739/100‑1 (MQ=255)
                                                                       ctCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGt                        >  1:2892252/1‑100 (MQ=255)
                                                                        tCGCGGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAAGCTGGCGAACACGTTCACGGGTg                       <  2:1943388/100‑1 (MQ=255)
                                                                         cgcgGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGa                      >  2:2595876/1‑100 (MQ=255)
                                                                         cgcgGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGa                      >  2:601845/1‑100 (MQ=255)
                                                                         cgcgGAACAGCGCTTCGATATTCAG‑CCACCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGa                      >  2:1214064/1‑100 (MQ=255)
                                                                           cgGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAgg                    <  1:1421361/100‑1 (MQ=255)
                                                                           cgGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAgg                    <  2:564215/100‑1 (MQ=255)
                                                                           cgGAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAgg                    <  2:1172158/100‑1 (MQ=255)
                                                                            ggAGCAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGc                   <  2:2431040/100‑1 (MQ=255)
                                                                            ggAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGc                   <  1:3098622/100‑1 (MQ=255)
                                                                            ggAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGc                   <  2:202578/100‑1 (MQ=255)
                                                                            ggAACAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGc                   >  2:508423/1‑100 (MQ=255)
                                                                                cAGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAAt               <  2:2317227/100‑1 (MQ=255)
                                                                                 aGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAAtt              >  1:134258/1‑100 (MQ=255)
                                                                                 aGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAAtt              <  2:2456144/100‑1 (MQ=255)
                                                                                 aGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAAtt              <  2:1297047/100‑1 (MQ=255)
                                                                                 aGCGCTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAAtt              >  2:298463/1‑100 (MQ=255)
                                                                                    gcTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCa           <  2:2208719/100‑1 (MQ=255)
                                                                                    gcTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCa           >  1:3027418/1‑100 (MQ=255)
                                                                                     cTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCAc          <  1:1507563/100‑1 (MQ=255)
                                                                                     cTTCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCAc          <  1:3198239/100‑1 (MQ=255)
                                                                                       tCGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGa        >  1:2945562/1‑100 (MQ=255)
                                                                                        cGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGAc       >  2:2284062/1‑100 (MQ=255)
                                                                                        cGATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGAc       >  2:2319170/1‑100 (MQ=255)
                                                                                         gATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGAcc      >  1:1941475/1‑100 (MQ=255)
                                                                                         gATATTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGAcc      <  2:2423147/100‑1 (MQ=255)
                                                                                          atatTCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCt     <  2:794140/100‑1 (MQ=255)
                                                                                             ttCAGCCC‑CCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACa  >  2:2079958/1‑100 (MQ=255)
                                                                                                  |                                                                                               
TCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAG‑CC‑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: ‑