Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,243,740 C→G 10.1% 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,7400CG10.1% 231.1 / 13.8 79intergenic (+29/+70)tsaA/rntperoxiredoxin/ribonuclease T
Reads supporting (aligned to +/- strand):  ref base C (37/34);  new base G (3/5);  total (40/39)
Fisher's exact test for biased strand distribution p-value = 4.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGTTCCACGAAGAGCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTG‑AAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGCTCAGCGGTTTTCTC  >  NC_002947/1243592‑1243887
                                                                                                                                                    |                                                                                                                                                      
cAGTTCCACGAAGAGCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTc                                                                                                                                                        <  8:10652/149‑1 (MQ=255)
     ccACGAAGAGCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTgg                                                                                                                                                    <  7:112569/148‑1 (MQ=255)
       aCGAAGAGCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCg                                                                                                                                                 >  5:130047/1‑149 (MQ=255)
           agagCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGttt                                                                                                                                             <  3:278368/149‑1 (MQ=255)
            gagCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATAAAAAAACCGGCCACGAAGGCCGGtttt                                                                                                                                            <  8:166575/149‑1 (MQ=255)
             agCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGttttt                                                                                                                                           >  7:116026/1‑149 (MQ=255)
                 cGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGtttttagat                                                                                                                                       >  5:162510/1‑145 (MQ=255)
                 cGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGtttt                                                                                                                                            <  6:162510/144‑1 (MQ=255)
                 cGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTAt                                                                                                                                       <  8:23049/149‑1 (MQ=255)
                  ggCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATg                                                                                                                                      <  6:130047/149‑1 (MQ=255)
                  ggCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATAAAAAAACCGGCCACGAAGGCCGGTTTTTTCATg                                                                                                                                      <  5:262059/149‑1 (MQ=38)
                  ggCGAAGTGTCCCCTGCCAACTCGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGAGCCCGGTTTTTTTATg                                                                                                                                      >  8:287032/1‑149 (MQ=255)
                          gtgCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCtg                                                                                                                              >  4:209825/1‑149 (MQ=255)
                              ccTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGcc                                                                                                                          >  8:41427/1‑149 (MQ=255)
                                  ccAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGc                                                                                                                       >  3:243066/1‑148 (MQ=255)
                                       tGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTc                                                                                                               >  4:24197/1‑149 (MQ=255)
                                               aGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGgcgc                                                                                                       >  8:95374/1‑149 (MQ=255)
                                                   gACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGccc                                                                                                   >  7:162949/1‑149 (MQ=255)
                                                   gACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGccc                                                                                                   >  1:297134/1‑149 (MQ=255)
                                                     cAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCcgc                                                                                                 <  8:29526/149‑1 (MQ=255)
                                                     cAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCcgc                                                                                                 <  8:162949/149‑1 (MQ=255)
                                                        gggCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATAAAAAAACCGGCCACGAAGGCCGGTTTTTTCATGc                                                                                                                                     <  4:69525/112‑1 (MQ=37)
                                                        gggCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATAAAAAAACCGGCCACGAAGGCCGGTTTTTTCATGc                                                                                                                                     <  4:65354/112‑1 (MQ=37)
                                                        gggCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATAAAAAAACCGGCCACGAAGGCCGGTTTTTTCATGc                                                                                                                                     >  3:65354/1‑112 (MQ=37)
                                                        gggCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATAAAAAAACCGGCCACGAAGGCCGGTTTTTTCATGc                                                                                                                                     >  3:69525/1‑112 (MQ=37)
                                                        gggCAGGAAGGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAAGGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGaa                                                                                              <  5:118849/149‑1 (MQ=255)
                                                         ggCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAAt                                                                                             >  6:180530/1‑149 (MQ=255)
                                                         ggCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAAt                                                                                             <  5:176677/149‑1 (MQ=255)
                                                           cATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCcgcg                                                                                                <  6:192466/144‑1 (MQ=255)
                                                           cATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGaagat                                                                                           >  5:192466/1‑146 (MQ=255)
                                                           cATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGc                                                                                           >  1:266611/1‑149 (MQ=255)
                                                           cATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGc                                                                                           <  3:147280/149‑1 (MQ=255)
                                                            aTGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGcc                                                                                          >  5:281556/1‑149 (MQ=255)
                                                            aTGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGcc                                                                                          >  2:193712/1‑149 (MQ=255)
                                                            aTGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGGTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGcg                                                                                          >  6:241336/1‑148 (MQ=255)
                                                              gAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCtt                                                                                        >  6:290077/1‑149 (MQ=255)
                                                               aaCGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCtt                                                                                        <  3:278678/148‑1 (MQ=255)
                                                                 cGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTAg                                                                                     <  7:41427/149‑1 (MQ=255)
                                                                    ttCGCCTGAAGGCGTTGCTGCCTACCTGAGCCAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTAg                                                                                     >  7:198636/1‑146 (MQ=255)
                                                                      cGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAg                                                                                  <  1:193712/147‑1 (MQ=255)
                                                                       gCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGcc                                                                                                                                                  >  5:46661/1‑84 (MQ=255)
                                                                       gCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGcc                                                                                                                                                  <  6:46661/84‑1 (MQ=255)
                                                                         cTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCaa                                                                             >  6:93412/1‑149 (MQ=255)
                                                                              ggCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGtt                                                                        <  4:37653/149‑1 (MQ=255)
                                                                               gCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGaagaggcc                                                                                        >  2:114623/1‑126 (MQ=255)
                                                                               gCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGaagaggcc                                                                                        <  1:114623/132‑7 (MQ=255)
                                                                                  ttGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTG‑aaa                                                                   <  8:116026/149‑1 (MQ=255)
                                                                                      tgcCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTG‑a                                                                     >  7:290336/1‑143 (MQ=255)
                                                                                      tgcCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTG‑a                                                                     <  8:290336/143‑1 (MQ=255)
                                                                                       gcCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTG‑AAATcgcg                                                              >  4:169671/1‑149 (MQ=255)
                                                                                        ccTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTG‑AAATcgcg                                                              <  3:209825/148‑1 (MQ=255)
                                                                                        ccTACCTGAGCGAGAACGAGGGCAAGCTGTAATTGGTGCGTGCATGAAAAAACCGGCCTTCGTGGCCTGTTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTAAGCAGATCAATCGTAG‑AAATcgcg                                                              <  4:215720/148‑1 (MQ=255)
                                                                                             cTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTG‑AAATCGCGCCAgcc                                                        >  8:146590/1‑149 (MQ=255)
                                                                                                  cGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTG‑AAATCGCGCCAGCCGCCCa                                                   >  2:272031/1‑149 (MQ=255)
                                                                                                   gagaACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTG‑AAATCGCGCCAGCCGCCCAt                                                  >  1:268121/1‑149 (MQ=255)
                                                                                                   gagaACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTG‑AAATCGCGCCAACCGCCCAt                                                  <  7:156972/149‑1 (MQ=255)
                                                                                                        cGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGgcgcgc                                                                                                     >  4:14930/1‑94 (MQ=255)
                                                                                                        cGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGgcgcgc                                                                                                     <  3:14930/94‑1 (MQ=255)
                                                                                                           tGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGAGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCCCGAATGCCTTTAGCACATCAATCGTTG‑AAAACGAGCCAGCCGCCCAatactttc                                           >  3:248145/1‑148 (MQ=255)
                                                                                                                  cTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTG‑AAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGtt                                   <  5:93412/149‑1 (MQ=255)
                                                                                                                  cTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTG‑AAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGtt                                   <  7:146590/149‑1 (MQ=255)
                                                                                                                   tGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTG‑AAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGtt                                   >  2:108204/1‑148 (MQ=255)
                                                                                                                    gTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTAg                                                                                     <  2:491/98‑1 (MQ=255)
                                                                                                                    gTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTAg                                                                                     >  1:491/1‑98 (MQ=255)
                                                                                                                     tAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGa                                 <  1:53988/149‑1 (MQ=255)
                                                                                                                     tAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTG‑AAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGtt                                   <  1:272031/146‑1 (MQ=255)
                                                                                                                             gcGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGAt                                                                                >  3:57025/1‑94 (MQ=255)
                                                                                                                             gcGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGAt                                                                                <  4:57025/94‑1 (MQ=255)
                                                                                                                               gTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCTGGCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTG‑AAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCg                         >  1:84051/1‑148 (MQ=255)
                                                                                                                                      aaaaaaCCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTG‑AAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGCt               >  1:85213/1‑149 (MQ=255)
                                                                                                                                      aaaaaaCCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCCTTG‑AAATCGAGCCAGCCGCCCATTTCTTTACAGAGGTTGAAGATGCAGCAGAACAGc                >  4:265999/1‑148 (MQ=255)
                                                                                                                                       aaaaaCCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTAGc                                                                                    >  8:14475/1‑80 (MQ=255)
                                                                                                                                       aaaaaCCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTAGc                                                                                    <  7:14475/80‑1 (MQ=255)
                                                                                                                                       aaaaaCCGGCCACGAAGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTG‑AAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGc                           >  4:49053/1‑136 (MQ=255)
                                                                                                                                       aaaaaCCGGCCACGAAGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTG‑AAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGc                           <  3:49053/136‑1 (MQ=255)
                                                                                                                                        aaaaCCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTG‑AAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGCTCa             <  7:95374/149‑1 (MQ=255)
                                                                                                                                           aCCGTCCTTCGTGTCCGGTTTTTTTTTGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGAGAATGCCTTTAGGAGATCAATCGTTG‑AAATCGCGCCAGCCGTCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGCTCAggg          <  8:198636/149‑3 (MQ=255)
                                                                                                                                              ggCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTG‑AAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGCTCAGCGGtt       <  3:169671/149‑1 (MQ=255)
                                                                                                                                               gCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTG‑AAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGCTCAGCGGttt      >  4:109479/1‑149 (MQ=255)
                                                                                                                                                ccTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTG‑AAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGCTCAGCGGtttt     <  5:290077/149‑1 (MQ=255)
                                                                                                                                                    cGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTG‑AAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGCTCAGCGGTTTtctc  <  6:281556/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
CAGTTCCACGAAGAGCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTG‑AAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGCTCAGCGGTTTTCTC  >  NC_002947/1243592‑1243887

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑