Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,243,746 C→G 7.7% intergenic (+29/+70) tsaA → / ← rnt peroxiredoxin/ribonuclease T

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,243,7460CG 7.7% 578.7 / 27.9 180intergenic (+29/+70)tsaA/rntperoxiredoxin/ribonuclease T
Reads supporting (aligned to +/- strand):  ref base C (76/90);  new base G (7/7);  total (83/97)
Fisher's exact test for biased strand distribution p-value = 7.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

Only 100 of 184 total aligned reads displayed.

GCAGTTCCACGAAGAGCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGCTCAGCGGTTTTCTCG  >  NC_002947/1243597‑1243894
                                                                                                                                                     |                                                                                                                                                    
gCAGTTCCACGAAGAGCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTc                                                                                                                                                      >  2:5070767/1‑150 (MQ=255)
  aGTTCCACGAAGAGCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGt                                                                                                                                                    >  1:3825752/1‑150 (MQ=255)
  aGTTCCACGAAGAGCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGt                                                                                                                                                    <  2:566752/150‑1 (MQ=255)
   gTTCCACGAAGAGCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTg                                                                                                                                                   <  1:1937052/150‑1 (MQ=255)
       cACGAAGAGCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCg                                                                                                                                               <  2:2502810/150‑1 (MQ=255)
        aCGAAGAGCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCgg                                                                                                                                              >  1:3352347/1‑150 (MQ=255)
        aCGAAGAGCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCgg                                                                                                                                              >  2:135243/1‑150 (MQ=255)
           aaGAGCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGttt                                                                                                                                           <  1:135243/150‑1 (MQ=255)
            agagCACGGCGACGTGGGCCCGGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTCCCTGAGCGAGAACGCTGGCAAGCTGTAAGTGCTGCGTGCATGAAAAAACCGGCCACGACGGCCGGtttt                                                                                                                                          <  2:202128/150‑1 (MQ=255)
            agagCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGtttt                                                                                                                                          >  1:1251690/1‑150 (MQ=255)
             gagCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGttttt                                                                                                                                         <  1:1456748/150‑1 (MQ=255)
             gagCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCACGAAGGCCGGttttt                                                                                                                                         <  2:4730396/150‑1 (MQ=255)
              agCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGtttttt                                                                                                                                        <  2:83812/150‑1 (MQ=255)
               gCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGttttttt                                                                                                                                       >  1:4877795/1‑150 (MQ=255)
                cACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCACGAAGGCCGGttttttca                                                                                                                                      <  1:4011535/150‑3 (MQ=255)
                cACGGCGAAGTGGGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATAAAAAAACCGGCCACGAAGGCCGGttttttca                                                                                                                                      <  2:3242269/150‑3 (MQ=38)
                 aCGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTAt                                                                                                                                     <  1:5354104/150‑1 (MQ=255)
                   ggCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGc                                                                                                                                   >  2:3732948/1‑150 (MQ=255)
                   ggCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGc                                                                                                                                   <  2:3825752/150‑1 (MQ=255)
                    gCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGcc                                                                                                                                  <  2:4877795/150‑1 (MQ=255)
                      gAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCtg                                                                                                                                >  1:5160247/1‑149 (MQ=255)
                       taGTGTGCCCCGCCCACTGGAAAAAAGGCGCCCAGGGCATTAACGCGTCGCCTGCAGGCGGTGCGGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCACGAAGGCCGGTTTTTTCATGcacgc                                                                                                                               <  2:3321296/149‑5 (MQ=25)
                        aGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGGGCATAAAAAAACCGGCCACGAAGGCCGGTTTTTTCATGcacgca                                                                                                                              <  2:5131148/150‑6 (MQ=37)
                         gtgtGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCt                                                                                                                             <  2:4445179/150‑1 (MQ=255)
                             gCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGc                                                                                                                         >  1:2957160/1‑150 (MQ=255)
                              cccTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGcc                                                                                                                        <  1:3948621/150‑1 (MQ=255)
                              cccTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGcc                                                                                                                        <  2:2113020/150‑1 (MQ=255)
                                 tGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGc                                                                                                                     >  2:4472307/1‑150 (MQ=255)
                                        tGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGTCCTCTTcg                                                                                                              >  2:825786/1‑150 (MQ=255)
                                         ggAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTcgc                                                                                                             <  2:3917159/150‑1 (MQ=255)
                                          gAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTcgcg                                                                                                            >  1:596503/1‑150 (MQ=255)
                                           aaaaaaGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCgg                                                                                                           <  2:3093850/150‑1 (MQ=255)
                                            aaaaaGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCggg                                                                                                          >  2:47369/1‑150 (MQ=255)
                                            aaaaaGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCggg                                                                                                          >  2:4840861/1‑150 (MQ=255)
                                              aaaGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGgcg                                                                                                        >  2:5041487/1‑150 (MQ=255)
                                               aaGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGgcgc                                                                                                       >  1:846298/1‑150 (MQ=255)
                                                aGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGgcgcg                                                                                                      >  1:4361012/1‑150 (MQ=255)
                                                 ggCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGgcgcgc                                                                                                     <  2:3162598/150‑1 (MQ=255)
                                                     aCAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCcgc                                                                                                 <  2:1701917/150‑1 (MQ=255)
                                                       aaGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGa                                                                                               <  1:3481192/150‑1 (MQ=255)
                                                         gggCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGaag                                                                                             <  1:4840861/150‑2 (MQ=255)
                                                          ggCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTATTCGCGGGCGCGCCCGCGaaga                                                                                            <  2:308335/150‑3 (MQ=255)
                                                           gCATGAACGCTTCGCCTGATGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGGTGCGTGAATGAAAAAACCGGCCTTCTTGGCGGGTTTTTTTATGCCTTCCTGTGCCGGCCTATTCGCGGGCGCGCCCGCGaagag                                                                                           <  1:4111022/150‑4 (MQ=255)
                                                              tGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCtt                                                                                        <  1:3086096/150‑1 (MQ=255)
                                                              tGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCtt                                                                                        <  1:2342752/150‑1 (MQ=255)
                                                               gAACGCTTCGCCTGAAGGCGTTGCTGCATACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACAGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCCTCTGCCttt                                                                                       >  1:3693190/1‑150 (MQ=255)
                                                                 aCGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAg                                                                                     <  2:2689963/150‑1 (MQ=255)
                                                                 aCGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAg                                                                                     <  1:5070767/150‑1 (MQ=255)
                                                                   gCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCCAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCa                                                                                   <  2:3559129/150‑1 (MQ=255)
                                                                   gCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCa                                                                                   <  2:756547/150‑1 (MQ=255)
                                                                     ttCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGa                                                                                 >  2:2949005/1‑150 (MQ=255)
                                                                     ttCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGa                                                                                 >  2:732992/1‑150 (MQ=255)
                                                                     ttCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGa                                                                                 <  2:2329914/150‑1 (MQ=255)
                                                                           tGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTGTAGCAGATCAATc                                                                           >  1:3142358/1‑150 (MQ=255)
                                                                            gAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCg                                                                          >  2:705159/1‑150 (MQ=255)
                                                                              aGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGtt                                                                        <  1:5434190/150‑1 (MQ=255)
                                                                              aGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGtt                                                                        <  2:596503/150‑1 (MQ=255)
                                                                              aGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGtt                                                                        <  1:47369/150‑1 (MQ=255)
                                                                               ggCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTg                                                                       <  1:3179318/150‑1 (MQ=255)
                                                                                  gTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCGTTAGCAGATCAATAGTTGtaa                                                                    >  1:667480/1‑150 (MQ=255)
                                                                                      ctgcCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGTGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATcgc                                                                >  1:4998527/1‑150 (MQ=255)
                                                                                      ctgcCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATcgc                                                                <  1:3218707/150‑1 (MQ=255)
                                                                                         ccTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGcc                                                             >  2:5371665/1‑150 (MQ=255)
                                                                                               tGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAg                                                                                     >  2:4175615/1‑120 (MQ=255)
                                                                                               tGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAg                                                                                     <  1:4175615/120‑1 (MQ=255)
                                                                                                gAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAgccgcc                                                      >  2:1452330/1‑150 (MQ=255)
                                                                                                   cGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCAt                                                   <  1:2694108/150‑1 (MQ=255)
                                                                                                   cGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCACTTAGCGGGCGCGCCCCCGACTGACTGTAGCAGATAAATCGATTAAATAACGCCCGccccccct                                                   >  1:42942/1‑148 (MQ=255)
                                                                                                        aCGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGGGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCAtttctt                                              <  1:3041066/150‑1 (MQ=255)
                                                                                                            tGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCa                                          >  2:3658907/1‑150 (MQ=255)
                                                                                                            tGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCa                                          >  2:1418146/1‑150 (MQ=255)
                                                                                                            tGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCa                                          >  1:1773895/1‑150 (MQ=255)
                                                                                                            tGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTGTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATAAATAGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCa                                          >  2:1418389/1‑150 (MQ=255)
                                                                                                               cAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCg                                       <  2:297275/150‑1 (MQ=255)
                                                                                                               cAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCGTTAGCAGATCAATAGTTTAAATAGCGCCACCCGCCCCATTATTTACAGCg                                       >  1:935876/1‑150 (MQ=255)
                                                                                                                 aGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGt                                     <  1:718489/150‑1 (MQ=255)
                                                                                                                 aGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGt                                     <  2:1251690/150‑1 (MQ=255)
                                                                                                                  gCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGtt                                    <  2:1146921/150‑1 (MQ=255)
                                                                                                                    tGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGa                                  <  1:825786/150‑1 (MQ=255)
                                                                                                                    tGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGa                                  >  2:103193/1‑150 (MQ=255)
                                                                                                                    tGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGa                                  <  1:2728676/150‑1 (MQ=255)
                                                                                                                      tAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACg                                >  2:1592810/1‑150 (MQ=255)
                                                                                                                                gTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGTTCAATCGTTGAAATAGCGCCAGCCGCCCATTTAGTTGCAGCGGTTGACGGTGCCGCAGa                      >  1:2877418/1‑150 (MQ=255)
                                                                                                                                   cATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGCAAGCCTTTAGCAGATGAATCGGTGAAATAGCGCCAGCCCCCCCATTCTTTCCAGCGGTTGACGATGCCGCCGaaaa                   >  1:1021769/1‑148 (MQ=255)
                                                                                                                                   cATAAAAAAACCGGCCACGAAGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACa                   >  2:3522529/1‑150 (MQ=255)
                                                                                                                                    ataaaaaaaCCGGCCACGAAGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGa                                                                      <  1:1802415/95‑1 (MQ=37)
                                                                                                                                    ataaaaaaaCCGGCCACGAAGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGa                                                                      >  2:1802415/4‑98 (MQ=37)
                                                                                                                                    ataaaaaaaCCGGCCACGAAGGCCGGTTTTTTCATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGCATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAg                  >  1:1776133/4‑150 (MQ=38)
                                                                                                                                    aTGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAg                  <  2:22261/150‑1 (MQ=255)
                                                                                                                                     taaaaaaaCCGGCCACGAAGGCCGGTTTTTTCATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGc                 >  2:1792066/3‑150 (MQ=39)
                                                                                                                                      gAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGGCCGCGAATGCCTTTAGCAGATCAAACCTTGAACACGCGGCAGCCGCCCATTTGTTTCCAGCGGTTGACGAGTCCGCAGAACAGCt                >  1:3195318/1‑150 (MQ=255)
                                                                                                                                      gAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCGTTAGCAGATCAATCGTTGAAATAGCGCCAGCCGCCCCTTTCTTTCCAGCGGTTGACGATGCCGCCGAACAGcg                >  1:1643964/1‑149 (MQ=255)
                                                                                                                                       aaaaaaCCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGCTc               >  2:2886436/1‑150 (MQ=255)
                                                                                                                                       aaaaaaCCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGCTc               >  2:160153/1‑150 (MQ=255)
                                                                                                                                        aaaaaCCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGCTCa              >  2:3976382/1‑150 (MQ=255)
                                                                                                                                         aaaaCCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGCTCAg             >  2:2525050/1‑150 (MQ=255)
                                                                                                                                         aaaaCCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAAGGCCGTTAGCAGATCAAACGGTGGAATCGCGCCACCCGCCCAATTATTTCCACCGGGTGACGATGCCGCCGAACAGCACAg             >  1:1035567/1‑150 (MQ=255)
                                                                                                                                           aaCCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCGGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGCTCAGCg           >  2:4806314/1‑150 (MQ=255)
                                                                                                                                                  cTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTATTTCCAGCGGTTGACGATGCCGCCGAACAGCTCAGCGGttttat    >  1:5348751/1‑148 (MQ=255)
                                                                                                                                                   ttCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGCTCAGCGGTTTtctc   >  2:456126/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GCAGTTCCACGAAGAGCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGCTCAGCGGTTTTCTCG  >  NC_002947/1243597‑1243894

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: