Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,243,741 G→A 10.0% intergenic (+30/+69) tsaA → / ← rnt peroxiredoxin/ribonuclease T

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,243,7410GA10.0% 220.6 / 10.7 80intergenic (+30/+69)tsaA/rntperoxiredoxin/ribonuclease T
Reads supporting (aligned to +/- strand):  ref base G (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 = 7.15e-01

CCACGAAGAGCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTG‑AAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGCTCAGCGGTTTTCTCG  >  NC_002947/1243597‑1243888
                                                                                                                                                |                                                                                                                                                      
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:162949/149‑1 (MQ=255)
                                                cAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCcgc                                                                                                  <  8:29526/149‑1 (MQ=255)
                                                   gggCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATAAAAAAACCGGCCACGAAGGCCGGTTTTTTCATGc                                                                                                                                      <  4:69525/112‑1 (MQ=37)
                                                   gggCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATAAAAAAACCGGCCACGAAGGCCGGTTTTTTCATGc                                                                                                                                      <  4:65354/112‑1 (MQ=37)
                                                   gggCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATAAAAAAACCGGCCACGAAGGCCGGTTTTTTCATGc                                                                                                                                      >  3:69525/1‑112 (MQ=37)
                                                   gggCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATAAAAAAACCGGCCACGAAGGCCGGTTTTTTCATGc                                                                                                                                      >  3:65354/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                                                                                           >  2:193712/1‑149 (MQ=255)
                                                       aTGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGcc                                                                                           >  5:281556/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                                                                                                                                                   <  6:46661/84‑1 (MQ=255)
                                                                  gCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGcc                                                                                                                                                   >  5:46661/1‑84 (MQ=255)
                                                                    cTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCaa                                                                              >  6:93412/1‑149 (MQ=255)
                                                                         ggCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGtt                                                                         <  4:37653/149‑1 (MQ=255)
                                                                          gCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGaagaggcc                                                                                         <  1:114623/132‑7 (MQ=255)
                                                                          gCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGaagaggcc                                                                                         >  2:114623/1‑126 (MQ=255)
                                                                             ttGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTG‑aaa                                                                    <  8:116026/149‑1 (MQ=255)
                                                                                 tgcCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTG‑a                                                                      <  8:290336/143‑1 (MQ=255)
                                                                                 tgcCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTG‑a                                                                      >  7:290336/1‑143 (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)
                                                                                                                                               cGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTG‑AAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGCTCAGCGGTTTtctc   <  6:281556/148‑1 (MQ=255)
                                                                                                                                                gTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTG‑AAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGCTCAGCGGTTTTCTCg  <  2:268121/148‑1 (MQ=255)
                                                                                                                                                gTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CATCGCGGGCGCTCCCGCGAATGCCTTTAGCAGATCAATGGTTG‑AAATCGCGGCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGCTCAGCGGTTTTCTCg  <  4:243066/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                      
CCACGAAGAGCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTG‑AAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGCTCAGCGGTTTTCTCG  >  NC_002947/1243597‑1243888

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