Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,632,086	A→C	5.7%	E327A (GAG→GCG)	PP_5590 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,632,086	0	A	C	5.7%	206.0 / 2.3	66	E327A (GAG→GCG)	PP_5590	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (28/34);  new base C (2/2);  total (30/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 = 1.00e+00

GGCCACCGCCGGCGTCATGAGCTTCGATGAGCTGTACGTGTTCAACGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCC‑ATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATGCCGCGTAGCACCTACTACCTCGACGAAGACACC  >  NC_002947/3631946‑3632225                                                                                                                                             |                                                                                                                                             ggCCACCGCCGGCGTCATGAGCTTCGATGAGCTGTACGTGTTCAACGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAGCGCCTCGGCgttgc                                                                                                                                      <  4:89172/149‑1 (MQ=255)     aCCGCCGGCGTCATGAGCTTCGATGAGCTGTACGTGTTCAACGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAGCGCCTCGGCGTTGCTGGc                                                                                                                                  <  2:114071/149‑1 (MQ=255)     aCCGCCGGCGTCATGAGCTTCGATGAGCTGTACGTGTTCAACGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGc                                                                                                                                  >  6:135526/1‑149 (MQ=255)              gTCATGAACTTCGATGAGCTGTACGTGTTCAACGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCC‑ATCa                                                                                                                        <  6:282678/149‑1 (MQ=255)                 aTGAGCTTCGATGAGCTGTACGTGTTCAACGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCC‑ATCACCt                                                                                                                     <  5:325944/149‑1 (MQ=255)                         cGATGAGCTGTACGTGTTCAACGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCC‑ATCACCTGAAcccc                                                                                                              <  8:6363/148‑1 (MQ=255)                                  gtcggtgtTAAACGGTCGCCTCGACCGCTGCGAGTGGAAGGTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGAGTTAACCGCCCTCAACCCCCCGGCGCTGCTGGCCCCCC‑AGCACCTGCACCCCGAAGCCATcc                                                                                                    >  2:216254/5‑149 (MQ=255)                                      gtgtTCAACGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTTGGCAAGAAGGAGATGTACATCCCCTAAAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGGTGCTGGCCGCCC‑ATCACCTGAACCCCGACGCCATCCGCTg                                                                                                <  1:216254/149‑1 (MQ=255)                                      gtgtTCAACGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCC‑ATCACCTGAACCCCGACGCCATCCGCTg                                                                                                >  1:298305/1‑149 (MQ=255)                                       tgtTCAACGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCC‑ATCACCTGAACCCCGACGCCATCCGCTgg                                                                                               <  6:345365/149‑1 (MQ=255)                                            aaCGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCC‑ATCACCTGAACCCCGACGCCATCCGCTGGGAGCt                                                                                          >  7:36580/1‑149 (MQ=255)                                                    ccTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCC‑ATCACCTGAACCCCGa                                                                                                            <  8:306035/123‑1 (MQ=255)                                                    ccTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCC‑ATCACCTGAACCCCGa                                                                                                            >  7:306035/1‑123 (MQ=255)                                                    ccTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCC‑ATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCgtg                                                                                  <  4:5818/149‑1 (MQ=255)                                                    ccTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCC‑ATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCgtg                                                                                  <  4:1150/149‑1 (MQ=255)                                                      tCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCC‑ATCGCCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCgtgtg                                                                                >  4:169512/1‑149 (MQ=255)                                                      tCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCC‑ATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCgtgtg                                                                                >  6:278665/1‑149 (MQ=255)                                                        gACCGCTTCGAGTGGTAGCTGGTGGGCAAGAATGAGATGTACATCCGCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCC‑ATCACCTGAACCCCGACGCCATCCGCTGGGAGCTACACCgtgtgt                                                                               <  7:130673/148‑1 (MQ=255)                                                        gACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGGGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCC‑ATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGtgtgtg                                                                              <  3:57665/149‑1 (MQ=255)                                                          ccGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCC‑ATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTgtggg                                                                            >  3:45018/1‑149 (MQ=255)                                                             cTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCC‑ATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGc                                                                                        >  3:122711/1‑134 (MQ=255)                                                             cTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCC‑ATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGc                                                                                        <  4:122711/134‑1 (MQ=255)                                                              ttCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCC‑ATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGggtggt                                                                        >  2:246839/1‑149 (MQ=255)                                                                cGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCC‑ATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGgtggtg                                                                       <  8:36580/148‑1 (MQ=255)                                                                  aGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCC‑ATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAg                                                                    <  1:277625/149‑1 (MQ=255)                                                                   gTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCC‑ATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAgg                                                                   >  4:313618/1‑149 (MQ=255)                                                                    tGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCC‑ATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGc                                                                  <  7:180325/149‑1 (MQ=255)                                                                     ggAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCGCC‑ATCAACTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGc                                                                  >  3:66024/1‑148 (MQ=255)                                                                     ggAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCC‑ATCACCTGAACCCCGAAGCCATCCGCTGGGAGCTGGACCGTGTGTGGGTGGTGGAGGc                                                                  >  3:70527/1‑148 (MQ=255)                                                                         gCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCC‑ATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGc                                                             <  5:278665/149‑1 (MQ=255)                                                                         gCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCC‑ATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGc                                                             <  5:209382/149‑1 (MQ=255)                                                                         gCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCC‑ATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGc                                                             >  5:202889/1‑149 (MQ=255)                                                                             gtggGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCC‑ATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTg                                                           >  8:11423/1‑147 (MQ=255)                                                                              ttggCAAGAATGAGATGTACATCCACTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCC‑ATCACCTGAACCCCGACCCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTgaag                                                        <  4:70527/147‑1 (MQ=255)                                                                              ttggCAAGAAGTATATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGAAACGCCGAGGCGCTGCTGGCCCCCC‑ATCACCTGAACCCCGACGCCATCCGCTGGTAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTgaag                                                        <  4:66024/147‑1 (MQ=255)                                                                                   aagaagGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCC‑ATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATgg                                                   <  4:45018/149‑1 (MQ=255)                                                                                             tGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCC‑ATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGACCAGCGCCAt                                         <  2:298305/149‑1 (MQ=255)                                                                                              gTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCC‑ATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTgaa                                                         >  6:136629/1‑132 (MQ=255)                                                                                              gTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCC‑ATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTgaa                                                         <  5:136629/132‑1 (MQ=255)                                                                                               tACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCC‑ATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCAtgt                                       <  5:135526/149‑1 (MQ=255)                                                                                                     cccTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCC‑ATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATTcc                                 <  7:11423/149‑1 (MQ=255)                                                                                                      ccTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCC‑ATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGAt                                    <  6:32859/145‑1 (MQ=255)                                                                                                      ccTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCC‑ATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGAt                                    >  5:32859/1‑145 (MQ=255)                                                                                                      ccTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCC‑ATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATGCcg                                >  5:250106/1‑149 (MQ=255)                                                                                                      ccTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCC‑ATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATGCcg                                >  1:197530/1‑149 (MQ=255)                                                                                                        tACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCC‑ATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATGc                                  >  4:86365/1‑145 (MQ=255)                                                                                                        tACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCC‑ATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATGc                                  <  3:86365/145‑1 (MQ=255)                                                                                                          cAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCC‑ATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATGCcgcg                              >  7:247947/1‑147 (MQ=255)                                                                                                            aCGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCC‑ATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATGCCGCGTAGc                          <  4:146587/149‑1 (MQ=255)                                                                                                             cGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCC‑ATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATGCCGCGTAGCa                         >  5:315378/1‑149 (MQ=255)                                                                                                              gCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCC‑ATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATGCCGCGTAGCAc                        >  3:116324/1‑149 (MQ=255)                                                                                                                 aacaaGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCC‑ATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATGCCGCGTAGCACcta                     <  4:91151/149‑1 (MQ=255)                                                                                                                 aacaaGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCC‑ATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATGCCGCGTAGCACcta                     <  2:197530/149‑1 (MQ=255)                                                                                                                 aacaaGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCC‑ATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATGCCGCGTAGCACcta                     >  1:246767/1‑149 (MQ=255)                                                                                                                  acaaGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCC‑ATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGGGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATGCCGCGTAGCACctac                    <  5:113835/149‑1 (MQ=255)                                                                                                                              ccGCCAGCAACGCCGAGGCGCTGCTGGCCCCCC‑ATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATGCCGCATAGCACCTACTACCTCGACGaa        <  1:209321/149‑1 (MQ=255)                                                                                                                                gCCAGCAGCGCCTCGGCGTTGCTGGCCCCCCAATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATGCCGCGTAGCACCTACTACCTcgacg          >  1:36339/1‑146 (MQ=255)                                                                                                                                  cAGCAACGCCGAGGCGCTGCTGGCCCCCC‑ATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATGCCGCGTAGCACCTACTACCTCGACGAAGaca    <  5:212767/149‑1 (MQ=255)                                                                                                                                  aaGCAACGCCGAGGCGCTGCTGGCCCCCC‑ATAACCTGACCACCGACGCCATCCGCTGGGAGATGCACCGTGTGTGGGTGGTGGAGGCGCCGCTGAAGAATGGCCAGCGCGAGGTGATGCCGCGTATCAACTACTACCGCGACGAAGaca    >  5:126811/2‑149 (MQ=255)                                                                                                                                    gCAACGCCGAGGCGCTGCTGGCCCCCC‑ATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATGCCGCGTAGCACCTACTACCTCGACGAAGACAcc  >  3:293911/1‑149 (MQ=255)                                                                                                                                      aaCGCCGAGGCGCTGCTGGCCCCCC‑ATCACCTGAACCCCGACGCCATCCGCTGGGAg                                                                                            <  7:257565/57‑1 (MQ=255)                                                                                                                                      aaCGCCGAGGCGCTGCTGGCCCCCC‑ATCACCTGAACCCCGACGCCATCCGCTGGGAg                                                                                            >  8:257565/1‑57 (MQ=255)                                                                                                                                        cGCCGAGGCGCTGCTGGCCCCCC‑ATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATGCCGCGTAGCACcta                     >  7:133838/1‑126 (MQ=255)                                                                                                                                        cGCCGAGGCGCTGCTGGCCCCCC‑ATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATGCCGCGTAGCACcta                     <  8:133838/126‑1 (MQ=255)                                                                                                                                           cGAGGCGCTGCTGGCCCCCC‑ATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATGCCGCGTAGCACCTACTACCTCGACGaa        <  2:31116/136‑1 (MQ=255)                                                                                                                                           cGAGGCGCTGCTGGCCCCCC‑ATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATGCCGCGTAGCACCTACTACCTCGACGaa        >  1:31116/1‑136 (MQ=255)                                                                                                                                             |                                                                                                                                             GGCCACCGCCGGCGTCATGAGCTTCGATGAGCTGTACGTGTTCAACGGTCGCCTCGACCGCTTCGAGTGGAAGCTGGTGGGCAAGAAGGAGATGTACATCCCCTACAACGCCAACAAGGTGTTCACCGCCAGCAACGCCGAGGCGCTGCTGGCCCCCC‑ATCACCTGAACCCCGACGCCATCCGCTGGGAGCTGCACCGTGTGTGGGTGGTGGAGGCGACGCTGAAGAATGGCCAGCGCCATGTGATGCCGCGTAGCACCTACTACCTCGACGAAGACACC  >  NC_002947/3631946‑3632225

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