Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	535,048	T→A	5.9%	intergenic (+33/‑87)	nusG → / → rplK	transcription antipausing factor NusG/50S ribosomal protein L11

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	535,048	0	T	A	5.9%	271.1 / 2.8	92	intergenic (+33/‑87)	nusG/rplK	transcription antipausing factor NusG/50S ribosomal protein L11
Reads supporting (aligned to +/- strand):  ref base T (37/48);  new base A (2/4);  total (39/53)
Fisher's exact test for biased strand distribution p-value = 6.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAATGGTAGTGTCGAAGAAGTTAACTACGAAAAGAGTCGCCTGCAGGTTGCGGTGCTCATTTTCGGTCGCTCTACTCCGGTGGAGCTCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGC‑ATCCCATACCCCGCAGCCCTGTGTGCT‑‑GCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCGGGGAGCCATCAGGCGTTCGTACCCGATTTTGGAGTAGCTCATGGCTAAGAAGATTCAGGCTTACATCAAGCTGCAAGTTAAGGCCGGCCAGGCCAACCC  >  NC_002947/534901‑535193                                                                                                                                                     |                                                                                                                                                    cAATGGTAGTGTCGAAGAAGTTAACTACGAAAAGAGTCGCCTGCAGGTTGCGGTGCTCATTTTCGGTCGCTCTACTCCGGTGGAGCTCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGC‑ATCCCATACCCCGCAGCCCtgtg                                                                                                                                                    >  8:326246/1‑149 (MQ=255) cAATGGTAGTGTCGAAGAAGTTAACTACGAAAAGAGTCGCCTGCAGGTTGCGGTGCTCATTTTCGGTCGCTCTACTCCGGTGGAGCTCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGC‑ATCCCATACCCCGCAGCCCtgt                                                                                                                                                     <  6:279891/148‑1 (MQ=255)  aaTGGTAGTGTCGAAGAAGTTAACTACGAAAAGAGTCGCCTGCAGGTTGCGGTGCTCATTTTCGGTCGCTCTACTCCGGTGGAGCTCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGC‑ATCCCATACCCCGCAGCCCtgtgt                                                                                                                                                   >  4:88785/1‑149 (MQ=255)         gtgtCGAAGAAGTTAACTACGAAAAGAGTCGCCTGCAGGTTGCGGTGCTCATTTTCGGTCGCTCTACTCCGGTGGAGCTCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGC‑ATCCCATACCCCGCAGCCCTGTGTGCT‑‑GCg                                                                                                                                           <  3:272714/148‑1 (MQ=255)          tgtCGAAGAAGTTAACTACGAAAAGAGTCGCCTGCAGGTTGCGGTGCTCATTTTCGGTCGCTCTACTCCGGTGGAGCTCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGC‑ATCCCATACCCCGCAGCCCTGTGTGCT‑‑GCggg                                                                                                                                         <  5:46981/149‑1 (MQ=255)           ttctaagaagTTAACTACGAAAAGATTCGCCTGCAGGTTGCGTTGCTCATTTTCGGTCTCTCTACTCCGGTTGAGTTTGAGTTTAGTCATTTCGAAAAGGTCTATGCGGACGAAGC‑ATCCCATACCCCGCAGCACACAGGGCT‑‑GCggg                                                                                                                                         <  7:56458/144‑1 (MQ=255)           gtCGAAGAAGTTAACTACGAAAAGAGTCGCCTGCAGGTTGCGGTGCTCATTTTCGGTCGCTCTACTCCGGTGGAGCTCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGC‑ATCCCATACCCCGCAGCCCTGTGTGCT‑‑GCggg                                                                                                                                         <  2:341585/148‑1 (MQ=255)           cccGAAGAAGCTACCTCCGCAACGAGGCGCCCGCCCGCTGCGGGGCTCATTTTCGGTCGCTCTACTCCGGTGGAGCTCGAGTTCAGCCAGGTCGAAAAGGTCTCGGCGGACCAAGC‑ACCCCAAACCCCGCACCACCCCGGGCT‑‑GCgggg                                                                                                                                        <  4:41344/147‑1 (MQ=255)            tCGAAGAAGTTAACTACGAAAAGAGTCGCCTGCAGGTTGCGGTGCTCATTTTCGGTCGCTCTACTCCGGTGGAGCTCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGC‑ATCCCATACCCCGCAGCCCTGTGTGCT‑‑GCgg                                                                                                                                          >  5:4814/1‑146 (MQ=255)            tCGAAGAAGTTAACTACGAAAAGAGTCGCCTGCAGGTTGCGGTGCTCATTTTCGGTCGCTCTACTCCGGTGGAGCTCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGC‑ATCCCATACCCCGCAGCCCTGTGTGCT‑‑GCgg                                                                                                                                          <  6:4814/146‑1 (MQ=255)            tCGAAGAAGTTAACTACGAAAAGAGTCGCCTGCAGGTTGCGGTGCTCATTTTCGGTCGCTCTACTCCGGTGGAGCTCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGC‑ATCCCATACCCCGCAGCACACAGGGCT‑‑GCgggg                                                                                                                                        <  4:21460/148‑1 (MQ=255)              gaagaaGTTAACTACGAAAAGAGTCGCCTGGAGGTTGCGGTGCTCATTTTCGGTCGCTCTACTCCGGTGGAGCTCGAGTTCAGCCATGTCGAAAAGGTCTAGGCGGACGAAGC‑ATCCCATACCCCGCAGCCCTGTGTGCT‑‑GCGGGGtt                                                                                                                                      <  8:84568/148‑1 (MQ=255)              gaagaaGTTAACTACGAAAAGAGTCGCCTGCAGGTTGCGGTGCTCATTTTCGGTCGCTCTACTCCGGTGGAGCTCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGC‑ATCCCATACCCCGCAGCACACAGGGCT‑‑GCGGGGtat                                                                                                                                     <  6:280588/149‑3 (MQ=255)               aagaagTTAACTACGAAAAGAGTCGCCTGCAGGTTGCGGTGCTCATTTTCGGTCGCTCTACTCCGGTGGAGCTCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGC‑ATCCCATACCCCGCAGCCCTGTGTGCT‑‑GCGGGGtttt                                                                                                                                    >  1:240051/1‑149 (MQ=255)                agaagTTAACTACGAAAAGAGTCGCCTGCAGGTTGCGGTGCTCATTTTCGGTCGCTCTACTCCGGTGGAGCTCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGC‑ATCCCATACCCCGCAGCCCTGTGTGCT‑‑GCGGGGTTTTg                                                                                                                                   <  7:326246/149‑1 (MQ=255)                agaagTTAACTACGAAAAGAGTCGCCTGCAGGTTGCGGTGCTCATTTTCGGTCGCTCTACTCCGGTGGAGCTCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGC‑ATCCCATACCCCGCAGCACACAGGGCT‑‑GCGGGGtatgg                                                                                                                                   <  6:383981/149‑5 (MQ=255)                 gaagTTAACTACGAAAAGAGTCGCCTGCAGGTTGCGGTGCTCATTTTCGGTCGCTCTACTCCGGTGGAGCTCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGC‑ATCCCATACCCCGCAGCCCTGTGTGCT‑‑GCGGGGTTTTgt                                                                                                                                  >  3:123986/1‑149 (MQ=255)                   agTTAACTACGAAAAGAGTCGCCTGCAGGTTGCGGTGCTCATTTTCGGTCGCTCTACTCCGGTGGAGCTCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGC‑ATCCCATACCCCGCAGCCCTGTGTGCT‑‑GCGGGGTTTTgtc                                                                                                                                 <  2:86598/148‑1 (MQ=255)                    gTTAACTACGAAAAGAGTCGCCTGCAGGTTGCGGTGCTCATTTTCGGTCGCTCTACTCCGGTGGAGCTCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGC‑ATCCCATACCCCGCAGCCCTGTGTGCT‑‑GCGGGGTTTTgtcgt                                                                                                                               >  3:146563/1‑149 (MQ=255)                       aaCTACGAAAAGAGTCGCCTGCAGGTTGCGGTGCTCATTTTCGGTCGCTCTACTCCGGTGGAGCTCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGC‑ATCCCATACCCCGCAGCCCTGTGTGCT‑‑GCGGGGTTTTGTCGTCAc                                                                                                                            >  6:264166/1‑149 (MQ=255)                         cTACGAAAAGAGTCGCCTGCAGGTTGCGGGGCTCATTTTCGGTCGCTCTACTCCGGTGGAGCTCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGC‑ATCCCATACCCCGCAGCCCTGTGTGCT‑‑GCGGGGTTTTGTCGTCACTg                                                                                                                          <  5:215068/149‑1 (MQ=255)                           aCGAAAAGAGTCGCCTGCAGGTTGCGGTGCTCATTTTCGGTCGCTCTACTCCGGTGGAGCTCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGC‑ATCCCATACCCCGCAGCCCTGAGTGCT‑‑GCGGGGTTTTGTCGTCACTgg                                                                                                                         >  2:107462/1‑148 (MQ=255)                            cGAAAAGAGTCGCCTGCAGGTTGCGGTGCTCATTTTCGGTCGCTCTACTCCGGTGGAGCTCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGC‑ATCCCATACCCCGCAGCCCTGTGTGCT‑‑GCGGGGTTTTGTCGTCACTGGGa                                                                                                                       <  1:337331/149‑1 (MQ=255)                            cGAAAAGAGTCGCCTGCAGGTTGCGGTGCTCATTTTCGGTCGCTCTACTCCGGTGGAGCTCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGC‑ATCCCATACCCCGCAGCCCTGTGTGCT‑‑GCGGGGTTTTGTCGTCACTGGGa                                                                                                                       <  4:244082/149‑1 (MQ=255)                             gAAAAGAGTCGCCTGCAGGTTGCGGTGCTCATTTTCGGTCGCTCTACTCCGGTGGAGCTCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGC‑ATCCCATACCCCGCAGCCCTGTGTGCT‑‑GCGGGGTTTTGTCGTCACTGGGa                                                                                                                       >  1:100565/1‑148 (MQ=255)                             gAAAAGAGTCGCCTGCAGGTTGCGGTGCTCATTTTCGGTCGCTCTACTCCGGTGGAGCTCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGC‑ATCCCATACCCCGCAGCCCTGTGTGCT‑‑GCGGGGTTTTGTCGTCACTGGGa                                                                                                                       >  6:181975/1‑148 (MQ=255)                              aaaaGAGTCGCCTGCAGGTTGCGGTGCTCATTTTCGGTCGCTCTACTCCGGTGGAGCTCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGC‑ATCCCATACCCCGCAGCCCTGTGTGCT‑‑GCGGGGTTTTGTCGTCACTGGGAt                                                                                                                      >  2:316686/1‑148 (MQ=255)                              aaaaGAGTCGCCTGCAGGTTGCGGTGCTCATTTTCGGTCGCTCTACTCCGGTGGAGCTCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGC‑ATCCCATACCCCGCAGCCCTGTGTGCT‑‑GCGGGGTTTTGTCGTCACTGGGAt                                                                                                                      >  2:303504/1‑148 (MQ=255)                                       gCCTGCAGGTTGCGGTGCTCATTTTCGGTCGCTCTACTCCGGTGGAGCTCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGC‑ATCCCATACCCCGCAGCCCTGTGTGCT‑‑GCGGGGTTTTGTCGTCACTGGGATAAAACGCaa                                                                                                             <  6:166704/148‑1 (MQ=255)                                       gCCTGCAGGTTGCGGTGCTCATTTTCGGTCGCTCTACTCCGGTGGAGCTCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGC‑ATCCCATACCCCGCAGCCCTGTGTGCT‑‑GCGGGGTTTTGTCGTCACTGGGATAAAACGCAAg                                                                                                            <  2:271345/149‑1 (MQ=255)                                          tGCAGGTTGCGGTGCTCATTTTCGGTCGCTCTACTCCGGTGGAGCTCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGC‑ATCCCATACCCCGCAGCCCTGTGTGCT‑‑GCGGGGtt                                                                                                                                      >  3:180367/1‑120 (MQ=255)                                          tGCAGGTTGCGGTGCTCATTTTCGGTCGCTCTACTCCGGTGGAGCTCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGC‑ATCCCATACCCCGCAGCCCTGTGTGCT‑‑GCGGGGtt                                                                                                                                      <  4:180367/120‑1 (MQ=255)                                            cAGGTTGCGGTGCTCATTTTCGGTCGCTCTACTCCGGTGGAGCTCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGC‑ATCCCATACCCCGCAGCCCTGTGTGCT‑‑GCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATc                                                                                                       >  2:303361/1‑149 (MQ=255)                                                ttGCGGTGCTCATTTTCGGTCGCTCTACTCCGGTGGAGCTCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGC‑ATCCCATACCCCGCAGCCCTGTGTGCT‑‑GCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCggg                                                                                                   <  3:326415/149‑1 (MQ=255)                                                ttGCGGTGCTCATTTTCGGTCGCTCTACTCCGGTGGAGCTCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGC‑ATCCCATACCCCGCAGCCCTGTGTGCT‑‑GCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCggg                                                                                                   <  1:316686/149‑1 (MQ=255)                                                        cTCATTTTCGGTCGCTCTACTCCGGTGGAGCTCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGC‑ATCCCATACCCCGCAGCCCTGTGTGCT‑‑GCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCGGGGAGCCAt                                                                                            <  8:197159/148‑1 (MQ=255)                                                         tCATTTTCGGTCGCTCTACTCCGGTGGAGCTCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGC‑ATCCCATACCCCGCAGCCCTGTGTGCT‑‑GCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCGGGGAGCCAt                                                                                            >  6:140233/1‑147 (MQ=255)                                                           aTTTTCGGTCGCTCTACTCCGGTGGAGCTCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGC‑ATCCCATACCCCGCAGCCCTGTGTGCT‑‑GCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCGGGGAGCCATCAgg                                                                                        >  1:348122/1‑149 (MQ=255)                                                           aTTTTCGGTCGCTCTACTCCGGTGGAGCTCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGC‑ATCCCATACCCCGCAGCCCTGTGTGCT‑‑GCGGGGTTTTGTCGTCACTGGCATAAAACGCAAGTCATACGGGGAGCCATCAgg                                                                                        >  1:349515/1‑149 (MQ=255)                                                             tttCGGTCGCTCTACTCCGGTGGAGCTCGAGTTCAGACAGGTCGAAAAGCTCTAGGCCGACGAAGC‑ATCCCATACCCAGCAGCCCTGTCTGCT‑‑GCGGCGTATTGTCGTCACAGGGATAAAACGCAAGTCATCCCGGGAGCCATCAGGCg                                                                                      >  7:263008/1‑149 (MQ=255)                                                                 ggTCGCTCTACTCCGGTGGAGCTCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGC‑ATCCCATACCCCGCAGCCCTGTGTGCT‑‑GCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCGGGGAGCCATCAGGCGTTc                                                                                   <  4:123986/148‑1 (MQ=255)                                                                 ggTCGCTCTACTCCGGTGGAGCTCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGC‑ATCCCATACCCCGCAGCCCTGTGTGCT‑‑GCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCGGGGAGCCATCAGGCGTTc                                                                                   <  6:57597/148‑1 (MQ=255)                                                                  gTCGCTCTACTCCGGTGGAGCTCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGC‑ATCCCATACCCCGCAGCCCTGTGTGCT‑‑GCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCGGGGAGCCATCAGGCGTTCg                                                                                  >  7:94581/1‑148 (MQ=255)                                                                       tctACTCCGGTGGAGCTCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGC‑ATCCCATACCCCGCAGCCCTGTGTGCT‑‑GCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCGGGGAGCCATCAGGCGTTCGTAccc                                                                             >  3:211861/1‑148 (MQ=255)                                                                            tCCGGTGGAGCTCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGC‑ATCCCATACCCCGCAGCCCTGTGTGCT‑‑GCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGACATCCGGGGAGCCATCAGGCGTTCTTACCCGAtatt                                                                       >  2:128815/1‑149 (MQ=255)                                                                              cGGTGGAGCTCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGC‑ATCCCATACCCCGCAGCCCTGTGTGCT‑‑GCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCGGGGAGCCATCAGGCGTTCGTACCCGATTTTgg                                                                     <  5:264166/149‑1 (MQ=255)                                                                              cGGTGGAGCTCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGC‑ATCCCATACCCCGCAGCCCTGTGTGCT‑‑GCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCGGGGAGCCATCAGGCGTTCGTACCCGATTTTg                                                                      <  2:240051/148‑1 (MQ=255)                                                                              cGGTGGAGCTCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGC‑ATCCCATACCCCGCAGCCCTGTGTGCT‑‑GCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCGGGGAGCCATCAGGCGTTCGTACCCGATTTTg                                                                      <  6:326774/148‑1 (MQ=255)                                                                              cGGTGGAGCTCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGC‑ATCCCATACCCCGCAGCCCTGTGTGCT‑‑GCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCGGGGAGCCATCAGGCGTTCGTACCCGATTTTg                                                                      <  2:348122/148‑1 (MQ=255)                                                                              cGGTGGAGCTCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGC‑ATCCCATACCCCGCAGCCCTGTGTGCT‑‑GCGGGGTGTTGTCGTCACTGGGATAAAACGCAAGTCATCCGGGGAGACATCAGGCGTTCGTACGCGAGTTTg                                                                      >  8:327033/1‑148 (MQ=255)                                                                              cGGTGCAGCGCGAGTTCAGCCAGGGCGAAAAGGTCTAGGCGGACGAAGC‑ATCCCATACCCCGCAGCCCTGTGTGCT‑‑GCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCGGGGAGCCATCAGGCGTTCGTACCCGAtttt                                                                       <  4:51899/147‑1 (MQ=255)                                                                               ggTGGAGCTCGATTTCAGCCCGGGGTAAAAGGTCTTGGGGCACTAACA‑ATCCCATACCCCGCAGCCCTGTTTGCTGGGCCCCGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCGGGGAGCCATCAGGCGTTCGTACCCGATTTTg                                                                      <  1:303504/149‑1 (MQ=255)                                                                               ggTGGAGCTCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGC‑ATCCCATACCCCGCAGCCCTGTGTGCT‑‑GCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCGGGGAGCCATCAGGCGTTCGTACCCGATTTTGGa                                                                    <  6:141384/149‑1 (MQ=255)                                                                               ggTGGAGCTCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGC‑ATCCCATACCCCGCAGCCCTGTGTGCT‑‑GCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCGGGGAGCCATCAGGCGTTCGTACCCGATTTTGGa                                                                    <  8:94581/149‑1 (MQ=255)                                                                               ggTGGAGCTCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGC‑ATCCCATACCCCGCAGCCCTGTGTGCT‑‑GCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCGGGGAGCCATCAGGCGTTCGTACCCGATTTTGGa                                                                    <  7:262324/149‑1 (MQ=255)                                                                               ggTGGAGCTCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGC‑ATCCCATACCCCGCAGCCCTGTGTGCT‑‑GCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCGGGGAGCCATCAGGCGTTCGTACCCGATTTTGGa                                                                    <  1:136532/149‑1 (MQ=255)                                                                               ggTGGAGCTCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGC‑ATCCCATACCCCGCAGCCCTGTGTGCT‑‑GCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCGGGGAGCCATCAGGCGTTCGTACCCGATTTTGGa                                                                    <  2:100565/149‑1 (MQ=255)                                                                               ggTGGAGCTCGAGTTCAGCAAGGTCGAAAAGGTCTAGGCGGACGAAGC‑ATCCCATACCCCGCAGCCCTGTGTGCT‑‑GCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCGGGGAGCCATCAGGCGTTCGTACCCGATTTTGGa                                                                    >  4:306281/1‑149 (MQ=255)                                                                                   gAGCTCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGC‑ATCCCATACCCCGCAGCCCTGTGTGCT‑‑GCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCGGGGAGCCATCAGGCGTTCGTACCCGATTTTGGAGTAg                                                                >  3:108860/1‑149 (MQ=255)                                                                                     gCTCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGC‑ATCCAATACCCCGCAGCCCTGTGTGCT‑‑GCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCGGGGAGCCATCAGGCGTTCGTACCCGATTTTGGAGTAGc                                                               >  7:203689/1‑148 (MQ=255)                                                                                        cGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGC‑ATCCCATACCCCGCAGCCCTGTGTGCT‑‑GCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCGGGGAGCCATCAGGCGTTCGTACCCGAtttt                                                                       >  8:43396/1‑137 (MQ=255)                                                                                        cGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGC‑ATCCCATACCCCGCAGCCCTGTGTGCT‑‑GCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCGGGGAGCCATCAGGCGTTCGTACCCGAtttt                                                                       <  7:43396/137‑1 (MQ=255)                                                                                          aGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGC‑ATCCCATACCCCGCAGCCCTGTGTGCT‑‑GCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCGGGGAGCCATCAGGCGTTCGTACCCGATTTTGGAGTAGCTCATgg                                                         <  3:88785/149‑1 (MQ=255)                                                                                           gTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGCGATCCCATACCCCGCAGCCCTGTGTGCT‑‑GCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCGGGGAGCCATCAGGCGTTCGTACCCGATTTTGGAGTAGCTCATgg                                                         >  1:249545/1‑149 (MQ=255)                                                                                               aGCCAGGTCGAAAAGGTCTAGGCGGACGAAGC‑ATCCCATACCCCGCAGCCCTGTGTGCT‑‑GCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCGGGGAGCCATCAGGCGTTCGTACCCGATTTTg                                                                      <  6:129148/131‑1 (MQ=255)                                                                                               aGCCAGGTCGAAAAGGTCTAGGCGGACGAAGC‑ATCCCATACCCCGCAGCCCTGTGTGCT‑‑GCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCGGGGAGCCATCAGGCGTTCGTACCCGATTTTg                                                                      >  5:129148/1‑131 (MQ=255)                                                                                                 ccAGGTCGAAAAGGTCTAGGCGGACGAAGC‑ATCCCATACCCCGCAGCCCTGTGTGCT‑‑GCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCGGGGAGCCATCAGGCGTTCGTACCCGATTTTGGAGTAGCTCATGGCTaag                                                    >  7:221158/1‑147 (MQ=255)                                                                                                 ccAGGTCGAAAAGGTCTAGGCGGACGAAGC‑ATCCCATACCCCGCAGCCCTGTGTGCT‑‑GCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCGGGGAGCCATCAGGCGTTCGTACCCGATTTTGGAGTAGCTCATGGCTaag                                                    <  8:221158/147‑1 (MQ=255)                                                                                                 ccAGGTCGAAAAGGTCTAGGCGGACGAAGC‑ATCCCATACCCCGCAGCCCTGTGTGCT‑‑GCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCGGGGAGCCATCAGGCGTTCGTACCCGATTTTGGAGTAGCTCATGGCTaag                                                    >  3:333966/1‑147 (MQ=255)                                                                                                  cAGGTCGAAAAGGTCTAGGCGGACGAAGC‑ATCCCATACCCCGCAGCCCTGTGTGCT‑‑GCGGGGTTTTGTCGGCACTGGGATAAAACGCAAGTCATCCGGGGAGCCCTCAGGCGTTCGTAACCGATTTTGGAGTAGCTCATTGGTaag                                                    >  7:250294/1‑146 (MQ=255)                                                                                                   aGGTCGAAAAGGTCTAGGCGGACGAAGC‑ATCCCATACCCCGCAGCCCTGTGTGCT‑‑GCGGGGTTTTGTCGTCACTGGGATAAAACGCGAGTCATCCGGGGAGCCATCAGGCGTTCGTACCCGATTTTGGAGTAGCTCATGGCTaag                                                    >  7:196925/1‑145 (MQ=255)                                                                                                     gTCGAAAAGGTCTAGGCGGACGAAGC‑ATCCCATACCCCGCAGCCCTGTGTGCT‑‑GCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCGGGGAGCCATCAGGCGTTCGTACCCGATTTTGGAGTAGCTCATGGCTAAGAAGAtt                                              <  2:68025/149‑1 (MQ=255)                                                                                                      tCGAAAAGGTCTAGGCGGACGAAGC‑ATCCCATACCCCGCAGCCCTGTGTGCT‑‑GCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCGGGGAGCCATCAGGCGTTCGTACTCGATTTTGGAGTAGCTCATGGCTAAGAAGAtt                                              >  3:363921/1‑148 (MQ=255)                                                                                                        gAAAAGGTCTAGGCGGACGAAGC‑ATCCCATACCCCGCAGCCCTGTGTGCT‑‑GCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCGGGGAGCCATCAGGCGTTCGTACCCGATTTTGGAGTAGCTCATGGCTAAGAAGATTCa                                            <  4:211861/148‑1 (MQ=255)                                                                                                         aaaaGGTCTAGGCGGACGAAGC‑ATCCCATACCCCGCAGCCCTGTGTGCT‑‑GCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCGGGGAGCCATCAGGCGTTCGTACCCGATTTTGGAGTAGCTCATGGCTAAGAAGATTCAgg                                          <  7:90824/149‑1 (MQ=255)                                                                                                          aaaGGTCTAGGCGGACGAAGC‑ATCCCGTACCCCGCAGCCCTGTGTGCT‑‑GCGGGGTTTTGTAGTCACTGGGATAAAACGCAAGTCATCCGGGGAGCCATCAGGCGTTCGTACCCGATTTTGGAGTAGCTCATGGCTAAGAAGATTCAGGc                                         >  1:9253/1‑149 (MQ=255)                                                                                                          aaaGGTCTAGGCGGACGAAGC‑ATCCCATACCCCGCAGCCCTGAGAGCT‑‑GCGGGGAATTGTCGTCACTGGGATAAAACGCAACTCATCCGTGGAGACATCAGGCGTTCGTACCCGATTTTGAAGTAGCTCATGGCTAAGAAGATTCATGc                                         >  1:57676/1‑149 (MQ=255)                                                                                                             ggTCTAGGCGGACGAAGC‑ATCCCATACCCCGCAGCCCTGTGTGCT‑‑GCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCGGGGAGCCATCAGGCGTTCGTACCCGATTTTGGAGTAGCTCATGGCTAAGAAGATTCAGGCTTa                                      <  8:196925/149‑1 (MQ=255)                                                                                                                 tAGGCGGACGAAGC‑ATCCCATACCCCGCAGCCCTGTGTGCT‑‑GCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCGGGGAGCCATCAGGCGTTCGTACCCGTTTTTGGAGTAGCTCATGGCTAAGAAGATTCAGGCTTACAt                                   <  2:349515/148‑1 (MQ=255)                                                                                                                        aCGAAGC‑ATCCCATACCCCGCAGCCCTGTGTGCT‑‑GCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCGGGGAGCCATCAGGCGTTCGTACCCGATTTTGGAGTAGCTCATGGCTAAGAAGATTCAGGCTTACATCAAGCTGc                           <  3:43521/149‑1 (MQ=255)                                                                                                                        aCGAAGC‑ATCCCATACCCCGCAGCCCTGTGTGCT‑‑GCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCGGGGAGCCATCAGGCGTTCGTACCCGATTTTGGAGTAGCTCATGGCTAAGAAGATTCAGGCTTACATCAAGCTGc                           <  4:324834/149‑1 (MQ=255)                                                                                                                         cGAAGC‑ATCCCATACCCCGCAGCCCTGTGTGCT‑‑GCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCGGGGAGCCATCAGGCGTTCGTACCCGATTTTGGAGTAGCTCATGGCTAAGAAGATTCAGGCTTACATCAAGCTGc                           >  6:209360/1‑148 (MQ=255)                                                                                                                          gAAGC‑ATCCCATACCCCGCAGCCCTGTGTGCT‑‑GCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCGGGGAGCCATCAGGCGTTCGTACCCGATTTTGGAGTAGCTCATGGCTAAGAAGATTCAGGCTTACATCAAGCTGCaa                         <  5:209360/149‑1 (MQ=255)                                                                                                                                aTCCAATACCCCGCAGCCCTGTGTGCT‑‑GCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCGGGGAGCCATCAGGCGTTCGTACCCGATTTTGGAGTAGCTCATGGCTAAGAAGATTCAGGCTTACATCAAGCTGCAAGTTaa                    <  8:203689/149‑1 (MQ=255)                                                                                                                                  cccATACCCCGCAGCCCTGTGTGCT‑‑GCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCGGGGAGCCATCAGGCGTTCGTACCCGATTTTGGAGTAGCTCATGGCTAAGAAGATTCAGGCTTACATCAAGCTGCAAGTTAAgg                  >  5:50801/1‑149 (MQ=255)                                                                                                                                   ccATACCCCTCAGGCCTGTGTGCT‑‑GTGGGGTTTTGTCGGCACTGGGATAAATCGCAAGTCATCCGGGGAGCCATCAGGCGTTCGTACCCGATTTTGGAGTAGCTCATGGCTAAGAAGATTCAGGCTTAGATCAAGCTGCAAGTTAAgg                  <  4:146563/148‑1 (MQ=255)                                                                                                                                   ccATACCCCGCAGCCCTGTGTGCT‑‑GCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCGGGGAGCCATCAGGCGTTCGTACCCGATTTTGGAGTAGCTCATGGCTAAGAAGATTCAGGCTTACATCAAGCTGCAAGTTAAggc                 >  5:259104/1‑149 (MQ=255)                                                                                                                                        ccccGCAGCCCTGTGTGCT‑‑GCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCGGGGAGCCATCAGGCGTTCGTACCCGATTTTGGAGTAGCTCATGGCTAAGAAGATTAAGGCTTACATAAAGCTGCAAGTTAAgtccggc             >  6:344879/1‑148 (MQ=255)                                                                                                                                         cccGCAGCCCTGTGTGCT‑‑GCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCGGGGAGCCATCAGGCGTTCGTACCCGATTTTTTAGTAGCTCATGGCTAAGAAGATTCAGGCTTACATCAAGCTGCAAGTTAAGGCCggcca           <  2:57676/149‑1 (MQ=255)                                                                                                                                         cccGCAGCCCTGTGTGCT‑‑GCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCGGGGAGCCATCAGGCGTTCGTACCCGATTTTGGAGTAGCTCATGGCTAAGAAGATTCAGGCTTACATCAAGCTGCAAGTTAAggccggc             >  7:18726/1‑147 (MQ=255)                                                                                                                                         cccGCAGCCCTGTGTGCT‑‑GCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCGGGGAGCCATCAGGCGTTCGTACCCGATTTTGGAGTAGCTCATGGCTAAGAAGATTCAGGCTTACATCAAGCTGCAAGTTAAggccg               <  7:329102/145‑1 (MQ=255)                                                                                                                                         cccGCAGCCCTGTGTGCT‑‑GCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCGGGGAGCCATCAGGCGTTCGTACCCGATTTTGGAGTAGCTCATGGCTAAGAAGATTCAGGCTTACATCAAGCTGCAAGTTAAggccg               >  8:329102/1‑145 (MQ=255)                                                                                                                                          acGCAGCACTGTGTGCT‑‑GCGGGGTTTTGTCGTCACTGGGATAAAACTCAAGTCATCCGGGGAGCCATCAGTCGTTCGTACCCGATTTTGGAGTAGCTCATGGCTAAGAAGATTCAGGCTTACATCAAGCTGCAAGTTAAGGCCggccag          <  5:344879/148‑1 (MQ=255)                                                                                                                                           cGCAGCCCTGTGTGCT‑‑GCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCGGGGAGCCATCAGGCGTTCGTACCCGATTTTGGAGTAGCTCATGGCTAAGAAGATTCAGGCTTACATCAAGCTGCAAGTTAAGGCCggccagg         <  2:9253/149‑1 (MQ=255)                                                                                                                                                  cTGTGTGCT‑‑GCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCGGGGAGCCATCAGGCGTTCGTACCCGATTTTGGAGTAGCTCATGGCTAAGAAGATTCAGGCTTACATCAAGCTGCAAGTTAAGGCCtggccggccttaac  <  5:181975/149‑6 (MQ=255)                                                                                                                                                     |                                                                                                                                                    CAATGGTAGTGTCGAAGAAGTTAACTACGAAAAGAGTCGCCTGCAGGTTGCGGTGCTCATTTTCGGTCGCTCTACTCCGGTGGAGCTCGAGTTCAGCCAGGTCGAAAAGGTCTAGGCGGACGAAGC‑ATCCCATACCCCGCAGCCCTGTGTGCT‑‑GCGGGGTTTTGTCGTCACTGGGATAAAACGCAAGTCATCCGGGGAGCCATCAGGCGTTCGTACCCGATTTTGGAGTAGCTCATGGCTAAGAAGATTCAGGCTTACATCAAGCTGCAAGTTAAGGCCGGCCAGGCCAACCC  >  NC_002947/534901‑535193

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