Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,031,920	T→A	6.1%	L5* (TTG→TAG) ‡	PP_0889 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,031,920	0	T	A	6.1%	203.5 / 2.7	65	L5* (TTG→TAG) ‡	PP_0889	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (27/34);  new base A (2/2);  total (29/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

AATATCTCGGCCACCGCACGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGACCCGCAGTGGCAGAACGACCCTGCCACCCTCGACTGG  >  NC_002947/1031772‑1032065                                                                                                                                                     |                                                                                                                                                   aaTATCTCGGCCACCGCACGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCtt                                                                                                                                                    <  8:256405/149‑1 (MQ=255) aaTATCTCGGCCACCGCACGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCtt                                                                                                                                                    >  8:136115/1‑149 (MQ=255)  atatCTCGGCCACCGCACGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCATCTGTGAGCGCGGGCTTg                                                                                                                                                   >  3:310505/1‑149 (MQ=255)     tctcGGCCACCGCACGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCAAGccc                                                                                                                                                >  4:308916/1‑149 (MQ=255)       tcGGCCACCGCACGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCcg‑c                                                                                                                                             >  6:185528/1‑149 (MQ=255)         ggCCACCGCACGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCAAGCCcg‑cgc                                                                                                                                           <  3:308916/149‑2 (MQ=255)           ccACCGCACGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCcg‑c                                                                                                                                             <  1:312932/145‑1 (MQ=255)           ccACCGCACGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCcg‑c                                                                                                                                             >  2:312932/1‑145 (MQ=255)           ccACCGCACGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAg                                                                                                                                         <  5:185528/149‑1 (MQ=255)              ccGCACGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGc                                                                                                                                      >  5:250448/1‑149 (MQ=255)                 cACGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTg                                                                                                                                    >  7:331706/1‑148 (MQ=255)                  aCGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCa                                                                                                                                  <  8:302267/149‑1 (MQ=255)                                  cACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCaa                                                                                                                  <  7:136115/149‑1 (MQ=255)                                   aCCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAg                                                                                                                 >  6:126663/1‑149 (MQ=255)                                    ccGACGGACGTTGCAGCGGAAGCTGCAGAAGCGGCCGGTTCGGCGCTGTTTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGGGCGCGCGGGCTTGGCCGACGAAGGGCTGCAAAGCAGTCCCAATCAAg                                                                                                                 <  1:26453/149‑1 (MQ=255)                                     cGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCAAGCCcg‑cgc                                                                                                                                           <  6:22524/121‑2 (MQ=255)                                     cGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCAAGCCcg‑cgc                                                                                                                                           >  5:22524/1‑120 (MQ=255)                                         ggacgTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGc                                                                                                           <  7:61893/149‑1 (MQ=255)                                             gTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGa                                                                                                       <  8:331706/149‑1 (MQ=255)                                                gCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATg                                                                                                     <  6:250448/148‑1 (MQ=255)                                                  aGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCATCTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAg                                                                                                  <  4:310505/149‑1 (MQ=255)                                                           tGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTCTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTg                                                                                         >  1:110550/1‑149 (MQ=255)                                                             cAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTg                                                                                         <  7:168815/147‑1 (MQ=255)                                                             cAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTg                                                                                         >  8:168815/1‑147 (MQ=255)                                                                    gtccGGTTCGGCGCTTATTCCTTATCGTTGATGTTGGCCTGTTCGCTGGGGTATCCGCTCCTGCAACTGTTAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTc                                                                                                         <  3:134296/122‑1 (MQ=38)                                                                     gCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAg                                                                                                                               >  1:36890/1‑101 (MQ=255)                                                                     gCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAg                                                                                                                               <  2:36890/101‑1 (MQ=255)                                                                     gCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGa                                                                                                             >  6:57627/1‑119 (MQ=255)                                                                     gCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGa                                                                                                             <  5:57627/119‑1 (MQ=255)                                                                     gCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGa                                                                                                    <  1:347953/128‑1 (MQ=255)                                                                     gCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGa                                                                                                    >  2:347953/1‑128 (MQ=255)                                                                      ccGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTgc                                                                              <  5:126663/149‑1 (MQ=255)                                                                             ggCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGtt                                                                       >  6:345355/1‑149 (MQ=255)                                                                               cgcTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGc                                                                     <  6:155632/149‑1 (MQ=255)                                                                                 cTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAg                                                                   <  7:324980/149‑1 (MQ=255)                                                                                 cTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAg                                                                   <  7:320334/149‑1 (MQ=255)                                                                                 cTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAg                                                                   <  6:171185/149‑1 (MQ=255)                                                                                         ttATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATc                                                                                                                    <  6:215048/92‑1 (MQ=255)                                                                                         ttATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATc                                                                                                                    >  5:215048/1‑92 (MQ=255)                                                                                         ttATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTg                                                                                         <  3:146655/119‑1 (MQ=255)                                                                                         ttATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTg                                                                                         >  4:146655/1‑119 (MQ=255)                                                                                              gatgatGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTg                                                                                         <  2:98463/114‑1 (MQ=255)                                                                                              gatgatGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTg                                                                                         >  1:98463/1‑114 (MQ=255)                                                                                              gatgatGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGc                                                      >  8:165694/1‑149 (MQ=255)                                                                                              gatgatGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGc                                                      >  7:295907/1‑149 (MQ=255)                                                                                               atgatgTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTgc                                                                              >  6:48998/1‑124 (MQ=255)                                                                                               atgatgTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTgc                                                                              <  5:48998/124‑1 (MQ=255)                                                                                               atgatgTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCt                                                     <  8:295907/149‑1 (MQ=255)                                                                                                     ttGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCACCAAAGCTGAAGAc                                               >  7:95822/1‑149 (MQ=255)                                                                                                        gCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTg                                             <  2:110550/148‑1 (MQ=255)                                                                                                         ccTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAgcgc                                                            <  4:103885/132‑1 (MQ=255)                                                                                                         ccTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAgcgc                                                            >  3:103885/1‑132 (MQ=255)                                                                                                               cgcgGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGAc                                      <  5:10215/148‑1 (MQ=255)                                                                                                                   ggggAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGaa                                                  >  4:129321/1‑132 (MQ=255)                                                                                                                   ggggAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGaa                                                  <  3:129321/132‑1 (MQ=255)                                                                                                                       aaTCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGACCCGCAGTgg                             <  7:165694/149‑1 (MQ=255)                                                                                                                            gCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGa                                                >  2:120154/1‑125 (MQ=255)                                                                                                                            gCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGa                                                <  1:120154/125‑1 (MQ=255)                                                                                                                                cTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCaaa                                                        <  4:279314/113‑1 (MQ=255)                                                                                                                                cTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCaaa                                                        >  3:279314/1‑113 (MQ=255)                                                                                                                                cTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGACCCGCAGTGGCAGAACGAc                    <  5:345355/149‑1 (MQ=255)                                                                                                                                       tgtgAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGc                                           >  5:194288/1‑119 (MQ=255)                                                                                                                                       tgtgAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGc                                           <  6:194288/119‑1 (MQ=255)                                                                                                                                       tgtgAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGc                                           >  1:136663/1‑119 (MQ=255)                                                                                                                                       tgtgAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGc                                           <  2:136663/119‑1 (MQ=255)                                                                                                                                            gcgcGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGACCCGAAGTGGCAGAACGACCCTGCCACCCTc        >  6:363431/1‑149 (MQ=255)                                                                                                                                                  gCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGACCCGCAGTGGCAGAACGACCCTGCCACCCTCGACTgg  >  6:358241/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                   AATATCTCGGCCACCGCACGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGACCCGCAGTGGCAGAACGACCCTGCCACCCTCGACTGG  >  NC_002947/1031772‑1032065

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