Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,243,725 G→A 5.8% intergenic (+14/+85) tsaA → / ← rnt peroxiredoxin/ribonuclease T

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,243,7250GA 5.8% 265.2 / 3.2 86intergenic (+14/+85)tsaA/rntperoxiredoxin/ribonuclease T
Reads supporting (aligned to +/- strand):  ref base G (38/43);  new base A (3/2);  total (41/45)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.55e-01

TGACGCTCTGCAGTTCCACGAAGAGCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGC  >  NC_002947/1243582‑1243873
                                                                                                                                               |                                                                                                                                                        
tGACGCTCTGCAGTTCCACGAAGAGCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑aaaa                                                                                                                                                     >  4:182885/1‑148 (MQ=255)
   cGCTCTGCAGTTCCACGAAGAGCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAAcc                                                                                                                                                 >  2:422666/1‑149 (MQ=255)
    gCTCTGCAGTTCCACGAAGAGCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCg                                                                                                                                                <  5:91652/149‑1 (MQ=255)
        tGCAGTTCCACGAAGAGCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATA‑AAAAAACCGGcca                                                                                                                                            >  5:194906/1‑148 (MQ=255)
          cAGTTCCACGAAGAGCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTc                                                                                                                                          >  2:47618/1‑149 (MQ=255)
            gTTCCACGAAGAGCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATAAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGt                                                                                                                                        >  4:36473/1‑149 (MQ=255)
             ttCCACGAAGAGCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTg                                                                                                                                       <  8:199152/149‑1 (MQ=255)
             ttCCACGAAGAGCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTg                                                                                                                                       <  7:258232/149‑1 (MQ=255)
               ccACGAAGAGCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGc                                                                                                                                     >  4:29192/1‑149 (MQ=255)
                cACGAAGAGCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGcc                                                                                                                                    >  8:188570/1‑149 (MQ=255)
                cACGAAGAGCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGcc                                                                                                                                    <  7:393177/149‑1 (MQ=255)
                 aCGAAGAGCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAAACCGGCCTTCGTGGcc                                                                                                                                    <  4:386770/149‑1 (MQ=255)
                 aCGAAGAGCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCg                                                                                                                                   >  3:148525/1‑149 (MQ=255)
                 aCGAAGAGCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCg                                                                                                                                   <  2:214189/149‑1 (MQ=255)
                     agagCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑ttt                                                                                                                              >  3:298837/1‑149 (MQ=255)
                      gagCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCGGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCACGAAGGCCGG‑tttt                                                                                                                             <  6:303874/149‑1 (MQ=255)
                       agCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑ttttt                                                                                                                            >  5:231142/1‑149 (MQ=255)
                        gCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑tttttt                                                                                                                           >  2:247845/1‑149 (MQ=255)
                        gCACGGCGAAGGGGGCCCTGCCAACTGGAAAACAGGCGACAAGGGCATGAACGCGTCGCCTGAAGGCGGTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑tttttt                                                                                                                           <  3:34576/149‑1 (MQ=255)
                         cACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑tttttt                                                                                                                           >  1:195683/1‑148 (MQ=255)
                         cACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑tttttt                                                                                                                           >  6:315553/1‑148 (MQ=255)
                         cACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATA‑AAAAAACCGGCCACGAAGGCCGG‑ttttttc                                                                                                                          <  1:419461/149‑2 (MQ=39)
                          aCGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCg                                                                                                                                   <  3:49794/140‑1 (MQ=255)
                          aCGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCg                                                                                                                                   >  4:49794/1‑140 (MQ=255)
                           cGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTAt                                                                                                                        <  1:2504/149‑1 (MQ=255)
                           cGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTAt                                                                                                                        <  6:437085/149‑1 (MQ=255)
                             gCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAAAGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATg                                                                                                                       >  7:206200/1‑148 (MQ=255)
                                 aGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTc                                                                                                                  <  4:391456/149‑1 (MQ=255)
                                 aGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTc                                                                                                                  <  1:99372/149‑1 (MQ=255)
                                   tgtgCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCt                                                                                                                <  2:327585/149‑1 (MQ=255)
                                        ccTGCCAACTGGAAAAAAGGCGACAAGGCCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGAAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGcc                                                                                                           <  5:339678/149‑1 (MQ=255)
                                         cTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑TAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCg                                                                                                          <  7:257873/149‑1 (MQ=255)
                                            ccAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGcc                                                                                                       >  7:335645/1‑149 (MQ=255)
                                              aaCTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTc                                                                                                                  >  1:368022/1‑136 (MQ=255)
                                              aaCTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTc                                                                                                                  <  2:368022/136‑1 (MQ=255)
                                              aaCTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTACATG‑AAAAAACCGGCCTTCCTGGACGG‑TTTTATTATGCCTTCCTGTGCCGTccac                                                                                                     >  3:272033/1‑147 (MQ=255)
                                                    aaaaaaGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTcgcg                                                                                               <  4:148525/149‑1 (MQ=255)
                                                     aaaaaGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCgg                                                                                              <  8:58202/149‑1 (MQ=255)
                                                     aaaaaGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCgg                                                                                              <  2:396788/149‑1 (MQ=255)
                                                             gACAAGGGCATAAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑Gccc                                                                                    <  3:36473/149‑1 (MQ=255)
                                                              aCAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCcg                                                                                   >  6:49445/1‑149 (MQ=255)
                                                               cAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCcgc                                                                                  <  1:310368/149‑1 (MQ=255)
                                                                   ggCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAAt                                                                              <  1:422666/149‑1 (MQ=255)
                                                                      aTGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑ttttttt                                                                                                                          <  3:413528/104‑1 (MQ=255)
                                                                      aTGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑ttttttt                                                                                                                          >  4:413528/1‑104 (MQ=255)
                                                                      aTGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGcc                                                                           >  7:305269/1‑149 (MQ=255)
                                                                        gAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCtt                                                                         >  3:129118/1‑149 (MQ=255)
                                                                         aaCGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCtt                                                                         <  6:231142/148‑1 (MQ=255)
                                                                          aCGCTTCGCCTGAAGGCGTTTCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCTTTa                                                                       <  4:73810/149‑1 (MQ=255)
                                                                          aCGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCTTTa                                                                       <  6:129443/149‑1 (MQ=255)
                                                                           cGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCTTTAg                                                                      <  3:182885/149‑1 (MQ=255)
                                                                                 gCCTGAAGGCGTTGCTGCCTACCTGCGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCTTTAGCAGATc                                                                <  1:203086/149‑1 (MQ=255)
                                                                                  ccTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCTTTAGCAGATc                                                                >  3:181567/1‑148 (MQ=255)
                                                                                    tGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTATTTATGCCAACCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCTTTAGCAGATCAAt                                                             >  1:356262/1‑149 (MQ=255)
                                                                                        ggCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCTTTAGCAGATCAATCGtt                                                         >  4:167930/1‑149 (MQ=255)
                                                                                         gCGTTGCTGCCTACCTGAGCTAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCTTTAGCAGATCAATCGTTg                                                        <  1:247845/149‑1 (MQ=255)
                                                                                          cGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCTTTAGCAGATCAATCGTTg                                                        >  1:398707/1‑148 (MQ=255)
                                                                                           gTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAACCCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCTTTAGCAGATCAATCGTTGaa                                                      <  8:289195/149‑1 (MQ=255)
                                                                                           gTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGGTTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCTTTAGCAGATCAATCGTTg                                                        >  4:297052/1‑148 (MQ=255)
                                                                                             tgctgcCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCTTTAGCAGATCAATCGTTGa                                                       <  6:80531/146‑1 (MQ=255)
                                                                                             tgctgcCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCTTTAGCAGATCAATCGTTGa                                                       >  5:80531/1‑146 (MQ=255)
                                                                                                 gcCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATcgcg                                                >  1:106326/1‑149 (MQ=255)
                                                                                                 gcCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATcgcg                                                >  5:136908/1‑149 (MQ=255)
                                                                                                  ccTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATcgc                                                 <  2:195683/147‑1 (MQ=255)
                                                                                                     aCCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAg                                            >  2:136947/1‑149 (MQ=255)
                                                                                                        tGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAgccg                                         >  6:422786/1‑149 (MQ=255)
                                                                                                        tGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAgccg                                         >  2:39711/1‑149 (MQ=255)
                                                                                                                 aCGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCAtttct                                >  6:332688/1‑149 (MQ=255)
                                                                                                                     tGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCAtttctttc                             >  5:343273/1‑148 (MQ=255)
                                                                                                                       gCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCa                           <  5:422786/148‑1 (MQ=255)
                                                                                                                         aaGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCg                        <  4:181567/149‑1 (MQ=255)
                                                                                                                           gCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATcgc                                                 >  8:404987/1‑122 (MQ=255)
                                                                                                                           gCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATcgc                                                 <  7:404987/122‑1 (MQ=255)
                                                                                                                            cTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGtt                     <  5:332688/149‑1 (MQ=255)
                                                                                                                               tAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATg                                                                             <  3:7682/90‑1 (MQ=255)
                                                                                                                               tAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATg                                                                             >  4:7682/1‑90 (MQ=255)
                                                                                                                               tAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGa                   <  5:49445/148‑1 (MQ=255)
                                                                                                                                     ctgcGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGcc            <  4:298837/149‑1 (MQ=255)
                                                                                                                                     ctgcGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGcc            <  1:47618/149‑1 (MQ=255)
                                                                                                                                        cGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGcag           >  2:114695/1‑145 (MQ=255)
                                                                                                                                           gCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGaa      >  8:385392/1‑149 (MQ=255)
                                                                                                                                            cATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGaa      <  4:129118/148‑1 (MQ=255)
                                                                                                                                            cATA‑AAAAAACCGGCCACGAAGGCCGG‑TTTTTTCATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCa                           >  1:270062/1‑127 (MQ=38)
                                                                                                                                            cATA‑AAAAAACCGGCCACGAAGGCCGG‑TTTTTTCATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCa                           <  2:270062/127‑1 (MQ=38)
                                                                                                                                            cATA‑AAAAAACCGGCCACGAAGGCCGG‑TTTTTTCATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGaa      >  8:324327/1‑148 (MQ=38)
                                                                                                                                             aTG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACa    <  7:385392/149‑1 (MQ=255)
                                                                                                                                               g‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCAGGCCCGGGAATGCCGTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGaa      >  1:131320/1‑147 (MQ=255)
                                                                                                                                               g‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCTGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGc  >  3:310159/1‑149 (MQ=255)
                                                                                                                                               g‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAgg  >  5:133940/1‑148 (MQ=255)
                                                                                                                                               |                                                                                                                                                        
TGACGCTCTGCAGTTCCACGAAGAGCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGC  >  NC_002947/1243582‑1243873

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