Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,966,806	C→A	7.4%	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,806	0	C	A	7.4%	298.4 / 8.4	95	L319M (CTG→ATG) ‡	ydcP	peptidase
Reads supporting (aligned to +/- strand):  ref base C (44/44);  new base A (3/4);  total (47/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.98e-01

CCCTCAGTACGCCCGCGCGTCCAGTGGCCGTACCCAACACTTCTTCCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGACGGCCTCAATGTGCTGGTCAAGCGCGAAGTGGTGGGCTTCCGCGCCAACATCGCCGAACCGCGTGGCGAATTC  >  NC_002947/3966658‑3966955                                                                                                                                                     |                                                                                                                                                      cccTCAGTACGCCCGCGCGTCCAGTGGCCGTACCCAACACTTCTTCCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCt                                                                                                                                                      >  1:1574796/1‑150 (MQ=255) cccTCAGTACGCCCGCGCGTCCAGTGGCCGTACCCAACACTTCTTCCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCt                                                                                                                                                      >  1:1404057/1‑150 (MQ=255)       gTACGCCCGCGCGTCCAGTGGCCGTACCCAACACTTCTTCCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGt                                                                                                                                                <  1:1801320/150‑1 (MQ=255)         aCGCCCGCGCGTCCAGTGGCCGTACCCAACACTTCTTCCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTgg                                                                                                                                              <  1:2064674/150‑1 (MQ=255)           gCCCGCGCGTCCAGTGGCCGTACCCAACACTTCTTCCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGc                                                                                                                                            <  1:630166/150‑1 (MQ=255)           gCCCGCGCGTCCAGTGGCCGTACCCAACACTTCTTCCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCAGGCCGGTgaag                                                                                                                                            <  1:1270760/150‑4 (MQ=255)           gCCCGCGCGTCCAGTGGCCGTACCAAACACTTCTTCCTGCCCGACCCGTATAAACCCTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCAGGCCGGTgaag                                                                                                                                            <  1:1795675/150‑4 (MQ=255)            cccGCGCGTCCAGTGGCCGTACCCAACACTTCTTCCTGCCCGACCCGGCTAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCg                                                                                                                                           <  1:1144147/150‑1 (MQ=255)                cgcgTCCAGTGGCCGTACCCAACACTTCTTCCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGt                                                                                                                                       <  1:1030234/150‑1 (MQ=255)                 gcgTCCAGTGGCCGTACCCAACACTTCTTCCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTg                                                                                                                                      <  1:225525/150‑1 (MQ=255)                  cgTCCAGTGGCCGTACCCAACACTTCTTCCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTgg                                                                                                                                     <  1:965254/150‑1 (MQ=255)                     ccAGTGGCCGTACCCAACACTTCTTCCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGaaa                                                                                                                                  >  2:1072188/1‑150 (MQ=255)                      cAGTGGCCGTACCCAACACTTCTTCCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGaaaa                                                                                                                                 <  2:2369124/150‑1 (MQ=255)                      cAGTGGCCGTACCCAACACTTCTTCCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCAACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGaaaa                                                                                                                                 <  2:361658/150‑1 (MQ=255)                      cAGTGGCCGTACCCAACACTTCTTCCTGCCCCACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTCGTGGaaaa                                                                                                                                 <  2:837622/150‑1 (MQ=255)                       aGTGGCCGTACCCAACACTTCTTCCTGCCCGACCCGGATAAAACCTTCCACCGCGGCCGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAg                                                                                                                                >  1:272761/1‑150 (MQ=255)                          ggCCGTACCCAACACTTCTTCCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTCCTTAGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTc                                                                                                                             <  1:1956876/150‑1 (MQ=255)                           gCCGTACCCAACCCTTCTTCCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCg                                                                                                                            <  1:2153284/150‑1 (MQ=255)                           gCCGTACCCAACACTTCTTCCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCTATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCg                                                                                                                            <  1:1818679/150‑1 (MQ=255)                           gCCGTACCCAACACTTCTTCCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCg                                                                                                                            >  2:1114823/1‑150 (MQ=255)                               tatccaccaCTCCTTCCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCCCCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCaa                                                                                                                        <  1:866323/147‑1 (MQ=255)                                  ccAACACTTCTTCCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCg                                                                                                                     >  1:103827/1‑150 (MQ=255)                                  ccAACACTTCTTCCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCg                                                                                                                     >  2:905655/1‑150 (MQ=255)                                   cAACACTTCTTCCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCCCCGGCCTGCCGGTGGGCGTAGGGGAAAAGGGCAGCAcgcgc                                                                                                                    >  2:1338324/1‑149 (MQ=255)                                    aaCACTTCTTCCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTg                                                                                                                   >  2:2400692/1‑150 (MQ=255)                                        cttcttCCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACAt                                                                                                               >  1:2393116/1‑150 (MQ=255)                                         ttcttcCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGCCTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATg                                                                                                              >  2:1774601/1‑150 (MQ=255)                                          tcttcCTGCCCGACCCGGATAAAACCTTCCACCGCGGAAGCACCGACTCATTCGTCACCGCCCGCAAGGTCGATATCGGCGCTTTCGACTCTCCAACCTGCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGc                                                                                                             <  2:880449/150‑1 (MQ=255)                                            ttcCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCaa                                                                                                           <  2:505957/150‑1 (MQ=255)                                               cTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTg                                                                                                        <  2:1574796/150‑1 (MQ=255)                                                tGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTgg                                                                                                       >  1:1663570/1‑150 (MQ=255)                                                tGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTgg                                                                                                       <  1:1072188/150‑1 (MQ=255)                                                tGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTgg                                                                                                       <  2:1127851/150‑1 (MQ=255)                                                 gCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGt                                                                                                                                                <  1:1548588/108‑1 (MQ=255)                                                 gCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGt                                                                                                                                                >  2:1548588/1‑108 (MQ=255)                                                     gACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCAcc                                                                                                  >  1:1842124/1‑150 (MQ=255)                                                      aCCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCg                                                                                                 >  1:2288585/1‑150 (MQ=255)                                                        ccGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGAc                                                                                               >  2:554853/1‑150 (MQ=255)                                                        ccGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGAc                                                                                               <  1:2400692/150‑1 (MQ=255)                                                        ccGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAAACTTCACCGGCCTGCCGGTGGGCGTACTGGGATATGTCTGCAAGCGTGACATGCAAGTGGTCACCGAc                                                                                               >  2:1172176/1‑150 (MQ=255)                                                          ggATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGt                                                                                             >  2:1895740/1‑150 (MQ=255)                                                                  ccTTCCCGCGCAGCTACACCGACTACTTCGGGACCGCCCGCAAGGGCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGCAAAGGTGGGCAAGTGTTAGATGCAAGTGGTCACCGACGTACCGCTGa                                                                                     <  1:1318886/150‑1 (MQ=255)                                                                  ccTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGa                                                                                     >  1:1572398/1‑150 (MQ=255)                                                                  ccTTCCACCGCGGCAGCACCGACTACCTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGa                                                                                     >  2:54063/1‑150 (MQ=255)                                                                  ccTTCAACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCGGCACCGGCCTGCCGGTGGGCGGAGTGGAAAAGGTCGGCAAGCCTGACATGCAAGTGGTCACCGACGTACCGCTGa                                                                                     <  2:753941/150‑1 (MQ=255)                                                                   cTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACa                                                                                                                <  1:282417/122‑1 (MQ=255)                                                                   cTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACa                                                                                                                >  2:282417/1‑122 (MQ=255)                                                                   cTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGAc                                                                                    >  1:881710/1‑150 (MQ=255)                                                                        aCCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACg                                                                               <  1:336468/150‑1 (MQ=255)                                                                             ggCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAAACTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGAc                                                                          <  1:1172176/150‑1 (MQ=255)                                                                                  caCCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCg                                                                                         >  1:1026293/1‑130 (MQ=255)                                                                                  caCCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCg                                                                                         <  2:1026293/130‑1 (MQ=255)                                                                                     cGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCCAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGACGGCCTCaa                                                                  <  2:1404057/150‑1 (MQ=255)                                                                                     cGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGACGGCCTCaa                                                                  <  1:1114823/150‑1 (MQ=255)                                                                                       actactTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGACGGCCTCAAtc                                                                >  2:739490/1‑149 (MQ=255)                                                                                         tactTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCg                                                                                                                                                  <  1:126302/66‑1 (MQ=255)                                                                                         tactTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCg                                                                                                                                                  >  2:126302/1‑66 (MQ=255)                                                                                            ttCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGACGGCCTCAATGTGCTg                                                           <  2:1663570/150‑1 (MQ=255)                                                                                            ttCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGACGGCCTCAATGTGCTg                                                           <  1:725257/150‑1 (MQ=255)                                                                                            ttCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGACGGCCTCAATGTGCTg                                                           <  2:2288585/150‑1 (MQ=255)                                                                                            ctCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTTACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGACGGCCTCAATGTGCTg                                                           <  1:54063/149‑1 (MQ=255)                                                                                             tCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGACGGCCTCAATGTGCTgg                                                          >  2:1578347/1‑150 (MQ=255)                                                                                                  accgaccgCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGACGGCCTCAATGTGCTGGTCAAg                                                     >  2:933268/1‑150 (MQ=255)                                                                                                  accgaccgCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGACGGCCTCAATGTGCTGGTCAAg                                                     >  2:2348378/1‑150 (MQ=255)                                                                                                     gaccgCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGTCCAACGGCGACGGCCTCAATGTGCTGGTCAAgccc                                                  >  2:868614/1‑148 (MQ=255)                                                                                                        cgCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCCTCCCAGGCCTGCCGGTGGCCGTAGTGGAAAATGTCGGCAAGCGTGACATGCAAGTTGTCACCGACGTACCGCTGACCAACGGCGACGGCCTCAATGTGCTGGTCAAGCGCGaa                                               <  2:272761/150‑1 (MQ=255)                                                                                                         gCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAATCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGACGGCCTCAATGTGCTGGTCAAGCGCGAAg                                              >  1:2343390/1‑150 (MQ=255)                                                                                                          cAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGACGGCCTCAATGTGCTGGTCAAGCGCGAAgt                                             >  1:2239296/1‑150 (MQ=255)                                                                                                             ggTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGACGGCCTCAATGTGCTGGTCAAGCGCGAAgtggt                                          <  1:554853/150‑1 (MQ=255)                                                                                                              gTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGACGGCCTCAATGTGCTGGTCAAGCGCGAAgtggtg                                         <  1:652267/150‑1 (MQ=255)                                                                                                                 gATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGAc                                                                          <  1:2144287/114‑1 (MQ=255)                                                                                                                 gATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGAc                                                                          >  2:2144287/1‑114 (MQ=255)                                                                                                                    atCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGACGGCCTCAATGTGCTGGTCAAGCGCGAAGTGGTGGGCTTc                                   >  1:991025/1‑150 (MQ=255)                                                                                                                       ggCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGACGGCCTCAATGTGCTGGTCAAGCGCGAAGTGGTGGGCTTCcgc                                >  1:116026/1‑150 (MQ=255)                                                                                                                               cGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGACGGCCTCAATGTGCTGGTCAAGCGCGAAGTGGTGGGCTTCCGCGCCAACAt                        <  1:1895740/150‑1 (MQ=255)                                                                                                                                gACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGACGGCCTCAATGTGCTGGTCAAGCGCGAAGTGGTGGGCTTCCGCGCCAACATc                       >  2:1508841/1‑150 (MQ=255)                                                                                                                                   tCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGACGGCCTCAATGTGCTGGTCAAGCGCGAAGTGGTGGGCTTCCGCGCCAACATCGcc                    >  1:382237/1‑150 (MQ=255)                                                                                                                                     gCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGACGGCCTCAATGTGCTGGTCAAGCGCGAAGTGGTGGGCTTCCGCGCCAACATCGCCGa                  >  1:1665443/1‑150 (MQ=255)                                                                                                                                      ccAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGACGGCCTCAATGTGCTGGTCAAGCGCGAAGtggtgg                                        >  2:1322363/1‑127 (MQ=255)                                                                                                                                      ccAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGACGGCCTCAATGTGCTGGTCAAGCGCGAAGtggtgg                                        <  1:1322363/127‑1 (MQ=255)                                                                                                                                       cAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGACGGCCTCAATGTGCTGGTCAAGCGCGAAGTGGTGGGCTTCCGCGCCAACATCGCCGAAc                <  2:103827/150‑1 (MQ=255)                                                                                                                                           gcccACCGGCAGGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGACGGCCTCAATGTGCTGGTCAAGCGCGAAGTGGTGGGCTTCCGCGCCaa                           <  2:1762792/132‑1 (MQ=255)                                                                                                                                           gcccACCGGCAGGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGACGGCCTCAATGTGCTGGTCAAGCGCGAAGTGGTGGGCTTCCGCGCCaa                           >  1:1762792/4‑135 (MQ=255)                                                                                                                                           cTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGACGGCCTCAATGTGCTGGTCAAGCGCGAAgtggtg                                         <  1:636435/121‑1 (MQ=255)                                                                                                                                           cTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGACGGCCTCAATGTGCTGGTCAAGCGCGAAGTGGTGGGCTTCCGCGCCAACATCGCCGAACcgcg            >  1:1587716/1‑150 (MQ=255)                                                                                                                                           cTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACAGGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGACGGCCTCAATGTGCTGGTCAAGCGCGAAgtggtg                                         >  2:636435/1‑121 (MQ=255)                                                                                                                                            ttCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGACGGCCTCAATGTGCTGGTCAAGCGCGAAGTGGTGGGATGGCTCGCCAACATCGCCGAACCGCGt           >  1:611192/1‑150 (MQ=255)                                                                                                                                            cccACCGGCAGGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGACGGCCTCAATGTGCTGGTCAAGCGCGAAGTGGTGGGCTTCCGCGCCAACATCGcc                    <  2:1824282/139‑1 (MQ=255)                                                                                                                                            cccACCGGCAGGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGACGGCCTCAATGTGCTGGTCAAGCGCGAAGTGGTGGGCTTCCGCGCCAACATCGcc                    >  1:1824282/3‑141 (MQ=255)                                                                                                                                            cccACCGGCAGGCCGGTGGGCGAAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGACGGCCTCAATGTGCTGGTCAAGCGCGAAGTGGTGGGCTTCCGCGCCAACATCGCCGAACCGCGt           >  2:807505/3‑150 (MQ=255)                                                                                                                                               aCCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGACGGCCTCAATGTGCTGGTCAAGCGCGAAGTGGTGGGCTTCCGCGCCAACATCGCCGAACCGCGTGGc        <  2:1587716/150‑1 (MQ=255)                                                                                                                                                ccGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGACGGCCTCAATGTGCTGGTCAAGCGCGAAGTGGTGGGCTTCCgcgc                              >  1:139876/1‑127 (MQ=255)                                                                                                                                                ccGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGACGGCCTCAATGTGCTGGTCAAGCGCGAAGTGGTGGGCTTCCgcgc                              <  2:139876/127‑1 (MQ=255)                                                                                                                                                 cGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGACGGCCTCAATGTGCTGGTCAAGCGCGAAGTGGTGGGCTTCCGCGCCAACATCGCCGAACCGCGTGGCGa      >  2:984228/1‑150 (MQ=255)                                                                                                                                                  ggCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGGCATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGACGGCCTCAATGTGCTGGTCAAGCGCGAAGTGGTGGGCTTCCGCGCCAACATCGCCGAACCGCGTGGCGaa     >  2:775495/1‑150 (MQ=255)                                                                                                                                                     cTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGACGGCCTCAATGTGCTGGTCAAGCGCGAAGTGGTGGGCTTCCGCGCCAACATCGCCGAACCGCGTGGCGAATTc  <  2:116026/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                      CCCTCAGTACGCCCGCGCGTCCAGTGGCCGTACCCAACACTTCTTCCTGCCCGACCCGGATAAAACCTTCCACCGCGGCAGCACCGACTACTTCGTGACCGACCGCAAGGTCGATATCGGCGCCTTCGACTCGCCAACCTTCACCGGCCTGCCGGTGGGCGTAGTGGAAAAGGTCGGCAAGCGTGACATGCAAGTGGTCACCGACGTACCGCTGACCAACGGCGACGGCCTCAATGTGCTGGTCAAGCGCGAAGTGGTGGGCTTCCGCGCCAACATCGCCGAACCGCGTGGCGAATTC  >  NC_002947/3966658‑3966955

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑