Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	838,007:1	+C	5.8%	intergenic (‑43/+37)	pth ← / ← rplY	peptidyl‑tRNA hydrolase/50S ribosomal protein L25/general stress protein Ctc

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	838,007	1	.	C	5.8%	278.3 / 9.8	69	intergenic (‑43/+37)	pth/rplY	peptidyl‑tRNA hydrolase/50S ribosomal protein L25/general stress protein Ctc
Reads supporting (aligned to +/- strand):  ref base . (36/29);  new base C (2/2);  total (38/31)
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.87e-01

GCAAGGCGTTCAACGAAAAGAGCCCCTGCGTTATGCCGGGTCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTG‑CGCGCG‑GCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCTTTGTCATCACCGTGGGCCAGAGCGACGAACTCTACGCCTTTCGGA  >  NC_002947/837867‑838149                                                                                                                                               |                                                                                                                                               gCAAGGCGTTCAACGAAAAGAGCCCCTGCGTTATGCCGGGTCTGTTCGTTTTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTG‑CGCGCG‑GCACTGt                                                                                                                                         <  1:269055/148‑1 (MQ=255)  cAAGGCGTTCAACGAAAAGAGCCCCTGCGTTATGCCGGGTCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTG‑CGCGCG‑GCACTGTaa                                                                                                                                       <  6:86197/149‑1 (MQ=255)    aGGCGTTCAACGAAAAGAGCCCCTGCGTTATGCCGGGTCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTG‑CGCGCG‑GCACTGTAAcc                                                                                                                                     <  4:282175/149‑1 (MQ=255)             aCGAAAAGAGCCCCTGCGTTATGCCGGGTCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTG‑CGCGCG‑GCACTGTAACCCGATCAACg                                                                                                                            >  5:166387/1‑149 (MQ=255)             aCGAAAAGAGCCCCTGCGTTATGCCGGGTCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTG‑C‑CGCGCGCACTGTAAccca                                                                                                                                   <  2:84621/142‑2 (MQ=255)             aCGAAAAGAGCCCCTGCGTTATGCCGGGTCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTG‑C‑CGCGCGCACTGTAAccca                                                                                                                                   >  1:84621/1‑141 (MQ=255)               gAAAAGAGCCCCTGCGTTATGCCGGGTCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTG‑CGCGCG‑GCACTGTAACCCGATCAACGcc                                                                                                                          >  3:26294/1‑149 (MQ=255)                 aaaGAGCCCCTGCGTTATGCCGGGTCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTG‑CGCGCG‑GCACTGTAACCCGATCAACGCCgg                                                                                                                        >  2:23427/1‑149 (MQ=255)                       ccccTGCGTTATGCCGGGTCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTG‑CGCGCG‑GCACTGTAACCCGATCAACGCCGGCGTGc                                                                                                                   >  5:224638/1‑148 (MQ=255)                          cTGCGTTATGCCGGGTCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTG‑CGCGCG‑GCACTGTAACCCGATCAACGCCGGCGTGCGAGt                                                                                                               >  6:182074/1‑149 (MQ=255)                          cTGCGTTATGCCGGGTCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTG‑CGCGCG‑GCACTGTAACCCGATCAACGCCGGCGTGCGAGt                                                                                                               >  7:49266/1‑149 (MQ=255)                          cTGCGTTATGCCGGGTCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTG‑CGCGCG‑GCACTGTAACCCGATCAACGCCGGCGTGCGAGt                                                                                                               <  4:26294/149‑1 (MQ=255)                          cTGCGTTATGCCGGGTCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTG‑CGCGCG‑GCACTGTAACCCGATCAACGCCGGCGTGCGAGt                                                                                                               >  4:166667/1‑149 (MQ=255)                             cGTTATGCCGGGTCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTG‑CGCGCG‑GCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGa                                                                                                            <  7:157353/149‑1 (MQ=255)                                         tCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTG‑CGCGCG‑GCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAgc                                                                                                <  8:49266/149‑1 (MQ=255)                                             ttCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTG‑CGCGCG‑GCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCct                                                                                            >  4:121406/1‑149 (MQ=255)                                             ttCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTG‑CGCGCG‑GCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCct                                                                                            >  8:99215/1‑149 (MQ=255)                                             ttCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTG‑CGCGCG‑GCACTGTAACACGATCAACGACGGCGTGCGAGTGGATTACTCGGCAGCGCct                                                                                            >  7:206772/1‑149 (MQ=255)                                               cGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTG‑CGCGCG‑GCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCctt                                                                                           <  3:121406/148‑1 (MQ=255)                                                   ttCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTG‑CGCGCG‑GCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCcttctgct                                                                                      >  3:23640/1‑149 (MQ=255)                                                   tcgttggCCTGGGTAACCCCGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTG‑CGCGCG‑GCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCcttctgct                                                                                      >  3:269814/6‑149 (MQ=255)                                                     cGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTG‑CGCGCG‑GCACTGTAACCCGATCAACGCCGGCGTGCTAGTGGATTACTCGGCAGCGCcttctgctt                                                                                     <  1:23427/148‑1 (MQ=255)                                                      ggggCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTG‑CGCGCG‑GCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCcttc                                                                                          >  5:38602/1‑142 (MQ=255)                                                      ggggCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTG‑CGCGCG‑GCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCcttc                                                                                          <  6:38602/142‑1 (MQ=255)                                                       gggCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTG‑CGCGCG‑GCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCcttctgcttctg                                                                                  <  6:166387/149‑1 (MQ=255)                                                        ggCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTG‑CGCGCG‑GCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTgg                                                                                 <  6:224638/149‑1 (MQ=255)                                                              ggTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTG‑CGCGCG‑GCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAAc                                                                           >  7:140102/1‑149 (MQ=255)                                                                   cccAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTG‑CGCGCG‑GCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTg                                                                      <  8:140102/149‑1 (MQ=255)                                                                    ccAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTG‑CGCGCG‑GCACTGTAAc                                                                                                                                      >  5:243391/1‑84 (MQ=255)                                                                    ccAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTG‑CGCGCG‑GCACTGTAAc                                                                                                                                      <  6:243391/84‑1 (MQ=255)                                                                    ccAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTG‑CGCGCG‑GCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGAt                                                                                                           <  5:276701/111‑1 (MQ=255)                                                                    ccAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTG‑CGCGCG‑GCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGAt                                                                                                           >  6:276701/1‑111 (MQ=255)                                                                       ggCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTG‑CGCGCG‑GCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGc                                                                  >  5:306906/1‑149 (MQ=255)                                                                         ccAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTG‑CGCGCG‑GCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGt                                                                <  6:306906/149‑1 (MQ=255)                                                                          cAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTG‑CGCGCG‑GCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTc                                                                                                      <  8:203165/110‑1 (MQ=255)                                                                          cAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTG‑CGCGCG‑GCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTc                                                                                                      >  7:203165/1‑110 (MQ=255)                                                                          cAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTG‑CGCGCG‑GCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGGGCAACGCGTgg                                                                     >  6:123218/1‑143 (MQ=255)                                                                          cAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTG‑CGCGCG‑GCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGGGCAACGCGTgg                                                                     <  5:123218/143‑1 (MQ=255)                                                                           aaCGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTG‑CGCGCG‑GCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTAGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTg                                                               >  8:222662/1‑148 (MQ=255)                                                                                   aCTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTG‑CGCGCG‑GCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTgg                                                       >  3:101153/1‑148 (MQ=255)                                                                                   aCTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTG‑CGCGCG‑GCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTgg                                                       >  3:67028/1‑148 (MQ=255)                                                                                       gATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTG‑CGCGCG‑GCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTgg                                                       >  3:189470/1‑144 (MQ=255)                                                                                       gATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTG‑CGCGCG‑GCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTgg                                                       <  4:189470/144‑1 (MQ=255)                                                                                            cGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTG‑CGCGCG‑GCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCttt                                             <  4:269814/149‑1 (MQ=255)                                                                                            cGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTG‑CGCGCG‑GCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCttt                                             >  2:235221/1‑149 (MQ=255)                                                                                            cGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTG‑CGCGCG‑GCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCttt                                             <  4:23640/149‑1 (MQ=255)                                                                                                cACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTG‑CGCGCG‑GCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCTTTGtca                                         <  1:235221/149‑1 (MQ=255)                                                                                                  cGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTG‑CGCGCG‑GCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCcttctgctt                                                                                     <  7:274299/103‑1 (MQ=255)                                                                                                  cGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTG‑CGCGCG‑GCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCcttctgctt                                                                                     >  8:274299/1‑103 (MQ=255)                                                                                                   gTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTG‑CGCGCG‑GCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCTTTGtcatca                                      <  3:166667/149‑1 (MQ=255)                                                                                                    tCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTG‑CGCGCG‑GCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGGGCAACGCGTGGAGCGTGAACGTTGGCAACAGCTTTGTCATCAc                                     >  4:309045/1‑149 (MQ=255)                                                                                                         ggCTCTTCCTTGGAGTTGGTGGGTTACAGTG‑CGCGCG‑GCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCTTTGTCATCACCGTgg                                <  6:232041/149‑1 (MQ=255)                                                                                                              ttCCTTGGAGTTGGTGGGTTACAGTG‑CGCGCG‑GCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTgg                                                       <  2:92932/121‑1 (MQ=255)                                                                                                              ttCCTTGGAGTTGGTGGGTTACAGTG‑CGCGCG‑GCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTgg                                                       >  1:92932/1‑121 (MQ=255)                                                                                                                ccTTGGAGTTGGTGGGTTACAGTG‑CGCGCG‑GCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCTTTGTCATCACCGTGGGCCagag                         >  2:258681/1‑149 (MQ=255)                                                                                                                ccTTGGAGTTGGTGGGTTACAGTG‑CGCGCG‑GCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCTTTGTCATCACCGTGGGCCaga                          <  7:99215/148‑1 (MQ=255)                                                                                                                 cTTGGAGTTGGTGGGTTACAGTG‑CGCGCG‑GCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCTTTGTCATCACCGTGGGCCAGAGc                        >  4:260071/1‑149 (MQ=255)                                                                                                                 cTTGGAGTTGGTGGGTTACAGTG‑CGCGCG‑GCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCTATGTCATCACCGTGGGCCagacc                        >  5:163403/1‑147 (MQ=255)                                                                                                                     gAGTTGGTGGGTTACAGTG‑CGCGCG‑GCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCTTTGTCATCACCGTGGGCCAGAGcgac                     <  1:258681/148‑1 (MQ=255)                                                                                                                     gAGTTGGTGGGTTACAGTG‑CGCGCG‑GCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCTTTGTCATCACCGTGGGCCAGAGcgac                     >  3:270634/1‑148 (MQ=255)                                                                                                                      aGTTGGTGGGTTACAGTG‑CGCGCG‑GCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGtt                                                         <  1:271829/111‑1 (MQ=255)                                                                                                                      aGTTGGTGGGTTACAGTG‑CGCGCG‑GCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGtt                                                         >  2:271829/1‑111 (MQ=255)                                                                                                                      aGTTGGTGGGTTACAGTG‑CGCGCG‑GCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCTTTGtca                                         <  7:98900/127‑1 (MQ=255)                                                                                                                      aGTTGGTGGGTTACAGTG‑CGCGCG‑GCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCTTTGtca                                         >  8:98900/1‑127 (MQ=255)                                                                                                                           gtggGTTACAGTG‑CGCGCG‑GCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCTTTGTCATCACCGTGGGCCAGAGCGACGAActc               >  5:132123/1‑148 (MQ=255)                                                                                                                                ttACAGTG‑C‑CGCGCGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCttt                                             >  3:153416/1‑113 (MQ=255)                                                                                                                                ttACAGTG‑C‑CGCGCGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCttt                                             <  4:153416/113‑1 (MQ=255)                                                                                                                                   cAGTGCCGCGC‑‑GCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCTTTGTCATCACCGTGGGCCAGAGCGACGAACTCTACGCCttt      >  3:3893/1‑149 (MQ=255)                                                                                                                                      tG‑CGCGCG‑GCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCTTTGTCATCACCGTGGGCCAGAGCGACGAACTCTACGCCTTTCg    >  1:154836/1‑148 (MQ=255)                                                                                                                                        ccgcgc‑‑gCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCTTTGTCATCACCGTGGGCCAGAGCGACGAACTCTACGCCTTTCGGa  <  4:3893/147‑1 (MQ=255)                                                                                                                                         cgcgcg‑GCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCTTTGTCATCACCGTGGGCCAGAGCGACGAACTCTACGCCTTTCGGa  <  3:260071/148‑1 (MQ=255)                                                                                                                                               |                                                                                                                                               GCAAGGCGTTCAACGAAAAGAGCCCCTGCGTTATGCCGGGTCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTG‑CGCGCG‑GCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCTTTGTCATCACCGTGGGCCAGAGCGACGAACTCTACGCCTTTCGGA  >  NC_002947/837867‑838149

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