Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,187,984 C→A 5.6% R163R (CGC→CGA)  PP_4567 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,187,9840CA 5.6% 276.1 / 2.5 85R163R (CGC→CGA) PP_4567hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (42/38);  new base A (3/2);  total (45/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGATTTCGGCTGCGCCTGTGAAACCGAACAAGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGTTACCGAAGATCAGCAAGTCC  >  NC_002947/5187836‑5188132
                                                                                                                                                    |                                                                                                                                                    
tGGATTTCGGCTGCGCCTGTGAAACCGAACAAGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTgcgc                                                                                                                                                      >  1:134130/1‑149 (MQ=255)
  gATTTCGGCTGCGCCTGTGAAACCGAACAAGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGACGGCGTGTTGcgcca                                                                                                                                                    >  8:374047/1‑149 (MQ=255)
   aTTTCGGCTGCGCCTGTGAAACCGAACAAGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGcgccac                                                                                                                                                   <  6:154224/149‑1 (MQ=255)
    tttCGGCTGCGCCTGTGAAACCGAACAAGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGcgccacg                                                                                                                                                  >  3:114720/1‑149 (MQ=255)
          cTGCGCCTGTGAAACCGAACAAGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATcc                                                                                                                                            <  2:134130/149‑1 (MQ=255)
             cgcCTGTGAAACCGAACAAGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATccc                                                                                                                                           <  5:101140/147‑1 (MQ=255)
               ccTGTGAAACCGAACAAGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGaa                                                                                                                                       <  1:269095/149‑1 (MQ=255)
               ccTGTGAAACCGAACAAGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGaa                                                                                                                                       <  8:241449/149‑1 (MQ=255)
               ccTGTGAAACCGAACAAGCAGCGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCg                                                                                                                                         >  6:354962/1‑147 (MQ=255)
                   tgAAACCGAACAAGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGGGGATGGCGTGGCGCAACACGCCTTcct                                                                                                                                           >  2:124449/1‑140 (MQ=255)
                   tgAAACCGAACAAGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGGGGATGGCGTGGCGCAACACGCCTTcct                                                                                                                                           <  1:124449/141‑2 (MQ=255)
                   tgAAACCGAACAAGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAaga                                                                                                                                     <  4:354113/147‑1 (MQ=255)
                   tgAAACCGAACAAGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGt                                                                                                                                   >  4:376273/1‑149 (MQ=255)
                   tgAAACCGAACAAGCAGCGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAaga                                                                                                                                     <  5:354962/147‑1 (MQ=255)
                     aaaCCGAACAAGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTc                                                                                                                                 >  6:339690/1‑149 (MQ=255)
                     aaaCCGAACAAGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTc                                                                                                                                 >  5:310870/1‑149 (MQ=255)
                     aaaCCGAACAAGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGAGTGTTGCGCCACGCCATCCCCGAAGAGTTc                                                                                                                                 >  8:221504/1‑149 (MQ=255)
                       aCCGAACAAGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCg                                                                                                                               >  6:175132/1‑149 (MQ=255)
                          gAACAAGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTa                                                                                                                            <  3:376273/149‑1 (MQ=255)
                           aaCAAGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAg                                                                                                                           >  8:41571/1‑149 (MQ=255)
                           aaCAAGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAg                                                                                                                           <  2:110183/149‑1 (MQ=255)
                            aCAAGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAgg                                                                                                                          <  3:357426/149‑1 (MQ=255)
                            aCAAGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAgg                                                                                                                          <  8:186457/149‑1 (MQ=255)
                               aGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGcc                                                                                                                       >  7:359275/1‑149 (MQ=255)
                                gCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGATGAAGGCGTGTTGCGCCACTCCATCCCCGAAGAGTTCCGCTAGGGCCt                                                                                                                      <  7:221504/149‑1 (MQ=255)
                                gCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAgg                                                                                                                          <  4:114768/145‑1 (MQ=255)
                                gCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAgg                                                                                                                          >  3:114768/1‑145 (MQ=255)
                                    tGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGc                                                                                                                  >  4:363389/1‑149 (MQ=255)
                                      aaGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGaa                                                                                                                                       <  5:402070/126‑1 (MQ=255)
                                      aaGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGaa                                                                                                                                       >  6:402070/1‑126 (MQ=255)
                                           gTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGc                                                                                                            >  2:173010/1‑148 (MQ=255)
                                             tGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAg                                                                                                         >  2:221234/1‑149 (MQ=255)
                                              ggCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCt                                                                                                                             >  7:367220/1‑128 (MQ=255)
                                              ggCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCt                                                                                                                             <  8:367220/128‑1 (MQ=255)
                                               gCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGgc                                                                                                       <  8:359275/149‑1 (MQ=255)
                                                  gTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGgcgcg                                                                                                    >  1:112386/1‑149 (MQ=255)
                                                   tCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGgcgcgc                                                                                                   <  5:175132/149‑1 (MQ=255)
                                                       ggTCCCGCCGGATCGCTCGCTCGATACGCTGCTCTACTACCTGCTCAGCCAGCGCTAAGCCGGCGTGCTGGAGTTCGTTGAAGGCGTGTTGCGACACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGGTTCGGCCAGgcgctatgcg                                                                                               <  6:413853/149‑7 (MQ=255)
                                                        gTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGcggcgg                                                                                              <  5:75491/149‑1 (MQ=255)
                                                         tCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCggg                                                                                             <  6:99001/149‑1 (MQ=255)
                                                           ccgccgCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCt                                                                                           >  7:122686/1‑149 (MQ=255)
                                                             gccgcAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTg                                                                                         <  3:363389/149‑1 (MQ=255)
                                                              ccgcAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGc                                                                                        >  1:146002/1‑149 (MQ=255)
                                                                  aaCGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGaaa                                                                                    <  1:173010/149‑1 (MQ=255)
                                                                   aCGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAac                                                                                   <  4:114720/149‑1 (MQ=255)
                                                                         gctcgATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGaa                                                                              >  7:326853/1‑148 (MQ=255)
                                                                         gctcgATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGaa                                                                              >  7:326854/1‑148 (MQ=255)
                                                                         gctcgATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGaa                                                                              <  7:41571/148‑1 (MQ=255)
                                                                         gctcgATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGaa                                                                              <  2:242875/148‑1 (MQ=255)
                                                                               tACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAAcagc                                                                          <  3:192633/146‑1 (MQ=255)
                                                                               tACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAAcagc                                                                          >  4:192633/1‑146 (MQ=255)
                                                                               tACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCAAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACagcagc                                                                       <  2:64635/149‑1 (MQ=255)
                                                                                aCCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGcc                                                                      >  4:363966/1‑149 (MQ=255)
                                                                                aCCCTGCTCTACTACCTGCTCACCCAGCGCGAAGACGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCCCGGCGGGCTTGCGAAACACGAACAGCAGcc                                                                      >  5:282300/1‑149 (MQ=255)
                                                                                  ccTGCTCTACTACCTGCTCACCCAGCGCCAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGccc                                                                     >  1:226393/1‑148 (MQ=255)
                                                                                      ctctACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTAc                                                                  <  7:125254/147‑1 (MQ=255)
                                                                                      ctctACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCa                                                                <  1:221234/149‑1 (MQ=255)
                                                                                         tactacCTGCTCACCCAGCGCGGAGCCGGACTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCAACCACg                                                              >  7:256641/1‑148 (MQ=255)
                                                                                         tactacCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCTCGa                                                             >  3:255483/1‑149 (MQ=255)
                                                                                          actacCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCATGCGAAACACGAACAGCAGCCCTACCACGa                                                             >  8:183882/1‑148 (MQ=255)
                                                                                            tacCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGAt                                                            >  5:72018/1‑147 (MQ=255)
                                                                                              ccTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGc                                                        >  1:16269/1‑149 (MQ=255)
                                                                                               cTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGcc                                                       >  6:82802/1‑149 (MQ=255)
                                                                                                 gCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGccc                                                      >  1:280453/1‑148 (MQ=255)
                                                                                                     aCCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATg                                                  <  2:112386/148‑1 (MQ=255)
                                                                                                     aCCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGc                                                 >  5:63819/1‑149 (MQ=255)
                                                                                                      cccAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCACGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGcc                                                >  6:340820/1‑149 (MQ=255)
                                                                                                         aGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCcag                                             >  8:161854/1‑149 (MQ=255)
                                                                                                           cgcgAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGaaa                                                                                    >  3:117046/1‑108 (MQ=255)
                                                                                                           cgcgAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGaaa                                                                                    <  4:117046/108‑1 (MQ=255)
                                                                                                               aaGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCt                                       <  6:310870/149‑1 (MQ=255)
                                                                                                                 gCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCt                                       >  4:374154/1‑147 (MQ=255)
                                                                                                                      ccTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGt                                <  2:146002/149‑1 (MQ=255)
                                                                                                                                ttCGGGGATGGCGTGGCGCAACACGCCTTCCTCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGt                      >  7:4331/1‑149 (MQ=38)
                                                                                                                                ttCGGGGATGGCGTGGCGCAACACGCCTTCCTCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGt                      >  5:336140/1‑149 (MQ=38)
                                                                                                                                ttCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGt                      >  1:238988/1‑149 (MQ=255)
                                                                                                                                ttCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGt                      >  1:276576/1‑149 (MQ=255)
                                                                                                                                ttCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCACCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGt                      >  8:361150/1‑149 (MQ=255)
                                                                                                                                   gAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGc                         >  4:283511/1‑143 (MQ=255)
                                                                                                                                   gAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGc                         <  3:283511/143‑1 (MQ=255)
                                                                                                                                        aGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGTTACCGAAg              <  3:363966/149‑1 (MQ=255)
                                                                                                                                           cGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGTTACCGAAGATc           <  7:183882/149‑1 (MQ=255)
                                                                                                                                            gtgtTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGTTACCGAAGAt            >  1:211530/1‑147 (MQ=255)
                                                                                                                                            gtgtTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGTTACCGAAGATCa          <  6:72018/149‑1 (MQ=255)
                                                                                                                                            gtgtTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGACTGAGCGGTGCCAGGCGGTTACCGAAGATCa          <  8:122686/149‑1 (MQ=255)
                                                                                                                                             tgtTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGTTACCGAAGATCAg         <  2:276576/149‑1 (MQ=255)
                                                                                                                                              gtTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGTTACCGAAGATCAg         >  3:270429/1‑148 (MQ=255)
                                                                                                                                               ttGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGTTACCGAAg              <  6:19999/142‑1 (MQ=255)
                                                                                                                                               ttGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGTTACCGAAg              >  5:19999/1‑142 (MQ=255)
                                                                                                                                                 gcgcCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGTTACCGAAGATCAGc        >  6:310647/1‑146 (MQ=255)
                                                                                                                                                 gcgcCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGTTACCGAAGATCAGc        <  5:310647/146‑1 (MQ=255)
                                                                                                                                                  cgccacgccaTCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGTTACCGAAGATCAGCAAg     <  4:270429/148‑1 (MQ=255)
                                                                                                                                                  cgccacgccaTCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGTTACCGAAGATCAGCAAg     <  2:211530/148‑1 (MQ=255)
                                                                                                                                                   gccacgccaTCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGTTACCGAAGATCAGCAAGTc   <  7:306606/149‑1 (MQ=255)
                                                                                                                                                    ccacgccaTCCGCGAAGAGTTCCTCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGTTACCGAAGATCAGCAAGTcc  <  7:161854/149‑1 (MQ=255)
                                                                                                                                                    ccacgccaTCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGTTACCGAAGATCAGCAAGTcc  <  3:374154/149‑1 (MQ=255)
                                                                                                                                                    ccacgccaTCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGTTACCGAAGATCAGCAAGTcc  <  5:339690/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TGGATTTCGGCTGCGCCTGTGAAACCGAACAAGCAGTGAAGATGTTGGCGGTCGAGGTCCCGCCGCAACGCTCGCTCGATACCCTGCTCTACTACCTGCTCACCCAGCGCGAAGCCGGCCTGCTGGACTTCGAGGAAGGCGTGTTGCGCCACGCCATCCCCGAAGAGTTCCGCTAGGGCCTACGCTTCGGCCAGGCGCGCGGCGGGCTTGCGAAACACGAACAGCAGCCCTACCACGATCAGCCCCATGCCCAGCAGGCTGAGCGGTGCCAGGCGGTTACCGAAGATCAGCAAGTCC  >  NC_002947/5187836‑5188132

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