Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 406,031 C→G 7.6% S9T (AGC→ACC) ‡ PP_0337 ← sensory box protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947406,0310CG 7.6% 249.4 / 9.3 79S9T (AGC→ACC) ‡PP_0337sensory box protein
Reads supporting (aligned to +/- strand):  ref base C (33/40);  new base G (4/2);  total (37/42)
Fisher's exact test for biased strand distribution p-value = 4.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.65e-01

TCGCCGGTTTTCAGGCCGAACTGGCGCTCACAGGCTGGGTCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTGGCAAGCGGCCATGCTTGATGCTGTGGACAGGTTATCGG  >  NC_002947/405887‑406172
                                                                                                                                                  |                                                                                                                                               
tCGCCGGTTTTCAGGCCGAACTGGCACTCACAGGCTGGGTCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTggc                                                                                                                                             >  4:109884/1‑149 (MQ=255)
       ttttCAGGCCGAACTGGCGCTCACAGGCTGGGTCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCAGCCGTGTggcggccgag                                                                                                                                      <  1:185101/149‑5 (MQ=255)
       ttttCAGGCCGAACTGGCGCTCACAGGCTGGGTCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCg                                                                                                                                      >  4:283471/1‑149 (MQ=255)
            aGGCCGAACTGGCGCTCACCGGCTGGGTCAAGGTAGCGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCTGCATCCCGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGtt                                                                                                                                  <  1:276674/148‑1 (MQ=255)
            aGGCCGAACTGGCGCTCACAGGCTGGGTCAAGGTAGATGATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTGACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGtt                                                                                                                                  <  8:137104/148‑1 (MQ=255)
            aGGCCGAACTGGCGCTCACAGGCTGGGTCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGtt                                                                                                                                  <  6:7493/148‑1 (MQ=255)
            aGGCCGAACTGGCGCTCACAGGCTGGGTCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGtt                                                                                                                                  <  3:283471/148‑1 (MQ=255)
             ggCCGAACTGGCGCTCACAGGCTGGGTCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTgg                                                                                                                                >  2:83204/1‑149 (MQ=255)
                  aaCTGGCGCTCACAGGCTGGGTCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCtttt                                                                                                                           <  1:199637/149‑1 (MQ=255)
                         gcTCACAGGCTGGGTCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGaa                                                                                                                    <  1:125994/149‑1 (MQ=255)
                         gcTCACAGGCTGGGTCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGaa                                                                                                                    <  1:58133/149‑1 (MQ=255)
                         gcTCACAGGCTGGGTCAAGGTAGAGCATCGCCCAGGTGGCTTAGTTGAGACTTTCGAAACGCAGCATGCCGAGACGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACCCGGCTGGCGGCATCGGGTTCGCTTTTCATCGaa                                                                                                                    >  5:316284/1‑149 (MQ=255)
                           tCACAGGCTGGGTCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑Ga                                                                                                                 >  4:228202/1‑149 (MQ=255)
                             acaGGCTGGGTCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACt                                                                                                               >  3:149277/1‑149 (MQ=255)
                                      gTCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAag                                                                                                      <  8:98507/149‑1 (MQ=255)
                                       tCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGc                                                                                                            <  7:233007/142‑1 (MQ=255)
                                       tCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGc                                                                                                            >  8:233007/1‑142 (MQ=255)
                                         aaGGTAGAGCATGGCGCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAgg                                                                                                   <  3:5/149‑1 (MQ=255)
                                         aaGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAgg                                                                                                   <  4:234492/149‑1 (MQ=255)
                                             tAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑c                                                                                              <  1:83204/149‑1 (MQ=255)
                                               gagCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑cgc                                                                                            <  3:228202/149‑1 (MQ=255)
                                                 gCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCgg                                                                                          >  8:32213/1‑149 (MQ=255)
                                                 gCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCgg                                                                                          <  8:176709/149‑1 (MQ=255)
                                                  cATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGc                                                                                         <  1:233660/149‑1 (MQ=255)
                                                       cccAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGaa                                                                                                       <  5:333399/131‑1 (MQ=255)
                                                       cccAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGaa                                                                                                       >  6:333399/1‑131 (MQ=255)
                                                        ccAGGTGGGTTCGTTGAGCCTTTAGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCCCACGGCTGGCGGCATCGGGGTGGCTTTTAATCGAA‑GACTCGCTTGAAGAGGGCA‑CGCGGCCGtgtg                                                                                   >  2:110763/1‑149 (MQ=255)
                                                          aGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCgg                                                                                          >  6:36449/1‑140 (MQ=255)
                                                          aGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCgg                                                                                          <  5:36449/140‑1 (MQ=255)
                                                               gCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTtgt                                                                                    >  1:286565/1‑141 (MQ=255)
                                                               gCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTtgt                                                                                    <  2:286565/141‑1 (MQ=255)
                                                               gCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTTGTGGCTgcgg                                                                            <  7:32213/149‑1 (MQ=255)
                                                                    gTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTTGTGGCTGCGGCGGAt                                                                       >  3:205179/1‑149 (MQ=255)
                                                                      tGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTtgt                                                                                    >  2:235004/1‑134 (MQ=255)
                                                                      tGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTtgt                                                                                    <  1:235004/134‑1 (MQ=255)
                                                                        aGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTTGTGGCTGCGGCGGATTTg                                                                    <  1:42663/148‑1 (MQ=255)
                                                                         gCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCAGGGGTAACTGTGTCACG‑‑ACCTCGGCCGCGACACTGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTTGTGGCTGCGGCGGATTTGGc                                                                  <  2:93692/149‑1 (MQ=255)
                                                                          ccTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTTGTGGCTGCGGCGGATTTGGc                                                                  >  2:253161/1‑148 (MQ=255)
                                                                           cTTTCGAAACGCAGCATGCCGAGCCGCGCGGGCACGTGTAACTGTGTCACGTCACATCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTTGTGGCTGCGGCGGATTTGGc                                                                  <  1:110763/149‑1 (MQ=255)
                                                                            tttCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTTGTGGCTGCGGCGGATTTGGCAAg                                                               >  7:152422/1‑149 (MQ=255)
                                                                            tttCGAAACGCAGCATGCCGAGCCGCGAGGGCACGCGTAAGTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCGTTTAATGGAA‑GACGATCTTGAAAAGGGAA‑CGCTGCATTGTGCCTGAGGAGGATTTGCCAAg                                                               >  4:189327/1‑149 (MQ=255)
                                                                              tCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGt                                                             >  2:250470/1‑149 (MQ=255)
                                                                                gAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTg                                                           >  5:222882/1‑149 (MQ=255)
                                                                                    cGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCatg                                                       <  1:253161/149‑1 (MQ=255)
                                                                                          aTGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTcc                                                 <  3:109884/149‑1 (MQ=255)
                                                                                            gCCGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAACGACTCGCTTGAAGAGGGAA‑CGCGGCCTTGTGGCTGCGGCGGATTTg                                                                    <  1:245518/129‑1 (MQ=255)
                                                                                            gCCGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTTGTGGCTGCGGCGGATTTg                                                                    >  2:245518/1‑128 (MQ=255)
                                                                                              cGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAACCGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGa                                              <  1:250470/149‑1 (MQ=255)
                                                                                              cGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAg                                             >  8:315580/1‑149 (MQ=255)
                                                                                               gAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGa                                              <  4:285552/147‑1 (MQ=255)
                                                                                               gAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGa                                              >  3:285552/1‑147 (MQ=255)
                                                                                                   cgcgAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTg                                        >  5:253721/1‑149 (MQ=255)
                                                                                                          gCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTGGCAAGCg                                 <  3:166770/149‑1 (MQ=255)
                                                                                                          gCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTGGCAAGCg                                 <  6:253721/149‑1 (MQ=255)
                                                                                                          gCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTGGCAAGCg                                 <  4:205179/149‑1 (MQ=255)
                                                                                                           cACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTGGCAAGCgg                                >  1:161468/1‑149 (MQ=255)
                                                                                                            aCGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTGGCAAGCGGc                               <  8:285092/149‑1 (MQ=255)
                                                                                                                   aCTGTGTCACG‑‑ACCTCTGCCGCCACACGGCTGGCGGCATCGGGTTGTCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTTGTGGCTGCGGCGGTTTTGGCAAGGTTGCATTATGTCCCGAGGGCTGGCAAGCGGCCATGCt                         <  6:316284/148‑1 (MQ=255)
                                                                                                                   aCTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTGGCAAGCGGCCATGCtt                        >  7:269990/1‑149 (MQ=255)
                                                                                                                       tgtCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTGGCAAGCGGCCATGCTTGATg                    >  4:246680/1‑149 (MQ=255)
                                                                                                                         tCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTTGTGGCTGCGTCGGATTTGGCAAGGTTGCATCATGTAACGAGGGCTGGCAGGCGGCCCTGGTTgcttct                  >  8:223417/1‑144 (MQ=255)
                                                                                                                         tCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTGGCAAGCGGCCATGCTTGATGCt                  <  8:152422/149‑1 (MQ=255)
                                                                                                                         tCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTGGCAAGCGGCCATGCTTGATGCt                  >  2:147347/1‑149 (MQ=255)
                                                                                                                          ccgctctccTCGGCCGCCACACGGCTGGCGGGATCGGGTTGGCTGTTCATCGAA‑GACTCGCTTGCAGAGGGAA‑CGCGGCCTTGTGGCTGCGGCGGATTTGTCAAGGGTGCATGATGGCCCGAGGGCTGGCAAGCGGCCATGCTTGATg                    <  7:241907/141‑1 (MQ=255)
                                                                                                                           aCG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTGGCAAGCGGCCATGCtt                        <  1:209937/141‑1 (MQ=255)
                                                                                                                           aCG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTGGCAAGCGGCCATGCtt                        >  2:209937/1‑141 (MQ=255)
                                                                                                                           aCG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTGGCAAGCGGCCATGCTTGATGCtgt                >  7:72172/1‑149 (MQ=255)
                                                                                                                           aCG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTGGCAAGCGGCCATGCTTGATGCtgt                >  3:160547/1‑149 (MQ=255)
                                                                                                                            cG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTGGCAAGCGGCCATGCTTGATGCtgtg               <  7:315580/149‑1 (MQ=255)
                                                                                                                            cG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTGGCAAGCGGCCATGCTTGATGCtgtg               <  6:222882/149‑1 (MQ=255)
                                                                                                                            cG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTGGCAAGCGGCCATGCTTGATGCtgtg               <  3:246680/149‑1 (MQ=255)
                                                                                                                                 ccgatgccgccAGCCGTGTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTGGCAAGCGGCCATGCTTGATGCTGTGGa             >  6:58155/6‑148 (MQ=255)
                                                                                                                                 ccgatgccgccAGCCGTGTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTGGCAAGCGGCCATGCTTGATGCTGTGGa             <  5:58155/143‑1 (MQ=255)
                                                                                                                                  cTCGGCCGCCACACGGCTGGCGGCAGCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGATGGAA‑CGCGGCCTTGTGGCTGAGGCGGCTTTGGCAAGGTTGCATGATGGCCCGAGGGCTGGCACGAGGCCCTGCTAGCTGcagctcgg            >  8:241907/1‑141 (MQ=255)
                                                                                                                                   gatgccgccAGCCGTGTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GCCTCGCTTGAAGAGGGAA‑CGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTGGCAAGCGGCCATGCTTGATGCTGTGGACAg          >  3:291221/4‑149 (MQ=255)
                                                                                                                                   gatgccgccAGCCGTGTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTGGCAAGCGGCCATGCTTGATGCTGTGGACAg          >  3:108598/4‑149 (MQ=255)
                                                                                                                                    atgccgccAGCCGTGTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTGGCAAGCGGCCATGCTTGATGCTGTGGACAgg         >  3:192489/3‑149 (MQ=255)
                                                                                                                                          ccACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTGGCAAGCGGCCATGCTTGATGCTGTGGACAGGTTATCg   >  7:196789/1‑149 (MQ=255)
                                                                                                                                           cacaCGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTGGCAAGCGGCCATGCTTGATGCTGTGGACAGGTTATCgg  <  1:147347/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
TCGCCGGTTTTCAGGCCGAACTGGCGCTCACAGGCTGGGTCAAGGTAGAGCATGGCCCAGGTGGCTTCGTTGAGCCTTTCGAAACGCAGCATGCCGAGCCGCGAGGGCACGGGTAACTGTGTCACG‑‑ACCTCGGCCGCCACACGGCTGGCGGCATCGGGTTGGCTTTTCATCGAA‑GACTCGCTTGAAGAGGGAA‑CGCGGCCTTGTGGCTGCGGCGGATTTGGCAAGGTTGCATGATGTCCCGAGGGCTGGCAAGCGGCCATGCTTGATGCTGTGGACAGGTTATCGG  >  NC_002947/405887‑406172

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