Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,332	C→A	7.1%	R17R (CGG→CGT)	trmE ←	GTPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	6,332	0	C	A	7.1%	313.8 / 11.0	99	R17R (CGG→CGT)	trmE	GTPase
Reads supporting (aligned to +/- strand):  ref base C (49/42);  new base A (3/4);  total (52/47)
Fisher's exact test for biased strand distribution p-value = 7.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.46e-01

AGCAATCCCTTCGTCCAGCACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCT‑GGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAG  >  NC_002947/6188‑6478                                                                                                                                                 |                                                                                                                                                     agCAATCCCTTCGTCCAGCACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGc                                                                                                                                                  >  4:7747/1‑149 (MQ=255)      tCCCTTCGTCCAGCACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCT‑gg                                                                                                                                            <  3:7747/149‑1 (MQ=255)             gTCCAGCACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGGGCGGCCGGTGATCAACCGGCCGGCCTTGGCAGCCAGACGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCT‑GGGCCGTgg                                                                                                                                     <  7:227279/149‑1 (MQ=255)             gTCCAGCACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGGCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCT‑GGGCCGTg                                                                                                                                      <  7:67333/148‑1 (MQ=255)               ccAGCACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGACCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCT‑GGGCCGTGGc                                                                                                                                    >  8:49598/1‑148 (MQ=255)                cAGCACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCT‑GGGCCGTGGCg                                                                                                                                   <  2:30873/148‑1 (MQ=255)                cAGCACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCT‑GGGCCGTGGCGa                                                                                                                                  >  7:159713/1‑149 (MQ=255)                cAGCACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCT‑GGGCCGTGGCGa                                                                                                                                  <  7:33387/149‑1 (MQ=255)                  gCACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCATGGTGCGGCCGGTGATCAACAGGGCGGCCTTGGCAGCCAGAGGGCCGTACATCCTGACAATTCCCACACCGCCCCGGCCCT‑GGGCCGTGGCGATg                                                                                                                                <  1:278521/149‑1 (MQ=255)                    aCCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCT‑GGGCCGTGGCGATggc                                                                                                                              >  3:325303/1‑149 (MQ=255)                    aCCAGCCCTTCATCGGCGCGGAATCGACCTTACTGGGCATGCCGCGGGGTCTGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCACGGCCCA‑GGGCCGGGGCGGTgtg                                                                                                                              <  8:277901/149‑3 (MQ=255)                     ccAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCT‑GGGCCGTGGCGATggcg                                                                                                                             >  6:281070/1‑149 (MQ=255)                      cAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCT‑GGGCCGTGGCGATggcgg                                                                                                                            <  1:26769/149‑1 (MQ=255)                         cccTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGGCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCT‑GGGCCGTGGCGATGGCGGCGa                                                                                                                         <  6:341582/149‑1 (MQ=255)                           cTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACATCGCCACGGCCCA‑GGGCCGGGGCGGTgtgg                                                                                                                             <  5:85579/143‑4 (MQ=255)                           cTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACATCGCCACGGCCCA‑GGGCCGGGGCGGTgtgg                                                                                                                             >  6:85579/1‑140 (MQ=255)                           cTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCT‑GGGCCGTGGCGATGGCGGCGATg                                                                                                                       >  4:65266/1‑149 (MQ=255)                            ttCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCT‑GGGCCGTGGCGATGGCGGCGATgg                                                                                                                      >  7:138366/1‑149 (MQ=255)                              cATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCT‑GGGCCGTGGCGATGGCGGCGATgg                                                                                                                      <  4:93739/147‑1 (MQ=255)                                   tcgCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCT‑GGGCCGTGGCGATGGCGGCGATGGTTTCACg                                                                                                               >  8:170364/1‑149 (MQ=255)                                     gcgGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCACGGCCCA‑GGGCCGGGGCGGTgtg                                                                                                                              <  6:109539/132‑3 (MQ=255)                                     gcgGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCACGGCCCA‑GGGCCGGGGCGGTgtg                                                                                                                              >  5:109539/1‑130 (MQ=255)                                       ggAATGGACCGTAATGGGCATGCGGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCT‑GGGCCGTGGCGATGGCGGCGATGGTTTCACGcaca                                                                                                           <  6:36022/149‑1 (MQ=255)                                         aaTGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCT‑GGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAgt                                                                                                         >  1:225557/1‑149 (MQ=255)                                              aCCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCT‑GGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCa                                                                                                    <  7:170364/149‑1 (MQ=255)                                              aCCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCT‑GGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCa                                                                                                    <  3:215499/149‑1 (MQ=255)                                                cGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCT‑GGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATg                                                                                                  >  5:185134/1‑149 (MQ=255)                                                cGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCA‑GGGCCGGGGCGGTgtg                                                                                                                              >  8:160067/1‑119 (MQ=255)                                                 gTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCT‑GGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGc                                                                                                 <  8:159713/149‑1 (MQ=255)                                                  tAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCT‑GGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCt                                                                                                >  3:101658/1‑149 (MQ=255)                                                   aaTGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACATGACAATTCCCACACCGCCCCGGCCCT‑GGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTc                                                                                               >  8:228431/1‑149 (MQ=255)                                                             ccGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCT‑GGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCctc                                                                                      >  8:315077/1‑148 (MQ=255)                                                             ccGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCT‑GGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCctc                                                                                      <  2:40300/148‑1 (MQ=255)                                                               gcgGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCT‑GGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTAc                                                                                   <  5:281070/149‑1 (MQ=255)                                                               gcgGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCT‑GGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTAc                                                                                   <  8:45947/149‑1 (MQ=255)                                                               gcgGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCT‑GGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTAc                                                                                   <  3:21135/149‑1 (MQ=255)                                                                  gggTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCT‑GGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGAc                                                                                <  3:65266/149‑1 (MQ=255)                                                                       aGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCT‑GGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACaaaaa                                                                           <  7:228431/149‑1 (MQ=255)                                                                       aGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCT‑GGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACaaaaa                                                                           >  8:219990/1‑149 (MQ=255)                                                                       aGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCT‑GGGCCGTGGCGATGGCGGAGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACaaaaa                                                                           >  2:283812/1‑149 (MQ=255)                                                                         ggTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCT‑GGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAgtgt                                                                                                       >  3:212246/1‑119 (MQ=255)                                                                         ggTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCT‑GGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAgtgt                                                                                                       <  4:212246/119‑1 (MQ=255)                                                                          gTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCT‑GGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGa                                                                        <  7:219990/149‑1 (MQ=255)                                                                           tGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCT‑GGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACTAAAACGAc                                                                       <  4:39275/149‑1 (MQ=255)                                                                                ccGGTGATCAACAGGCCGACCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCT‑GGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATa                                                                  <  7:49598/149‑1 (MQ=255)                                                                                   gTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCT‑GGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGc                                                                >  4:7823/1‑148 (MQ=255)                                                                                    tGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCT‑GGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCaa                                                              <  4:101658/149‑1 (MQ=255)                                                                                     gATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCT‑GGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCaa                                                              >  3:181262/1‑148 (MQ=255)                                                                                       tCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCT‑GGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCaaaaa                                                           <  4:325303/149‑1 (MQ=255)                                                                                       tCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCT‑GGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCaaaaa                                                           >  5:182017/1‑149 (MQ=255)                                                                                        cAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCT‑GGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAAc                                                          >  1:216668/1‑149 (MQ=255)                                                                                           cAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCT‑GGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGcc                                                       <  2:114555/149‑1 (MQ=255)                                                                                              gccggccTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCT‑GGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCa                                                    >  1:243904/1‑149 (MQ=255)                                                                                                cggccTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCT‑GGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCctc                                                                                      <  4:288868/113‑1 (MQ=255)                                                                                                  gccTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCT‑GGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAg                                                >  7:244993/1‑149 (MQ=255)                                                                                                       ggCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCT‑GGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAg‑ggg                                            <  2:225557/147‑1 (MQ=255)                                                                                                       ggCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCT‑GGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGc                                          >  8:120455/1‑149 (MQ=255)                                                                                                       ggCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCT‑GGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCC‑CTAg‑tggg                                           <  8:138366/147‑1 (MQ=255)                                                                                                        gCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCT‑GGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGa                                                                                 >  1:139010/1‑110 (MQ=255)                                                                                                        gCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCT‑GGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGa                                                                                 <  2:139010/110‑1 (MQ=255)                                                                                                        gCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCT‑GGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAACCGACAGATAGCAAAAAAGCCCCACTAG‑GGGGCg                                         >  2:238643/1‑149 (MQ=255)                                                                                                        gCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCT‑GGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCg                                         >  6:210340/1‑149 (MQ=255)                                                                                                            ccAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCT‑GGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGtttt                                     >  3:124843/1‑149 (MQ=255)                                                                                                             cAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCT‑GGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGttttt                                    >  8:265029/1‑149 (MQ=255)                                                                                                               gagGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCT‑GGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGTTTTTTa                                  >  7:261132/1‑149 (MQ=255)                                                                                                               gagGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCT‑GGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑CGGGCGTTTTTTa                                  <  2:216668/149‑1 (MQ=255)                                                                                                                  ggCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCT‑GGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAATAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGTTTTTTATTc                               <  1:238643/149‑1 (MQ=255)                                                                                                                     cGGACAACCTGACAATTCCCACACCGCCCCGGCCCT‑GGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCC‑CTAGTGGGGCGttttt                                    >  4:213091/1‑141 (MQ=255)                                                                                                                     cGGACAACCTGACAATTCCCACACCGCCCCGGCCCT‑GGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCC‑CTAGTGGGGCGttttt                                    <  3:213091/141‑1 (MQ=255)                                                                                                                         cAACCTGACAATTCCCACACCGCCCCGGCCCT‑CGGCCGAGGCGATGGCGGCGATGGTTTCACGCACAGTGTTAATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGAGTTTATTATTAACAGGc                         >  3:67362/1‑148 (MQ=255)                                                                                                                            ccTGACAATTCCCACACCGCCCCGGCCCTAGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGTTTTTTATTCACAGGCTAg                      >  1:198424/1‑149 (MQ=255)                                                                                                                            ccTGACAATTCCCACACCGCCCCGGCCCT‑GGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGTTTTTTATTCACAGGCTAGc                     <  7:120455/149‑1 (MQ=255)                                                                                                                            ccTGACAATTCCCACACCGCCCCGGCCCT‑GGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGTTTTTTATTCACAGGCTAGc                     >  3:89307/1‑149 (MQ=255)                                                                                                                              tGACAATTCCCACACCGCCCCGGCCCT‑GGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGTTTTTTATTCACAGGCTAGCg                    >  3:159014/1‑148 (MQ=255)                                                                                                                              tGACAATTCCCACACCGCCCCGGCCCT‑GGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGATCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGTTTTTTATTCACAGGCTAGCg                    >  8:274597/1‑148 (MQ=255)                                                                                                                               gACAATTCCCACACCGCCCCGGCCCT‑GGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGTTTTTTATTCACAGGCTAGCg                    >  7:188541/1‑147 (MQ=255)                                                                                                                               gACAATTCCCACACCGCCCCGGCCCT‑GGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGTTTTTTATTCACAGGCTAGCGTa                  >  2:337948/1‑149 (MQ=255)                                                                                                                               gACAATTCCCACACCGCCCCGGCCCT‑GGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCC‑CTAGTGGGGCGtttttgcta                                >  1:315053/1‑131 (MQ=255)                                                                                                                               gACAATTCCCACACCGCCCCGGCCCT‑GGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCC‑CTAGTGGGGCGtttttgcta                                <  2:315053/135‑5 (MQ=255)                                                                                                                                  aaTTCCCACACCGCCCCGGCCCT‑GGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGGTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGTTTTTTATTCACAGGCTAGCGTAt                 >  8:252764/1‑147 (MQ=255)                                                                                                                                    ttCCCACACCGCCCCGGCCCT‑GGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGTTTTTTATTCACAGGCTAGCGTATTAGc             >  2:253554/1‑149 (MQ=255)                                                                                                                                     tCCCACACCGCCCCGGCCCT‑GGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAACGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGTTTTTTATTCACAGGCTAGCGTATTAGcc            >  5:237343/1‑149 (MQ=255)                                                                                                                                       ccACTCCGCCGCGGCGCT‑GGGCCGTGGCGATGGCGGCGTTGGTTTCACGCACTGTGTTCATGCTCGAAAGCCTCTACGTCAAAATCGACAGTTAGCAAATACGCCCCACTTG‑GTGGCGTTTTTTATTGACTGGCTAGCGTa                  <  4:276222/141‑1 (MQ=255)                                                                                                                                       ccACACCGCCCCGGCCCT‑GGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGTTTTTTATTCACAGGCTAGCGTa                  >  3:276222/1‑141 (MQ=255)                                                                                                                                       ccACACCGCCCCGGCCCT‑GGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGTTTTTTATTCACAGGCTAGCGTATTAGccc           <  4:89307/148‑1 (MQ=255)                                                                                                                                       ccACACCGCCCCGGCCCT‑GGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGTTTTTTATTCACAGGCTAGCGTATTAGccc           <  2:198424/148‑1 (MQ=255)                                                                                                                                       ccACACCGCCCCGGCCCT‑GGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCg          >  7:23476/1‑149 (MQ=255)                                                                                                                                       ccACACCGCCCCGGCCCT‑GGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCg          >  2:219095/1‑149 (MQ=255)                                                                                                                                        cacaCCGCCCCGGCCCT‑GGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGt         >  5:302907/1‑149 (MQ=255)                                                                                                                                         ccatcGCCACGGCCCA‑GGGCCGGGGCGGTGTGGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTc        <  7:22846/145‑1 (MQ=255)                                                                                                                                           acCGCCCCGGCCCT‑GGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAGACGCCCCACTAG‑GGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTc        >  6:303548/1‑147 (MQ=255)                                                                                                                                           acCGCCCCGGCCCT‑GGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGGTTTTGATTCACAGGCTAGCGTATTAGCCCGTc        >  6:146713/1‑147 (MQ=255)                                                                                                                                            ccGCCCCGGCCCT‑GGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAag     <  4:124843/149‑1 (MQ=255)                                                                                                                                            ccGCCCCGGCCCT‑GGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAag     <  7:265029/149‑1 (MQ=255)                                                                                                                                             cGCCCCGGCCCT‑GGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAagc    <  8:244993/149‑1 (MQ=255)                                                                                                                                              gCCCCGGCCCT‑GGGCCGTGGCGATGGCTGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGTTTTTTATTGACAGGCTAGCGTATTAGCCCGTCAagca   <  4:67362/149‑1 (MQ=255)                                                                                                                                              gCCCCGGCCCT‑GGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGTTTTTTATTcaca                            <  4:238725/124‑1 (MQ=255)                                                                                                                                              gCCCCGGCCCT‑GGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGTTTTTTATTcaca                            >  3:238725/1‑124 (MQ=255)                                                                                                                                              gCCACGGCCCA‑GGGCCGGGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAag     >  3:226433/1‑147 (MQ=255)                                                                                                                                               ccccGGCCCT‑GGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAagcag  >  1:329236/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                                     AGCAATCCCTTCGTCCAGCACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCT‑GGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAG  >  NC_002947/6188‑6478

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