Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,248,693 T→A 5.6% intergenic (+27/+72) tsaA → / ← rnt peroxiredoxin/ribonuclease T

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,248,6930TA 5.6% 1086.0 / 39.9 355intergenic (+27/+72)tsaA/rntperoxiredoxin/ribonuclease T
Reads supporting (aligned to +/- strand):  ref base T (166/160);  new base A (6/14);  total (174/181)
Fisher's exact test for biased strand distribution p-value = 1.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.46e-01

Only 100 of 361 total aligned reads displayed.

CTGCAGTTCCACGAAGAGCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGCTCAGCGGTTTTCTC  >  NC_002947/1248544‑1248842
                                                                                                                                                     |                                                                                                                                                     
  gCAGTTCCACGAAGAGCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTc                                                                                                                                                     <  1:1356426/150‑1 (MQ=255)
     gTTCCACGAAGAGCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTg                                                                                                                                                  <  2:3244463/150‑1 (MQ=255)
        ccACGAAGAGCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGcc                                                                                                                                               >  1:2598669/1‑150 (MQ=255)
          aCGAAGAGCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCgg                                                                                                                                             >  1:789703/1‑150 (MQ=255)
             aaGAGCACGGCGAAGTGTGCCCTGCCAACTGGAAAATAGGCGACAAGGGCATGAACGCTTCGCCTGAAGTCGTTGATGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAAACTGCCTTCGTGGCCGGttt                                                                                                                                          >  1:6579728/1‑150 (MQ=255)
             aaGAGCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGttt                                                                                                                                          >  2:536565/1‑150 (MQ=255)
             aaGAGCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCGTCGTGGCCGGttt                                                                                                                                          >  2:4838096/1‑150 (MQ=255)
              agagCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGtttt                                                                                                                                         <  1:3742242/150‑1 (MQ=255)
                agCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATAAAAAAACCGGCCACGAAGGCCGGtttttt                                                                                                                                       <  1:1883556/150‑1 (MQ=255)
                 gCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGttttttt                                                                                                                                      <  1:5851848/150‑1 (MQ=255)
                  cACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTa                                                                                                                                     <  1:6433893/150‑1 (MQ=255)
                   aCGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTAt                                                                                                                                    <  1:6957602/150‑1 (MQ=255)
                   aCGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTAt                                                                                                                                    <  2:4417801/150‑1 (MQ=255)
                   aCGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTAt                                                                                                                                    <  2:1949770/150‑1 (MQ=255)
                    cGGCGAAGTGTGCGCTGCCAACTGGAAAAAAGGCGACAAGGGCATGCACGCTTCGCCTGAAGGCGTTGCTGCCTCCCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATAAAAAAACCGGCCACGAAGGCCGGTTTTTTCATg                                                                                                                                   <  2:6134462/150‑1 (MQ=38)
                       cGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCt                                                                                                                                >  2:1318004/1‑150 (MQ=255)
                         aaGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTc                                                                                                                              <  1:2431012/150‑1 (MQ=255)
                               gCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGc                                                                                                                        <  1:5256517/150‑1 (MQ=255)
                                cccTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGcc                                                                                                                       <  2:2104201/150‑1 (MQ=255)
                                  cTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCgg                                                                                                                     >  2:5008491/1‑150 (MQ=255)
                                   ggCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGACCGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGc                                                                                                                    <  2:577734/149‑1 (MQ=255)
                                       aaCTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGc                                                                                                                    <  2:2657581/146‑1 (MQ=255)
                                         cTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTc                                                                                                              <  2:2734892/150‑1 (MQ=255)
                                          tGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTcg                                                                                                             >  2:4110257/1‑150 (MQ=255)
                                          tGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTcg                                                                                                             >  1:6159347/1‑150 (MQ=255)
                                            gAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTcgcg                                                                                                           >  2:6627108/1‑150 (MQ=255)
                                            gAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTcgcg                                                                                                           <  2:2531999/150‑1 (MQ=255)
                                             aaaaaaGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCgg                                                                                                          <  1:5806024/150‑1 (MQ=255)
                                             aaaaaaGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCgg                                                                                                          <  2:7054384/150‑1 (MQ=255)
                                               aaaaGGCGTCAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTTCCTGAGCGAGAACGCTGGCAAGCTGTATTTTCTGAGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGgc                                                                                                        <  1:3475745/150‑1 (MQ=255)
                                               aaaaGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGgc                                                                                                        >  1:6936299/1‑150 (MQ=255)
                                                 aaGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGgcgc                                                                                                      >  1:2926638/1‑150 (MQ=255)
                                                 aaGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGgcgc                                                                                                      >  1:6672579/1‑150 (MQ=255)
                                                      gACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCcg                                                                                                 >  2:1264011/1‑150 (MQ=255)
                                                      gACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCcg                                                                                                 <  2:2287785/150‑1 (MQ=255)
                                                      gACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCcg                                                                                                 >  2:6704794/1‑150 (MQ=255)
                                                      gACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTGTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCcg                                                                                                 >  2:6555113/1‑150 (MQ=255)
                                                       ccacgggCATGCACGCTTCGCCTGACGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCACGCTGGAAGTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTGCGCGGGCGCGCCcgc                                                                                                <  2:1767698/146‑1 (MQ=255)
                                                       aCAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCcgc                                                                                                >  2:4554553/1‑150 (MQ=255)
                                                         aaGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGa                                                                                              >  2:2848300/1‑150 (MQ=255)
                                                          aGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGaa                                                                                             <  2:3074702/150‑1 (MQ=255)
                                                           gggCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGaag                                                                                            <  1:2265225/150‑2 (MQ=255)
                                                           gggCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAAt                                                                                            >  1:4037986/1‑150 (MQ=255)
                                                               aTGAACGCTTCGCCTGATGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGATTTATTTGCTGCGTGAATTAAAAAACCGGCCTTGGTGGCCGGTTTTTTTTTGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGaagaggc                                                                                        <  1:7198817/150‑6 (MQ=255)
                                                               aTGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTACTTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCGGCGCGGGGGGGCCcccgccgaggg                                                                                        >  1:2700569/1‑143 (MQ=255)
                                                                 gAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCttt                                                                                      >  1:1176694/1‑150 (MQ=255)
                                                                  aaCGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTa                                                                                     <  1:6109367/150‑1 (MQ=255)
                                                                   aCGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAg                                                                                    <  2:5222281/150‑1 (MQ=255)
                                                                    cGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTatc                                                                                   >  1:6011279/1‑148 (MQ=255)
                                                                     gCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCTTGAAAATACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCa                                                                                  >  2:2441569/1‑150 (MQ=255)
                                                                       ttCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGa                                                                                <  2:3496501/150‑1 (MQ=255)
                                                                       ttCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGa                                                                                <  2:119215/150‑1 (MQ=255)
                                                                       gtCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTGGTGGCCGGTTTTTTTATGCCTTCCCGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGa                                                                                <  2:6845370/149‑1 (MQ=255)
                                                                        tCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTTCTGCGTGCATGAAAAAACCGGCCGTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGAt                                                                               <  1:2127200/150‑1 (MQ=255)
                                                                           ccTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGCATGCCGTTAGCAGATCaa                                                                            >  1:1015484/1‑150 (MQ=255)
                                                                             tGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATTCCTTTAGCAGATCAATc                                                                          >  2:233534/1‑150 (MQ=255)
                                                                             tGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATc                                                                          >  2:6433355/1‑150 (MQ=255)
                                                                             tGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTGTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATc                                                                          >  2:3947190/1‑150 (MQ=255)
                                                                                aGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGtt                                                                       <  2:5301073/150‑1 (MQ=255)
                                                                                     ttGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGGGCGCCCGCGAATGCCGTTAGCAGATCAATCGGTGTAAt                                                                  >  1:742036/1‑150 (MQ=255)
                                                                                        ctgcCTACCTGAGCGAGATCGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATcgc                                                               <  2:3355879/150‑1 (MQ=255)
                                                                                        ctgcCTACCTGAGCGAGAACGCTGGCAAGCTGTAAGTGCTGCGTGCATGAAAAAACCGGCCGTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATcgc                                                               <  2:3803190/150‑1 (MQ=255)
                                                                                         tgcCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATcgcg                                                              >  2:7361351/1‑150 (MQ=255)
                                                                                          gcCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCgcgc                                                             <  2:2249985/150‑1 (MQ=255)
                                                                                                 tGAGCGAGAACGCTGGCATGCTGTAATTGCTGCGTGCATGCAAACACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTGTAGCAGATCATTCTTTGAAATCGCGCCAgccga                                                      <  2:5797578/150‑2 (MQ=255)
                                                                                                  gAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAgccgcc                                                     <  1:114246/150‑1 (MQ=255)
                                                                                                     cGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCAt                                                  >  2:1032698/1‑150 (MQ=255)
                                                                                                      gagaACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTGTAGCAGATCAATCGTTGAAATCGCGCCACCCGccccat                                                 >  1:6876610/1‑147 (MQ=255)
                                                                                                       agaACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCAttt                                                <  2:840494/150‑1 (MQ=255)
                                                                                                        gaACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCAtttc                                               >  2:464146/1‑150 (MQ=255)
                                                                                                            gCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCAtttctttc                                           <  1:1643289/150‑1 (MQ=255)
                                                                                                            gCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCAtttctttc                                           >  2:4335624/1‑150 (MQ=255)
                                                                                                             cTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTcc                                          <  1:3971374/150‑1 (MQ=255)
                                                                                                                gCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCAGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCACCAGCCGCCCATTTCTTTCCAGc                                       <  2:2014652/150‑1 (MQ=255)
                                                                                                                  aaGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCAATTCGTTCCAGCgg                                     >  2:6305357/1‑150 (MQ=255)
                                                                                                                    gCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGGGCGCCCGCGACTGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGtt                                   >  2:5648692/1‑150 (MQ=255)
                                                                                                                    gCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGtt                                   <  2:1733401/150‑1 (MQ=255)
                                                                                                                      tGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTACAGCGGTTGa                                 >  1:1059515/1‑150 (MQ=255)
                                                                                                                       gTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTATCAGAAAAATAGTTGAAATCGCGCCAGCCGCCCAATTGTTTCCAGCGGTTGAc                                >  1:5657455/1‑150 (MQ=255)
                                                                                                                            tgctgcGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGc                           <  1:4254188/150‑1 (MQ=255)
                                                                                                                               tgcGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGCATGCCGTTAGCAGATCAAACGTTTGAATCGCGCCAGCCGCCCATGTCTTTACAGCGGTTGGCGATGCCGc                        >  1:4538639/1‑150 (MQ=255)
                                                                                                                               tgcGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCAtttctttc                                           <  1:3299678/131‑1 (MQ=255)
                                                                                                                                gcGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTGTGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGTAATCGCGGCCGCCGCCCCTTTCGTTGCGGCGGTTGACGGTggggcg                       >  1:1252372/1‑145 (MQ=255)
                                                                                                                                gcGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCa                       <  1:2091898/150‑1 (MQ=255)
                                                                                                                                gcGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCa                       <  1:5036305/150‑1 (MQ=255)
                                                                                                                                      aTGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAg                 >  2:5892026/1‑150 (MQ=255)
                                                                                                                                      aTGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAg                 <  1:4812610/150‑1 (MQ=255)
                                                                                                                                      aTGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAAACGATGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAg                 >  2:7206752/1‑150 (MQ=255)
                                                                                                                                        gAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGTGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGCt               >  1:7350844/1‑150 (MQ=255)
                                                                                                                                        gAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGCt               >  2:1546120/1‑150 (MQ=255)
                                                                                                                                         aaaaaaCCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCCGGCGCGCCCGCCCATGCCGTTTGCAGATCAATCGGTTTAATCGCGCCAGCCGCCCCAATCTTTACAGCGGTTGAAGGTGCCGCAGAACAGCtg              >  1:6274536/1‑149 (MQ=255)
                                                                                                                                         aaaaaaCCGGCCACGAAGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAGTCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGGGGATGCCGCAGAACAGCTc              >  1:978154/1‑150 (MQ=255)
                                                                                                                                          aaaaaCCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGCTCa             >  2:170786/1‑150 (MQ=255)
                                                                                                                                          aaaaaCCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGCTCa             <  2:5962561/150‑1 (MQ=255)
                                                                                                                                          aaaaaCCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGCTCa             >  2:4679059/1‑150 (MQ=255)
                                                                                                                                           aaaaCCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGCTCAg            >  2:2982324/1‑150 (MQ=255)
                                                                                                                                            aaaCCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGCTCAGc           <  1:1522799/150‑1 (MQ=255)
                                                                                                                                                cGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCTCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGCTCAGCTGtt       <  2:2926743/150‑1 (MQ=255)
                                                                                                                                                cGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGCTCAGCGGtt       <  2:5085587/150‑1 (MQ=255)
                                                                                                                                                cGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAAAAGCGCCAGCCGCCCAAATCTTTCCAGCGGTTGACGATGCCGCAGAACCGCACAGCGGtt       >  1:3032834/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CTGCAGTTCCACGAAGAGCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGCTCAGCGGTTTTCTC  >  NC_002947/1248544‑1248842

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: