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.3% 423.0 / ‑0.0 118coding (964/993 nt)rpoSRNA polymerase, sigma S (sigma 38) factor
Reads supporting (aligned to +/- strand):  ref base . (2/0);  new base C (57/59);  total (59/59)
Fisher's exact test for biased strand distribution p-value = 4.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.10e-01

Only 100 of 124 total aligned reads displayed.

CGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACA  >  NC_000913/2866492‑2866682
                                                                                                 |                                                                                               
cGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGccc                                                                                             <  2:1372248/100‑1 (MQ=255)
     cTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATCTTCAGccccctgc                                                                                        >  1:1345169/1‑96 (MQ=255)
      ttttATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGccccctgcg                                                                                       >  1:2665850/1‑95 (MQ=255)
      ttttATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGccccctgcg                                                                                       <  2:2946862/100‑6 (MQ=255)
      ttttATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGccccctgcg                                                                                       >  2:824719/1‑95 (MQ=255)
       tttATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGt                                                                                       <  1:2464619/100‑1 (MQ=255)
        ttATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGtt                                                                                      <  2:924622/100‑1 (MQ=255)
         tATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGttt                                                                                     >  2:2341848/1‑100 (MQ=255)
            tgtgCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCa                                                                                  <  2:508110/100‑1 (MQ=255)
             gtgCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGctg                                                                                 <  1:2199475/100‑3 (MQ=255)
              tgCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAgg                                                                                >  2:1547745/1‑100 (MQ=255)
              tgCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGAAAAATAAATAAATAATAGAGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAgg                                                                                <  2:3147530/100‑1 (MQ=255)
               gCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGa                                                                               >  2:527630/1‑100 (MQ=255)
                  caGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGAttt                                                                            >  2:1804511/1‑100 (MQ=255)
                  caGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGAttt                                                                            >  2:2739750/1‑100 (MQ=255)
                  caGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGAttt                                                                            >  1:383907/1‑100 (MQ=255)
                  caGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGAttt                                                                            <  2:2385408/100‑1 (MQ=255)
                   cagaaaGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTc                                                                           >  1:407310/4‑100 (MQ=255)
                    gcaaaGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAAGAGCGCTTCGATATTCCGCCCCCTGCGTTTGCAGGATCTcg                                                                          >  2:2643697/3‑100 (MQ=255)
                    gAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTcg                                                                          >  2:1261754/1‑100 (MQ=255)
                    gAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTcg                                                                          <  1:1109536/100‑1 (MQ=255)
                    gAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGCTTGCAGGATTTcg                                                                          >  2:2765136/1‑100 (MQ=255)
                      aaaGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTcgcg                                                                        >  1:2947488/1‑100 (MQ=255)
                         ggCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGCTTTCGCGCaa                                                                     >  2:1134740/1‑100 (MQ=255)
                         ggCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCaa                                                                     >  1:393793/1‑100 (MQ=255)
                         ggCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCaa                                                                     >  1:757240/1‑100 (MQ=255)
                            cAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACg                                                                  >  1:208557/1‑100 (MQ=255)
                               ccTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGgcg                                                               >  2:177010/1‑100 (MQ=255)
                                cTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGgcgc                                                              >  2:940867/1‑100 (MQ=255)
                                cTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGgcgc                                                              >  1:291252/1‑100 (MQ=255)
                                cTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGgcgc                                                              <  2:2445032/100‑1 (MQ=255)
                                   gCTTGAGACTGGCCTTTCTGACAGATTCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCTCAAACGGCGCAgg                                                           <  2:2021269/100‑1 (MQ=255)
                                   gCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAgg                                                           <  2:407303/100‑1 (MQ=255)
                                   gCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAgg                                                           >  1:621147/1‑100 (MQ=255)
                                    cTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCAGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGc                                                          <  1:188210/100‑1 (MQ=255)
                                     ttGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGcc                                                         >  1:1094985/1‑100 (MQ=255)
                                     ttGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGcc                                                         <  1:449486/100‑1 (MQ=255)
                                      tGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCt                                                        <  1:2702397/100‑1 (MQ=255)
                                       gagaCTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCtt                                                       <  2:3078269/100‑1 (MQ=255)
                                       gagaCTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCtt                                                       <  2:1139375/100‑1 (MQ=255)
                                        agaCTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTc                                                      >  1:69807/1‑100 (MQ=255)
                                         gaCTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCa                                                     <  2:2302053/100‑1 (MQ=255)
                                           cTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAAc                                                   <  1:2755551/100‑1 (MQ=255)
                                           cTGGCCTTTCTGACAGATGCGTACTTACTCGCGGAACAGCGCTTCGATATTCAGCACCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAAc                                                   >  1:1606084/1‑100 (MQ=255)
                                            tGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAAcc                                                  <  1:2680511/100‑1 (MQ=255)
                                             ggCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCt                                                 >  1:2852072/1‑100 (MQ=255)
                                             ggCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCt                                                 >  1:3081735/1‑100 (MQ=255)
                                              gCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTg                                                >  2:736903/1‑100 (MQ=255)
                                               ccTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGa                                               >  2:2885140/1‑100 (MQ=255)
                                                 tttCTGACAGATGCTTACTTACTCGCGGCACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCCAACGGCGCAGGCCTTCAACCTGaa                                              <  2:641461/99‑1 (MQ=255)
                                                 tttCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAAt                                             <  1:1261768/100‑1 (MQ=255)
                                                 tttCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAAt                                             >  1:2630292/1‑100 (MQ=255)
                                                    cTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTg                                          >  2:2612515/1‑100 (MQ=255)
                                                     tGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCACCCTGAATCTgg                                         <  2:3174591/100‑1 (MQ=255)
                                                      gACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGc                                        >  2:840568/1‑100 (MQ=255)
                                                      gACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGc                                        <  2:289064/100‑1 (MQ=255)
                                                      gACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGc                                        <  1:2249455/100‑1 (MQ=255)
                                                      gACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGc                                        <  2:2064600/100‑1 (MQ=255)
                                                       aCAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCg                                       <  2:2707873/100‑1 (MQ=255)
                                                        cAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGa                                      <  2:69805/100‑1 (MQ=255)
                                                         aGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGaa                                     <  1:4829/100‑1 (MQ=255)
                                                         aGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGaa                                     <  2:535964/100‑1 (MQ=255)
                                                            tGCTTACTTACTCGCGGAACAGCGCTTCGATATACAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCACAGGCCTTCAACCTGAATCTGGCGAacac                                  >  2:1959983/1‑100 (MQ=255)
                                                             gCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACg                                 <  1:53573/100‑1 (MQ=255)
                                                              cttacttaCTCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGt                                <  1:2629143/100‑1 (MQ=255)
                                                                tacttactCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGCTTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTc                              >  2:1080746/1‑100 (MQ=255)
                                                                   ttactCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACg                           >  1:911676/1‑100 (MQ=255)
                                                                     actCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACggg                         <  1:840578/100‑1 (MQ=255)
                                                                     actCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACggg                         <  1:2254999/100‑1 (MQ=255)
                                                                      ctCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCGGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGt                        >  1:840250/1‑100 (MQ=255)
                                                                      ctCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGt                        >  1:1431405/1‑100 (MQ=255)
                                                                      ctCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGt                        >  1:2978987/1‑100 (MQ=255)
                                                                       tCGCGGAACAGCGCTTCGATATTCAGCCCCCTGCCTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACATGTTCACGGGTg                       <  2:1720442/100‑1 (MQ=255)
                                                                        cgcgGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGa                      >  1:3257517/1‑100 (MQ=255)
                                                                        cgcgGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGa                      <  1:1859399/100‑1 (MQ=255)
                                                                         gcgGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAg                     <  2:2914354/100‑1 (MQ=255)
                                                                          cgGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAgg                    <  1:764044/100‑1 (MQ=255)
                                                                          cgGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAgg                    <  2:125498/100‑1 (MQ=255)
                                                                          cgGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAg                     >  2:2205100/1‑99 (MQ=255)
                                                                          cgGAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAg                     <  2:1108876/99‑1 (MQ=255)
                                                                           ggAACAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGc                   <  1:2321169/100‑1 (MQ=255)
                                                                              aCAGCGCTTCGATATTTAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCaa                >  2:1862109/1‑100 (MQ=255)
                                                                              aCAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCaa                >  2:3075117/1‑100 (MQ=255)
                                                                               cAGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAAt               <  1:1207263/100‑1 (MQ=255)
                                                                                aGCGCTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAAtt              >  2:706183/1‑100 (MQ=255)
                                                                                   gcTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCa           >  2:1317373/1‑100 (MQ=255)
                                                                                   gcTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCa           <  1:824728/100‑1 (MQ=255)
                                                                                    cTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCAc          <  2:2947453/100‑1 (MQ=255)
                                                                                    cTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCAc          <  1:1960693/100‑1 (MQ=255)
                                                                                    cTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCAc          >  2:1361244/1‑100 (MQ=255)
                                                                                    cTTCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCAc          >  1:1460563/1‑100 (MQ=255)
                                                                                     ttCGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACg         >  2:3024651/1‑100 (MQ=255)
                                                                                       cGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGAc       <  1:1015322/100‑1 (MQ=255)
                                                                                       cGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGAc       <  2:46338/100‑1 (MQ=255)
                                                                                       cGATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGAc       <  1:527638/100‑1 (MQ=255)
                                                                                        gATATTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGAcc      <  2:208552/100‑1 (MQ=255)
                                                                                          tatTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTa    <  2:2103155/100‑1 (MQ=255)
                                                                                           atTCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTAc   <  2:3029172/100‑1 (MQ=255)
                                                                                            ttCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACa  <  2:393786/100‑1 (MQ=255)
                                                                                            ttCAGCCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACa  <  1:238290/100‑1 (MQ=255)
                                                                                                 |                                                                                               
CGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACA  >  NC_000913/2866492‑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: