Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	326,065	C→A	9.0%	*556Y (TAG→TAT)	PP_0269 ←	glutamate synthase large subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	326,065	0	C	A	9.0%	193.3 / 9.4	67	*556Y (TAG→TAT)	PP_0269	glutamate synthase large subunit
Reads supporting (aligned to +/- strand):  ref base C (30/31);  new base A (3/3);  total (33/34)
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.60e-01

GTCGATCAAGCTGCGTAGCTCGTTCTTGCGTACCAACAACGCCGTGCAGCAGAACGGTTACAACGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCC‑GATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGTGCAAAGCTGTCGCTACGCGCCATCCGCCACATCCGCGCATAAAACTCGCTTTCGATCGAGCCACTGAGTAA  >  NC_002947/325917‑326205                                                                                                                                                      |                                                                                                                                             gTCGATCAAGCTGCGTAGCTCGTTCTTGCGTACCAACAACGCCGTGCAGCAGAACGGTTACAACGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCC‑GATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGcc                                                                                                                                              >  5:123982/1‑149 (MQ=255) gTCGATCAAGCTGCGTAGCTCGTTCTTGCGTACCAACAACGCCATGCAGCAGAACGGTTACAACGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCC‑GATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGcc                                                                                                                                              >  1:219700/1‑149 (MQ=255)    gATCAAGCTGCGTAGCTCGTTCTTGCGTACCAACAACGCCGTGCAGCAGAACGGTTACAACGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCC‑GATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTAt                                                                                                                                           >  4:253608/1‑149 (MQ=255)    gATCAAGCTGCGTAGCTCGTTCTTGCGTACCAACAACGCCGTGCAGCAGAACGGTTACAACGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCC‑GATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTAt                                                                                                                                           >  6:35337/1‑149 (MQ=255)     aTCAAGCTGCGTAGCTCGTTCTTGCGTACCAACAACGCCGTGCAGCAGAACGGTTACAACGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCC‑GATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCATagg                                                                                                                                          >  7:273350/1‑147 (MQ=255)         aGCTGCGTAGCTCGTTCTTGCGTACCAACAACGCCGTGCAGCAGAACGGTTACAACGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCC‑GATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCgg                                                                                                                                      <  6:78089/149‑1 (MQ=255)           cTGCGTAGCTCGTTCTTGCGTACCAACAACGCCGTGCAGCAGAACGGTTACAACGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCC‑GATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCgggg                                                                                                                                    <  4:167386/149‑1 (MQ=255)              cGTAGCTCGTTCTTGCGTACCAACAACGCCGTGCAGCAGAACGGTTACAACGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCCAGATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCATAGGCCGGGGCt                                                                                                                                  <  1:90677/149‑1 (MQ=255)               gTAGCTCGTTCTTGCGTACCAACAACGCCGTGCAGCAGAACGGTTACAACGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCC‑GATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGt                                                                                                                                >  4:135565/1‑149 (MQ=255)                tAGCTCGTTCTTGCGTACCAACAACGCCGTGCAGCAGAACGGTTACAACGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCC‑GATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGtt                                                                                                                               >  1:110273/1‑149 (MQ=255)                   cTCGTTCTTGCGTACCAACAACGCCGTGCAGCAGAACGGTTACAACGACGACGCCAACGAAGTCCGCGTATTCGTCGAGTACCC‑GATCAGCATCTACTGATACTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCAGTTTCt                                                                                                                            >  1:112110/1‑149 (MQ=255)                      gTTCTTGCGTACCAACAACGCCGTGCAGCAGAACGGTTACAACGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCC‑GATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCtttt                                                                                                                         >  1:191093/1‑149 (MQ=255)                      gTTCTTGCGTACCAACAACGCCGTGCAGCAGAACGGTTACAACGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCC‑GATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCttt                                                                                                                          <  3:253608/148‑1 (MQ=255)                      gTTCTTGCGTACCAACAACGCCGTGCAGCAGAACGGTTACAACGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCC‑GATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCATAGGCCGGGGCTGttcgata                                                                                                                          <  8:273350/148‑6 (MQ=255)                          ttGCGTACCAACAACGCCGTGCAGCAGAACGGTTACAACGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCC‑GATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGc                                                                                                                     >  7:57448/1‑149 (MQ=255)                            gCGTACCAACAACGCCGTGCAGCAGAACGGTTACAACGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCC‑GATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGc                                                                                                                     <  5:35337/147‑1 (MQ=255)                              gTACCAACAACGCCGTGCAGCAGAACGGTTACAACGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCC‑GATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAAc                                                                                                                 >  6:126300/1‑149 (MQ=255)                                 ccAACAACGCCGTGCAGCAGAACGGTTACAACGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCC‑TATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGt                                                                                                              <  7:223053/149‑1 (MQ=255)                                   aacaacGCCGTGCAGCAGAACGGTTACAACGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCC‑GATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCg                                                                                                                       >  6:281820/1‑138 (MQ=255)                                   aacaacGCCGTGCAGCAGAACGGTTACAACGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCC‑GATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCg                                                                                                                       <  5:281820/138‑1 (MQ=255)                                   aacaacGCCGTGCAGCAGAACGGTTACAACGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCC‑GATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTa                                                                                                            >  5:173009/1‑149 (MQ=255)                                    acaacGCCGTGCAGCAGAACGGTTACAACGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCC‑GATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTa                                                                                                            <  3:137072/148‑1 (MQ=255)                                    acaacGCCGTGCAGCAGAACGGTTACAACGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCC‑GATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTa                                                                                                            >  4:137072/1‑148 (MQ=255)                                     caacGCCGTGCAGCAGAACGGTTACAACGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCC‑GATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTc                                                                                                                        >  2:81004/1‑135 (MQ=255)                                     caacGCCGTGCAGCAGAACGGTTACAACGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCC‑GATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTc                                                                                                                        <  1:81004/135‑1 (MQ=255)                                     caacGCCGTGCAGCAGAACGGTTACAACGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCC‑GATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACt                                                                                                          <  6:123982/149‑1 (MQ=255)                                        cGCCGTGCAGCAGAACGGTTACAACGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCC‑GATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGtt                                                                                                       <  6:173009/149‑1 (MQ=255)                                         gCCGTGCAGCAGAACGGTTACAACGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCC‑GATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGttt                                                                                                      <  8:57448/149‑1 (MQ=255)                                         cccGTGCAGCTTAACGGTTAGATCGACGACGTCAACTAAGTCCGTGTATTCTTCGAGTACCC‑GTTCAGCATCTTCTTATTCTGATCTGATGTATCGAACATCCAGGGCCTATGCGGGGGGTGTTTCTTTTCGGCGATCGGTTTCTGttt                                                                                                      <  2:112110/148‑1 (MQ=255)                                                  cagAACGGTTACAACGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCC‑GATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCtt                                                                                                                           <  5:292143/119‑1 (MQ=255)                                                  cagAACGGTTACAACGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCC‑GATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCtt                                                                                                                           >  6:292143/1‑119 (MQ=255)                                                  cagAACGGTTACAACGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCC‑GATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCATAGGCCGGGGCt                                                                                                                                  >  6:145815/1‑112 (MQ=255)                                                  cagAACGGTTACAACGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCC‑GATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCATAGGCCGGGGCt                                                                                                                                  <  5:145815/112‑1 (MQ=255)                                                     aaCGGTTACAACGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCC‑GATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTg                                                                                          >  6:227810/1‑149 (MQ=255)                                                     aaCGGTTACAACGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCC‑GATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTg                                                                                          >  6:151401/1‑149 (MQ=255)                                                          ttACAACGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCC‑GATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGaa                                                                                                                  >  4:215124/1‑120 (MQ=255)                                                          ttACAACGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCC‑GATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGaa                                                                                                                  <  3:215124/120‑1 (MQ=255)                                                            acacCGACGACGGCAACGCAGTCCGCGTATTCGTCGAGTACCC‑GATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTa                                                                                    <  3:185596/148‑1 (MQ=255)                                                            aCAACGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCC‑GATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTAc                                                                                   >  4:17461/1‑149 (MQ=255)                                                             cAACGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCC‑GATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTAc                                                                                   <  2:191093/148‑1 (MQ=255)                                                             cAACGACGACGGCAAAGAAGTCCGCGTTTTCGTCGAGTACCC‑GATCCGCATCTTCTGATTCTTATCTGATTTATCGAACACCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACg                                                                                  <  8:165546/149‑1 (MQ=255)                                                                    gacgGCAACGAAGTCCGCGTATTCGTCGAGTACCC‑GATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCTAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCg                                                                           <  2:219700/149‑1 (MQ=255)                                                                    gacgGCAACGAAGTCCGCGTATTCGTCGAGTACCC‑GATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCg                                                                           <  5:151401/149‑1 (MQ=255)                                                                       ggCAACGAAGTCCGCGTATTCGTCGAGTACCC‑GATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGTg                                                                        <  2:110273/149‑1 (MQ=255)                                                                            cGAAGTCCGCGTATTCGTCGAGTACCC‑GATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGTGCAAAg                                                                   <  5:126300/149‑1 (MQ=255)                                                                                    gcgTATTCGTCGAGTACCC‑GATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGTGCAAAGCTGTCGCt                                                           <  3:135565/149‑1 (MQ=255)                                                                                            gtcgAGTACCC‑GATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGTGCAAAGCTGTCGCTACGCGcc                                                    >  8:130394/1‑148 (MQ=255)                                                                                                 gTACCC‑GATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGTGCAACGCTGTCGCTACGCGCCATCCGc                                              >  6:248297/1‑149 (MQ=255)                                                                                                 gTACCC‑GATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGTGCAAAGCTGTCGCTACGCGCCATCCGc                                              >  3:282452/1‑149 (MQ=255)                                                                                                  tACCC‑GATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGTGCAAAGCTGTCGCTACGCGCCATCCGcc                                             <  1:210137/149‑1 (MQ=255)                                                                                                     cc‑GATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGtctc                                                                             <  2:118403/114‑1 (MQ=255)                                                                                                        gATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGTGCAAAGCTGTCGCTACGCGCCATCCGCCACATc                                        >  2:181920/1‑149 (MQ=255)                                                                                                             gCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGTGCAAAGCTGTCGCTACGCGCCATCCGCCACATCcgcg                                    >  5:130424/1‑148 (MQ=255)                                                                                                                tcttctGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTg                                                                                          <  6:246964/91‑1 (MQ=255)                                                                                                                tcttctGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTg                                                                                          >  5:246964/1‑91 (MQ=255)                                                                                                                    ctgattctgatCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGTGCAAAGCTGTCGCTACGCGCCATCCGCCACATCCGCGCATaaaa                             <  2:211418/148‑1 (MQ=255)                                                                                                                       attctgatCTGATGTATCGAACAGCCCCGGCCAATGCCCGGGCTGGTTGGTTTCGGCGATCGGTTACTGTTTGAGCCGGCGCAGGGCGTACGGTCTCCCGTGCAAACCTGTCGCTACGCGCCATCCCCCCCATACGCGCATAAAACTCg                         >  8:102689/1‑149 (MQ=255)                                                                                                                           tgatctgatGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGTGCAAAGCTGTCGCTACGCGCCATCCGCCACATCCGCGCATAAAACTCGCttt                     >  3:84537/1‑149 (MQ=255)                                                                                                                                  atGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGTGCAAAGCTGTCGCTACGCGCCATCCGCCACATCCGCGCATAAAACTCGCTTtcgatc                >  2:127206/1‑147 (MQ=255)                                                                                                                                     tATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGTGCAAAGCTGTCGCTACGCGCCATCCGCCACATCCGCGCATAAAACTCGCTTTCGATCGAGc            >  1:66263/1‑148 (MQ=255)                                                                                                                                         gAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGTGCAAAGCTGTCGCTACGCGCCATCCGCCACATCCGCGCATAAAACTCGCTTTCGATCGAGCCACt        <  3:17461/148‑1 (MQ=255)                                                                                                                                          aaCAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGTGCAAAGCTGTCGCTACGCGCCATCCGCCACATCCGCGCATAAAACTCGCTTTCGATCGAGCCACTGa      <  6:130424/149‑1 (MQ=255)                                                                                                                                          aaCAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGTGCAAAGCTGTCGCTACGCGCCATCCGCCACATCCGCGCATAAAACTCGCTTTCGATCGAGCCACTGa      <  6:226341/149‑1 (MQ=255)                                                                                                                                          aaCAGCCCCGGCATAGGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGTGCAAAGCTGTCGCTACGCGCCATCCGCCACATCCGCGCATAAAACTCGCTTTCGATCGAGCCACTGa      >  6:70180/1‑149 (MQ=255)                                                                                                                                            cAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGTGCAAAGCTGTCGCTACGCGCCATc                                                 <  5:199976/104‑1 (MQ=255)                                                                                                                                            cAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGTGCAAAGCTGTCGCTACGCGCCATc                                                 >  6:199976/1‑104 (MQ=255)                                                                                                                                              gCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGTGCAAAGCTGTCGCTACGCGCCATCCGCCACATCCGCGCATAAAACTCGCTTTCGATCGAGCCACTGAGTaa  >  5:149359/1‑149 (MQ=255)                                                                                                                                                 ccGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGTGCAAAGCTGTCGCTACGCGCCATCCGCCACATCCGCGCATAAAACTCGCTTtcgatc                >  7:76488/1‑132 (MQ=255)                                                                                                                                                 ccGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGTGCAAAGCTGTCGCTACGCGCCATCCGCCACATCCGCGCATAAAACTCGCTTtcgatc                <  8:76488/132‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                             GTCGATCAAGCTGCGTAGCTCGTTCTTGCGTACCAACAACGCCGTGCAGCAGAACGGTTACAACGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCC‑GATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGTGCAAAGCTGTCGCTACGCGCCATCCGCCACATCCGCGCATAAAACTCGCTTTCGATCGAGCCACTGAGTAA  >  NC_002947/325917‑326205

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