Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,966,763	T→G	6.1%	L319R (CTG→CGG) ‡	ydcP →	peptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,966,763	0	T	G	6.1%	208.0 / 4.1	65	L319R (CTG→CGG) ‡	ydcP	peptidase
Reads supporting (aligned to +/- strand):  ref base T (25/36);  new base G (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.68e-01

CTCAGTACGCCCGCGCGTCCAGTGGCCGTACCCAACACTTCTTCCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCG‑‑AC‑CGCAAGGTCGATAT‑CGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGACGGCCTCAATGTGCTGGTCAAGCGCGAAGTGGTGGGCTTCCGCGCCAACATCGCCGAACCGCGTGGCGAA  >  NC_002947/3966616‑3966908                                                                                                                                                        |                                                                                                                                                  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                                                                                                                                           <  4:121221/117‑1 (MQ=255)                                        tcttcCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCG‑‑AC‑CGCAAGGTCGATAT‑CGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTgg                                                                                                                                           >  3:121221/1‑117 (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                                                                  <  7:287355/149‑1 (MQ=255)                                                                                 aCCGACTACTTCGTGACCG‑‑AC‑CGCAAGGTCGATAT‑CGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGACGGCCt                                                                  <  1:202548/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)                                                                                                                                                        |                                                                                                                                                  CTCAGTACGCCCGCGCGTCCAGTGGCCGTACCCAACACTTCTTCCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCG‑‑AC‑CGCAAGGTCGATAT‑CGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGACGGCCTCAATGTGCTGGTCAAGCGCGAAGTGGTGGGCTTCCGCGCCAACATCGCCGAACCGCGTGGCGAA  >  NC_002947/3966616‑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: ‑