Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,243,738 T→A 9.7% 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,243,7380TA 9.7% 224.2 / 9.3 80intergenic (+27/+72)tsaA/rntperoxiredoxin/ribonuclease T
Reads supporting (aligned to +/- strand):  ref base T (37/35);  new base A (3/5);  total (40/40)
Fisher's exact test for biased strand distribution p-value = 7.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGTTCCACGAAGAGCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTG‑AAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGCTCAGCGGTTTT  >  NC_002947/1243592‑1243884
                                                                                                                                                  |                                                                                                                                                     
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                                <  7:146590/149‑1 (MQ=255)
                                                                                                                  cTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTG‑AAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGtt                                <  5:93412/149‑1 (MQ=255)
                                                                                                                   tGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTG‑AAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGtt                                >  2:108204/1‑148 (MQ=255)
                                                                                                                    gTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTAg                                                                                  >  1:491/1‑98 (MQ=255)
                                                                                                                    gTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTAg                                                                                  <  2:491/98‑1 (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)
                                                                                                                                                  |                                                                                                                                                     
CAGTTCCACGAAGAGCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTG‑AAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGCTCAGCGGTTTT  >  NC_002947/1243592‑1243884

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: ‑