Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,966,762	C→A	6.1%	L319M (CTG→ATG) ‡	ydcP →	peptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,966,762	0	C	A	6.1%	211.9 / 2.8	65	L319M (CTG→ATG) ‡	ydcP	peptidase
Reads supporting (aligned to +/- strand):  ref base C (25/36);  new base A (2/2);  total (27/38)
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.30e-01

CCCTCAGTACGCCCGCGCGTCCAGTGGCCGTACCCAACACTTCTTCCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCG‑‑AC‑CGCAAGGTCGATAT‑CGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGACGGCCTCAATGTGCTGGTCAAGCGCGAAGTGGTGGGCTTCCGCGCCAACATCGCCGAACCGCGTGGCGAA  >  NC_002947/3966614‑3966908                                                                                                                                                         |                                                                                                                                                   cccTCAGTACGCCCGCGCGTCCAGTGGCCGTACCCAACACTTCTTCCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCG‑‑AC‑CGCAAGGTCGATAT‑CGGCGCCTTCGACTCGCCAACCTTCACCGGcc                                                                                                                                                    >  1:256524/1‑149 (MQ=255) cccTCAGTACGCCCGCGCGTCCAGTGGCCGTACCCAACACTTCTTCCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCG‑‑AC‑CGCAAGGTCGATAT‑CGGCGCCTTCGACTCGCCAACCTTCACCGGcc                                                                                                                                                    >  3:38266/1‑149 (MQ=255)   cTCAGTACGCCCGCGCGTCCAGTGGCCGTACCCAACACTTCTTCCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCG‑‑AG‑CGCAAGGTCGATAT‑CGGCGCCTTCGACTCGCCAACCTTCACCGGCCTg                                                                                                                                                  <  5:277335/149‑1 (MQ=255)     cAGTACGCCCGCGCGTCCAGTGGCCGTACCCAACACTTCTTCCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCG‑‑AC‑CGCAAGGTCGATAT‑CGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGcc                                                                                                                                                >  7:77677/1‑149 (MQ=255)          cGCCCGCGCGTCCAGTGGCCGTCCCCAACACTTCTTCCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCG‑‑AC‑CGCAAGGTCGATAT‑CGGCGCCTTCGACTCGCCAACCTTCACCGGCAGGCCGGTga                                                                                                                                           <  5:277497/149‑2 (MQ=255)          cGCCCGCGCGTCCAGTGGCCGTACCCAACACTTCTTCCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCG‑‑AC‑CGCAAGGTCGATAT‑CGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTgg                                                                                                                                           <  7:205484/149‑1 (MQ=255)          cGCCCGCGCGTCCAGTGGCCGTACCCAACACTTCTTCCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCG‑‑AC‑CGCAAGGTCGATAT‑CGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTgg                                                                                                                                           <  4:117264/149‑1 (MQ=255)          cGCCCGCGCGTCCAGTGGCCGTACCCAACACTTCTTCCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCG‑‑AC‑CGCAAGGTCGATAT‑CGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTg                                                                                                                                            <  1:365890/148‑1 (MQ=255)          cGCCCGCGCGTCCAGTGGCCGTACCCAACACTTCTTCCTGCCCGACCCGGATAAAACCTACCACCGCGGCAGCACCGACTACTTCGTGACCG‑‑AC‑CGCAAGGTCGATAT‑CGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTgg                                                                                                                                           <  7:347922/149‑1 (MQ=255)            cccGCGCGTCCAGTGGCCGTACCCAACACTTCTTCCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCG‑‑AC‑CGCAAGGTCGATAT‑CGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGc                                                                                                                                         <  8:77677/149‑1 (MQ=255)             ccGCGCGTCCAGTGGCCGTACCCAACACTTCTTCCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCG‑‑AC‑CGCAAGGTCGATATCCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGc                                                                                                                                         <  1:71306/149‑1 (MQ=255)             ccGCGCGTCCAGTGGCCGTACCCAACACTTCTTCCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCG‑‑AC‑CGCAAGGTCGATAT‑CGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTgg                                                                                                                                           >  5:123359/1‑146 (MQ=255)             ccGCGCGTCCAGTGGCCGTACCCAACACTTCTTCCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCG‑‑AC‑CGCAAGGTCGATAT‑CGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTgg                                                                                                                                           <  6:123359/146‑1 (MQ=255)               gcgcgTCCAGTGGCCGTACCCAACACTTCTTCCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCG‑‑AC‑CGCAAGGTCGATAT‑CGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTa                                                                                                                                      <  6:276685/149‑1 (MQ=255)                cgcgTCCAGTGGCCGTACCCAACACTTCTTCCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCG‑‑AC‑CGCAAGGTCGATAT‑CGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAg                                                                                                                                     <  2:256524/149‑1 (MQ=255)                 gcgTCCAGTGGCCGTACCCAACACTTCTTCCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCG‑‑AC‑CGCAAGGTCGATAT‑CGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGt                                                                                                                                    <  7:260862/149‑1 (MQ=255)                 gcgTCCAGTGGCCGTACCCAACACTTCTTCCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCG‑‑AC‑CGCAAGGTCGATAT‑CGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGt                                                                                                                                    <  6:267937/149‑1 (MQ=255)                       aGTGGCCGTACCCAACACTTCTTCCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCG‑‑AC‑CGCAAGGTCGATAT‑CGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGaaaa                                                                                                                              <  2:188079/149‑1 (MQ=255)                       aGTGGCCGTACCCAACACTTCTTCCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCG‑‑AC‑CGCAAGGTCGATAT‑CGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGaa                                                                                                                                >  2:57678/1‑147 (MQ=255)                         tGGCCGTACCCAACACTTCTTCCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCG‑‑AC‑CGCAAGGTCGATAT‑CGGCGCCTTCGACTCGCTAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAgg                                                                                                                            >  8:251237/1‑149 (MQ=255)                           gCCGTACCCAACACTTCTTCCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCGCCACTCGCAATGTCGATAT‑CGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGaaaa                                                                                                                              <  1:300478/148‑1 (MQ=255)                           gCCGTACCCAACACTTCTTCCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCG‑‑AC‑CGCAAGGTCGATAT‑CGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTc                                                                                                                          <  8:184688/149‑1 (MQ=255)                                 cccAACACTTCTTCCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCG‑‑AC‑CGCAAGGTCGATAT‑CGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAg                                                                                                                    >  6:340735/1‑149 (MQ=255)                                  ccAACACTTCTTCCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCG‑‑AC‑CGCAAGGTCGATAT‑CGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGc                                                                                                                   >  5:315699/1‑149 (MQ=255)                                  ccAACACTTCTTCCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCG‑‑AC‑CGCAAGGTCGATAT‑CGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGc                                                                                                                   >  4:162560/1‑149 (MQ=255)                                    aaCACTTCTTCCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCG‑‑AC‑CGCAAGGTCGATAT‑CGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGt                                                                                                                 >  4:5125/1‑149 (MQ=255)                                     acacTTCTTCCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCG‑‑AC‑CGCAAGGTCGATAT‑CGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTg                                                                                                                <  1:57678/149‑1 (MQ=255)                                       acTTCTTCCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCG‑‑AC‑CGCAAGGTCGATAT‑CGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGa                                                                                                               >  4:56167/1‑148 (MQ=255)                                       acTTCTTCCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCG‑‑AC‑CGCAAGGTCGATAT‑CGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGAc                                                                                                              <  4:213992/149‑1 (MQ=255)                                         ttcttcCTGCCCTACCGGGATAAAACCTTCCACCGCGGAAGCACCGACTACTTCGTGACCG‑‑AC‑CGCAAGGTGGATAT‑CGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTc                                                                                                                          <  7:286344/135‑1 (MQ=255)                                          tcttcCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCG‑‑AC‑CGCAAGGTCGATAT‑CGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTgg                                                                                                                                           >  3:121221/1‑117 (MQ=255)                                          tcttcCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCG‑‑AC‑CGCAAGGTCGATAT‑CGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTgg                                                                                                                                           <  4:121221/117‑1 (MQ=255)                                           cttcCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCG‑‑AC‑CGCAAGGTCGATAT‑CGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATg                                                                                                           >  2:202548/1‑148 (MQ=255)                                            ttcCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCG‑‑AC‑CGCAAGGTCGATAT‑CGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGc                                                                                                          >  6:10468/1‑148 (MQ=255)                                              ccTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCG‑‑AC‑CGCAAGGTCGATAT‑CGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCaa                                                                                                        <  3:56167/148‑1 (MQ=255)                                                  cccGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCG‑‑AC‑CGCAAGGTCGATAT‑CGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGt                                                                                                   <  8:29755/149‑1 (MQ=255)                                                  cccGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCG‑‑AC‑CGCAAGGTCGATAT‑CGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGt                                                                                                   <  2:351501/149‑1 (MQ=255)                                                     gACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCG‑‑AC‑CGCAAGGTCGATAT‑CGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCAc                                                                                                <  3:5125/149‑1 (MQ=255)                                                       cccGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCG‑‑AC‑CGCAAGGTCGATAT‑CGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCAcc                                                                                               >  8:287355/1‑148 (MQ=255)                                                        ccGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCG‑‑AC‑CGCAAGGTCGATAT‑CGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTgg                                                                                                                                           <  5:195096/103‑1 (MQ=255)                                                        ccGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCG‑‑AC‑CGCAAGGTCGATAT‑CGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTgg                                                                                                                                           >  6:195096/1‑103 (MQ=255)                                                        ccGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCG‑‑AC‑CGCAAGGTCGATAT‑CGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGa                                                                                             <  6:315699/149‑1 (MQ=255)                                                         cGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCG‑‑AC‑CGCAAGGTCGATAT‑CGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGAc                                                                                            >  6:121735/1‑149 (MQ=255)                                                             tAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCG‑‑AC‑CGCAAGGTCGATAT‑CGGCGCCTTCGACTCGCTAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTAc                                                                                        <  7:251237/149‑1 (MQ=255)                                                             tAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCG‑‑AC‑CGCAAGGTCGATAT‑CGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTAc                                                                                        >  7:36314/1‑149 (MQ=255)                                                                    ttCCACCGCGGCAGCACCGACTACTTCGTGACCG‑‑AC‑CGCAAGGTCGATAT‑CGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGAc                                                                                 >  6:49671/1‑149 (MQ=255)                                                                      ccACCGCGGCAGCACCGACTACTTCGTGACCG‑‑AC‑CGCAAGGTCGATAT‑CGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGAc                                                                                 <  5:49671/147‑1 (MQ=255)                                                                              gcagcaCCGACTACTTCGTGACCG‑‑AC‑CGCAAGGTCGATAT‑CGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGAc                                                                       <  5:145178/149‑1 (MQ=255)                                                                               cagcaCCGACTACTTCGTGACCG‑‑AC‑CGCAAGGTCGATAT‑CGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCTTAGTGGAAAAGGTCGGCAAGCGTGACATGCAATTGGTCACCGACGTACAGCTGACCAACGGCGACg                                                                      >  8:149439/1‑149 (MQ=255)                                                                                   aCCGACTACTTCGTGACCG‑‑AC‑CGCAAGGTCGATAT‑CGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGACGGCCt                                                                  <  1:202548/149‑1 (MQ=255)                                                                                   aCCGACTACTTCGTGACCG‑‑AC‑CGCAAGGTCGATAT‑CGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGACGGCCt                                                                  <  7:287355/149‑1 (MQ=255)                                                                                       actactTCGTGACCG‑‑AC‑CGCAAGGTCGATAT‑CGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGACGGCCTCAAt                                                              >  4:109141/1‑149 (MQ=255)                                                                                                   ccg‑‑ac‑cgCAAGGTCGATAT‑CGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGACGGCCTCAATGTGCTGGTCAAg                                                  <  2:82475/149‑1 (MQ=255)                                                                                                    cg‑‑ac‑cgCAAGGTCGATAT‑CGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACg                                                                            >  1:216551/1‑122 (MQ=255)                                                                                                    cg‑‑ac‑cgCAAGGTCGATAT‑CGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACg                                                                            <  2:216551/122‑1 (MQ=255)                                                                                                            gCAAGGTCGATAT‑CGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGACGGCCTCAATGTGCTGGTCAAGCGCGaa                                            <  5:340735/149‑1 (MQ=255)                                                                                                                     atat‑CGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGACGGCCTCAATGTGCTGGTCAAGCGCGAAGTGGTGGGc                                   <  5:121735/149‑1 (MQ=255)                                                                                                                           ggCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGACGGCCTCAATGTGCTGGTCAAGCGCGAAGTGGTGGGCTTCcg                              <  3:162560/149‑1 (MQ=255)                                                                                                                                         gCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGACGGCCTCAATGTGCTGGTCAAGCGCGAAGTGGTGGGCTTCCGCGCCAACATCGcc                 >  4:117734/1‑148 (MQ=255)                                                                                                                                         gCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGACGGCCTCAATGTGCTGGTCAAGCGCGAAGTGGTGGGCTTCCGCGCCAACATCGcc                 >  2:26246/1‑148 (MQ=255)                                                                                                                                          ccAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCg                                                                                                                         <  4:119049/43‑1 (MQ=255)                                                                                                                                          ccAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGACGGCCTCAATGTGCTGGTCAAGCGCGAAGTCGTGGGCTTCCGCGCCAACATCGCCg                >  5:280965/1‑148 (MQ=255)                                                                                                                                             aCCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGACGGCCTCAATGTGCTGGTCAAGCGCGAAGTGGTGGGCTTCCGCGCCAACATCGCCGAAcc            >  7:49819/1‑149 (MQ=255)                                                                                                                                               gcccACCGGCAGGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGACGGCCTCAATGTGCTGGTCAAGCGCGAAGTGGTGGGCTTCCGCGCCAACATCGCCGAACcgc          >  6:148470/4‑149 (MQ=255)                                                                                                                                                  cACCGGCAGGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGACGGCCTCAATGTGCTGGTCAAGCGCGAAGTGGTGGGCTTCCGCGCCAACATCGc                  >  7:191770/1‑138 (MQ=255)                                                                                                                                                  cACCGGCAGGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGACGGCCTCAATGTGCTGGTCAAGCGCGAAGTGGTGGGCTTCCGCGCCAACATCGc                  <  8:191770/138‑1 (MQ=255)                                                                                                                                                       gCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGACGGCCTCAATGTGCTGGTCAAGCGCGAAGTGGTGGGCTTCCGCGCCAACATCGCCGAACCGCGTGGCGaa  >  4:36883/1‑149 (MQ=255)                                                                                                                                                         |                                                                                                                                                   CCCTCAGTACGCCCGCGCGTCCAGTGGCCGTACCCAACACTTCTTCCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCG‑‑AC‑CGCAAGGTCGATAT‑CGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGACGGCCTCAATGTGCTGGTCAAGCGCGAAGTGGTGGGCTTCCGCGCCAACATCGCCGAACCGCGTGGCGAA  >  NC_002947/3966614‑3966908

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