Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,560,381 G→T 6.8% intergenic (+28/+46) sbcB → / ← PP_1366 exodeoxyribonuclease I/transcriptional regulator MvaT

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,560,3810GT 6.8% 958.5 / 36.6 284intergenic (+28/+46)sbcB/PP_1366exodeoxyribonuclease I/transcriptional regulator MvaT
Reads supporting (aligned to +/- strand):  ref base G (144/120);  new base T (14/6);  total (158/126)
Fisher's exact test for biased strand distribution p-value = 2.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 286 total aligned reads displayed.

GGGCGACTTCGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCA  >  NC_002947/1560233‑1560529
                                                                                                                                                    |                                                                                                                                                     
gggCGACTTCGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAAtgt                                                                                                                                                      <  2:1692496/150‑1 (MQ=255)
    gACTTCGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCg                                                                                                                                                  >  1:5823139/1‑150 (MQ=255)
       ttCGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCg                                                                                                                                               >  2:1987235/1‑150 (MQ=255)
       ttCGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCg                                                                                                                                               >  2:4163419/1‑150 (MQ=255)
       ttCGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCg                                                                                                                                               >  1:2520181/1‑150 (MQ=255)
        tCGAAACGGCGCGCCAGCACGCCTGGGAGGGCGCCGACGTGGCAGGCTGGCGCGTGTTGGCGGCCTGGCAACTGCATGCTCGACAATGGCAGGAACGCTTTTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGt                                                                                                                                              <  1:5650022/150‑1 (MQ=255)
         cGAAAAGGCGCGCCAGCACGCCTGGGAGGGCGCCGACGCGGCAGGCAGGCGCGTGTTGGCGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGAGGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGtt                                                                                                                                             <  1:992866/150‑1 (MQ=255)
         cGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGCGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGAGGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGtt                                                                                                                                             <  2:1785354/150‑1 (MQ=255)
          gAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGTGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGttt                                                                                                                                            <  2:2307244/150‑1 (MQ=255)
          gAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGttt                                                                                                                                            >  1:1022036/1‑150 (MQ=255)
          gAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGttg                                                                                                                                            >  1:6520416/1‑149 (MQ=255)
           aaaaGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGtttt                                                                                                                                           >  2:3808940/1‑150 (MQ=255)
            aaaGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCACATTGCCGGCGttttt                                                                                                                                          <  2:2420180/150‑1 (MQ=255)
            aaaGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGttttt                                                                                                                                          >  2:3164817/1‑150 (MQ=255)
                 cgcgcCAGCACGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAACACGCCGGCAATGTGCCGGCGTTTTtgttgt                                                                                                                                     <  2:5646267/150‑1 (MQ=255)
                  gcgcCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTtgtg                                                                                                                                    >  1:4645390/1‑150 (MQ=255)
                  gcgcCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTtgtg                                                                                                                                    >  1:4116556/1‑150 (MQ=255)
                    gcCAGCACGCCGGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTGGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCg                                                                                                                                  <  2:954554/150‑1 (MQ=255)
                    gcCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCg                                                                                                                                  <  2:4682511/150‑1 (MQ=255)
                      cAGCACGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCCGGCTCGACAATTGCAGGAACGCTATTCCACCGGCTGATGGCAGATAAACAAACACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGaa                                                                                                                                <  1:174934/150‑1 (MQ=255)
                      cAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGaa                                                                                                                                <  2:4454311/150‑1 (MQ=255)
                          aaGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCaa                                                                                                                            >  1:776963/1‑150 (MQ=255)
                          aaGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCaa                                                                                                                            >  1:3790398/1‑150 (MQ=255)
                           aGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCCCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAag                                                                                                                           >  1:1712405/1‑150 (MQ=255)
                           aGCCGGGGTGGGCGCCGACGCGGCAGGCAGGGGCGTGTGGGCGGCCTGGAAACTGCTTGCTCGACAATTGCAGTAACGCTTTACAATCGGCTGTTGGATTATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAag                                                                                                                           <  1:2592453/150‑1 (MQ=255)
                              cTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAgg                                                                                                                        <  2:1745035/150‑1 (MQ=255)
                               tGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAggg                                                                                                                       >  1:5358097/1‑150 (MQ=255)
                                gggAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAACACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGc                                                                                                                      <  1:432148/150‑1 (MQ=255)
                                    gggCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTAt                                                                                                                  <  2:4019693/150‑1 (MQ=255)
                                      gcgcCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGc                                                                                                                <  1:7404644/150‑1 (MQ=255)
                                        gcCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGAAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAATCGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAATCAAGAGGGCGTATGCgg                                                                                                              >  2:6763946/1‑150 (MQ=255)
                                         cccacgaGGCAGGCAGGCGCGTGTTGGAGGCCGGCCAACTGCATGCTCGACAATTGCAGGAACGCTATGCAATCGGCTGAGGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGa                                                                                                             <  1:3880456/147‑1 (MQ=255)
                                           gacgaGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATa                                                                                                           >  1:3606448/1‑150 (MQ=255)
                                            acgaGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATaa                                                                                                          >  2:687726/1‑150 (MQ=255)
                                            acgaGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATaa                                                                                                          >  2:5725845/1‑150 (MQ=255)
                                              gaGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATc                                                                                                        >  2:2667847/1‑150 (MQ=255)
                                               aggcaggcaggcGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGACCGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCa                                                                                                       <  1:2717676/150‑1 (MQ=255)
                                               aggcaggcaggcGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCa                                                                                                       <  2:3828329/150‑1 (MQ=255)
                                               aggcaggcaggcGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGATATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCa                                                                                                       <  1:2284370/150‑1 (MQ=255)
                                                  caggctggcGGGTGTTGGAGGCCTGGCAACTGCATGCTCGACATGTGCAGGCAGGCTATGCACTCGGCTGATGGCAGATAAACAAAAAGGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGTGGGCGTATGCGGATAATCAGTc                                                                                                    <  2:5297340/150‑1 (MQ=255)
                                                  caggcaggcGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGTAACGCTATTAAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTc                                                                                                    <  1:4998517/150‑1 (MQ=255)
                                                  caggcaggcGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAGCGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTc                                                                                                    <  2:3384064/150‑1 (MQ=255)
                                                    ggcaggcGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCca                                                                                                  >  1:2170781/1‑150 (MQ=255)
                                                     gcaggcGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCcag                                                                                                 >  2:25352/1‑150 (MQ=255)
                                                       aggcGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCcagca                                                                                               >  1:4349410/1‑150 (MQ=255)
                                                         gcgcGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAgg                                                                                             >  1:4786492/1‑150 (MQ=255)
                                                          cgcgTTTTGGAGGCCTGGCAACTGCATGCCCGACAATTGCAGGAACGCTATACAATCGGCTGATGGCAGAAAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAggg                                                                                            <  1:133501/150‑1 (MQ=255)
                                                           gcgTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCcagctcagc                                                                                           >  1:2209051/1‑145 (MQ=255)
                                                               gtTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGcc                                                                                       >  2:4071331/1‑150 (MQ=255)
                                                                  ggAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCTCAg                                                                                    <  1:573707/150‑1 (MQ=255)
                                                                   gAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGc                                                                                   >  1:1257719/1‑150 (MQ=255)
                                                                     ggCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCATATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCtt                                                                                 <  1:6899766/150‑1 (MQ=255)
                                                                       ccTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTc                                                                               >  1:1461229/1‑150 (MQ=255)
                                                                        cTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCACTGTGCCGGCGTTTTTGTTGTGCGAAGCAAGTGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCTGCTTTct                                                                              >  2:7122226/1‑149 (MQ=255)
                                                                         tGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGCTAATCAGTCCAGCAGGGTTGCCCAGCTTTCaa                                                                             >  2:6583136/1‑150 (MQ=255)
                                                                         tGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCaa                                                                             >  2:1248424/1‑150 (MQ=255)
                                                                         tGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCaa                                                                             >  2:6687049/1‑150 (MQ=255)
                                                                             aaCTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCCGGGGTGCCCAGCTTTCAACCAc                                                                         >  1:5074968/1‑150 (MQ=255)
                                                                               cTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGt                                                                       >  2:1885858/1‑150 (MQ=255)
                                                                                  cATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGcc                                                                                       >  2:4330128/1‑131 (MQ=255)
                                                                                         gACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCa                                                            >  1:3645745/1‑150 (MQ=255)
                                                                                          aCAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCaa                                                           <  2:730151/150‑1 (MQ=255)
                                                                                             aTTGCAGGAACGATATTCAATCGGCTGTTGGCATATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTg                                                        <  1:3394800/150‑1 (MQ=255)
                                                                                              ttGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAg                                                                                                                               <  2:7172345/79‑1 (MQ=255)
                                                                                              ttGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTgg                                                       <  2:3663785/150‑1 (MQ=255)
                                                                                                 cAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCtt                                                    <  2:6279792/150‑1 (MQ=255)
                                                                                                  aGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCGTTCAACCACGTCGCCGCCCCAAGTGGCttt                                                   >  2:7097575/1‑150 (MQ=255)
                                                                                                  aGGAACGCGATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCttt                                                   <  1:7042634/150‑1 (MQ=255)
                                                                                                    gAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTcc                                                 >  2:5006153/1‑150 (MQ=255)
                                                                                                     aaCGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCa                                                <  2:684963/150‑1 (MQ=255)
                                                                                                        gCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCActc                                             >  1:1633997/1‑150 (MQ=255)
                                                                                                          tATTCAATCGGCTGATGGCAGATAAAGAAAAACGCCGGCAATGTGCCGGCGTTTGGGTTGTGCGAAGCAAGTGGGCGTATGGGGATAATCAGTCCAGCAGGGTGGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCtt                                           >  1:2960119/1‑150 (MQ=255)
                                                                                                          tATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCCACGTGGCTTTCCActct                                            >  2:6101603/1‑150 (MQ=255)
                                                                                                          tATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCTACTCtt                                           >  1:6073299/1‑150 (MQ=255)
                                                                                                          tATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCtt                                           >  2:2076372/1‑150 (MQ=255)
                                                                                                            ttCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTc                                         >  2:1269176/1‑150 (MQ=255)
                                                                                                             tCATTCGTCTGAGGGCATATAAACAAAAACGCCGGCAATGTGCCGGCTTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCa                                        <  1:5006381/150‑1 (MQ=255)
                                                                                                             tCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCa                                        <  2:2755104/150‑1 (MQ=255)
                                                                                                             tCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCa                                        <  2:6503624/150‑1 (MQ=255)
                                                                                                                    gCTGAGGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCAAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCtt                                 <  2:1423739/150‑1 (MQ=255)
                                                                                                                      tGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCGACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTtgt                               >  1:1131945/1‑150 (MQ=255)
                                                                                                                       taTGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTtgtg                              <  1:5581299/149‑1 (MQ=255)
                                                                                                                       gATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCAGGGTCTtgtg                              <  1:4478240/150‑1 (MQ=255)
                                                                                                                             aGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATTCGGCTAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCGTTCAGGGTCTTGTGGTTgcc                        >  1:3284194/1‑150 (MQ=255)
                                                                                                                             aGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTTCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTCgcc                        >  1:5167272/1‑150 (MQ=255)
                                                                                                                             aGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTgcc                        >  1:5520843/1‑150 (MQ=255)
                                                                                                                               cacaaCAAAAACGCCGGCACATTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTCCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTgccgc                      >  2:4833124/4‑150 (MQ=39)
                                                                                                                                   aCAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCAGGGTCTTGGGGTTGCCGCCttt                  >  1:1577798/1‑150 (MQ=255)
                                                                                                                                     aaaaaCGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCca                                                                                                  <  1:3971824/69‑1 (MQ=255)
                                                                                                                                     aaaaaCGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTggg                >  2:6200731/1‑150 (MQ=255)
                                                                                                                                         aCGCCGGCACATTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGTGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTCGCCGCCTTTGGTTTc            >  1:7174995/1‑150 (MQ=255)
                                                                                                                                         aCGCCGGCACATTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGttt             >  2:3236007/1‑150 (MQ=255)
                                                                                                                                         aCGCCGGCACATTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTc            >  2:391427/1‑150 (MQ=255)
                                                                                                                                         aCGCCGGCACATTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTATTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTc            >  2:6311826/1‑150 (MQ=255)
                                                                                                                                         aCGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTTCACTCTTTCAGGGTCTTGTGGTTGCCGCCGTTGGTTTc            >  1:7134309/1‑150 (MQ=255)
                                                                                                                                         aCGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTc            <  2:5111423/150‑1 (MQ=255)
                                                                                                                                         aCGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGGTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTc            >  1:5502034/1‑150 (MQ=255)
                                                                                                                                         aCGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGACCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTc            <  2:7018387/150‑1 (MQ=255)
                                                                                                                                           gCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGa          >  2:5841094/1‑150 (MQ=255)
                                                                                                                                           gCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTCCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGa          <  1:1885858/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GGGCGACTTCGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCAATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTTCGATCACTTCA  >  NC_002947/1560233‑1560529

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