Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,633,810	G→C	7.8%	A434G (GCA→GGA) ‡	uvrC ←	UvrABC system protein C

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,633,810	0	G	C	7.8%	191.0 / 7.3	64	A434G (GCA→GGA) ‡	uvrC	UvrABC system protein C
Reads supporting (aligned to +/- strand):  ref base G (29/30);  new base C (3/2);  total (32/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.83e-01

GCCCAGCTCCTGCATCACATCGCGGGCCATGTTCAGCTGGCCCTTGCCGCCGTCCACCAGCAGCACATCGGGCAGCTTGCCCTCGCCGTCCTTGATGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCGGGGCCGAACACCACGCAACTGGCCACGGTAGCTTCACCACTGGAGTGGCTGATGTCGTAGCATTCCAGGCGTTGCGGCACTTCATCCAGGCCC  >  NC_002947/4633662‑4633956                                                                                                                                                     |                                                                                                                                                   gCCCAGCTCCTGCATCACATCGCGGGCCATGTTCAGCTGGCCCTTGCCGCCGTCCACCAGCAGCACATCGGGCAGCTTGCCCTCGCCGTCCTTGATGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTg                                                                                                                                                    >  6:287373/1‑149 (MQ=255)   ccAGCTCCTGCATCACATCGCGGGCCATGTTCAGCTGGCCCTTGCCGCCGTCCACCAGCAGCACATCGGGCAGCTTGCCCTCGCCGTCCTTGATGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCg                                                                                                                                                  <  5:31525/149‑1 (MQ=255)   ccAGCCCCTGCATCACATCGCGGGCCATGTTCAGCTGGCCCTTGCCGCCGTCCACCAGCAGCACATCGGGCAGCTTGCCCTCGCCGTCCTTGATGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCg                                                                                                                                                  <  7:224754/149‑1 (MQ=255)          cTGCATCACATCGCGGGCCATGTTCAGCTGGCCCTTGCCGCCGTCCACCAGCAGCACATCGGGCAGCTTGCCCTCGCCGTCCTTGATGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGAc                                                                                                                                             <  2:220114/147‑1 (MQ=255)          cTGCATCACATCGCGGGCCATGTTCAGCTGGCCCTTGCCGCCGTCCACCAGCAGCACATCGGGCAGCTTGCCCTCGCCGTCCTTGATGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCGCAGGTGACga                                                                                                                                           <  3:193655/149‑2 (MQ=255)           tGCATCACATCGCGGGCCATGTTCAGCTGGCCCTTGCCGCCGTCCACCAGCAGCACATCGGGCAGCTTGCCCTCGCCGTCCTTTATGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACggc                                                                                                                                          >  7:241510/1‑147 (MQ=255)           tGCATCACATCGCGGGCCATGTTCAGCTGGCCCTTGCCGCCGTCCACCAGCAGCACATCGGGCAGCTTGCCCTCGCCGTCCTTGATGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTc                                                                                                                                          >  1:36057/1‑149 (MQ=255)           tGCATCACATCGCGGGCCATGTTCAGCTGGCCCTTGCCGCCGTCCACCAGCAGCACATCGGGCAGCTTGCCCTCGCCGTCCTTGATGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTc                                                                                                                                          >  6:194926/1‑149 (MQ=255)             cATCACATCGCGGGCCATGTTCAGCTGGCCCTTGCCGCCGTCCACCAGCAGCACATCGGGCAGCTTGCCCTCGCCGTCCTTGATGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCtt                                                                                                                                        <  7:286371/149‑1 (MQ=255)              aTCACATCGCGGGCCATGTTCAGCTGGCCCTTGCCGCCGTCCACCAGCAGCACATCGGGCAGCTTGCCCTCGCCGTCCTTGATGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTc                                                                                                                                       >  4:229407/1‑149 (MQ=255)              aTCACATCGCGGGCCATGTTCAGCTGGCCCTTGCCGCCGTCCACCAGCAGCACATCGGGCAGCTTGCCCTCGCCGTCCTTGATGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTc                                                                                                                                       <  2:249873/149‑1 (MQ=255)               tCACATCGCGGGCCATGTTCAGCTGGCCCTTGCCGCCGTCCACCAGCAGCACATCGGGCAGCTTGCCCTCGCCGTCCTTGATGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTc                                                                                                                                       >  4:220295/1‑148 (MQ=255)                 acaTCGCGGGCCATGTTCAGCTGGCCCTTGCCGCCGTCCACCAGCAGCACATCGGGCAGCTTGCCCTCGCCGTCCTTGATGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGAt                                                                                                                                    <  8:89324/149‑1 (MQ=255)                        gggCCATGTTCAGCTGGCCCTTGCCGCCGTCCACCAGCAGCACATCGGGCAGCTTGCCCTCGCCGTCCTTGATGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAg                                                                                                                             <  8:91447/149‑1 (MQ=255)                           ccATGTTCAGCTGGCCCTTGCCGCCGTCCACCAGCAGCACATCGGGCAGCTTGCCCTCGCCGTCCTTGATGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAgcgg                                                                                                                          >  3:227459/1‑149 (MQ=255)                             aTGTTCAGCTGGCCCTTGCCGCCGTCCACCAGCAGCACTTCGGGCAGCTTGCCCTCGCCGTCCTTGATGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGTAgcggcg                                                                                                                        <  2:95317/149‑1 (MQ=255)                                 tCAGCTGGCCCTTGCCGCCGTCCACCAGCAGCACATCGGGCAGCTTGCCCTCGCCGTCCTTGATGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAg                                                                                                                    >  5:141616/1‑149 (MQ=255)                                 tCAGCTGGCCCTTGCCGCCGTCCACCAGCAGCACATCGGGCAGCTTGCCCTCGCCGTCCTTGATGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAg                                                                                                                    >  1:360797/1‑149 (MQ=255)                                     cTGGCCCTTGCCGCCGTCCACCAGCAGCACATCGGGCAGCTTGCCCTCGCCGTCCTTGATGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCg                                                                                                                 >  8:149682/1‑148 (MQ=255)                                        gCCCTTGCCGCCGTCCACCAGCAGCACATCGGGCAGCTTGCCCTCGCCGTCCTTGATGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGTTGTTGAAGCGGCGGTAGTCGGACt                                                                                                             <  7:206734/149‑1 (MQ=255)                                        gCCCTTGCCGCCGTCCACCAGCAGCACATCGGGCAGCTTGCCCTCGCCGTCCTTGATGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGAc                                                                                                              <  3:229407/148‑1 (MQ=255)                                                cgccgTCCACCAGCAGCACATCGGGCAGCTTGCCCTCGCCGTCCTTGATGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGc                                                                                                     >  3:33739/1‑149 (MQ=255)                                                 gccgTCCACCAGCAGCACATCGGGCAGCTTGCCCTCGCCGTCCTTGATGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACt                                                                                                             <  5:358461/140‑1 (MQ=255)                                                 gccgTCCACCAGCAGCACATCGGGCAGCTTGCCCTCGCCGTCCTTGATGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACt                                                                                                             >  6:358461/1‑140 (MQ=255)                                                 gccgTCCACCAGCAGCACATCGGGCAGCTTGCCCTCGCCGTCCTTGATGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGcc                                                                                                    <  7:149682/149‑1 (MQ=255)                                                        accaGGAGCACATCGGGCAGCTTGCCCGCGCCGTCCTTGATGCGCCCATAGCGGCGGGGCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGCTGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCggg                                                                                             <  8:241510/149‑1 (MQ=255)                                                        accaGCAGCACATCGGGCAGCTTGCCCTCGCCGTCCTTGATGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCggg                                                                                             <  4:227459/149‑1 (MQ=255)                                                        accaGCAGCACATCGGGCAGCTTGCCCTCGCCGTCCTTGATGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCggg                                                                                             <  6:141616/149‑1 (MQ=255)                                                        accaGCAGCACATCGGGCAGCTTGCCCTCGCCGTCCTTGATGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCggg                                                                                             >  6:375807/1‑149 (MQ=255)                                                        accaGCAGCACATCGGGCAGCTTGCCCTCGCCGTCCTTGATGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCggg                                                                                             >  5:190464/1‑149 (MQ=255)                                                           agcagcACATCGGGCAGCTTGCCCTCGCCGTCCTTGATGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCGGGGcc                                                                                          <  5:331004/149‑1 (MQ=255)                                                               gcaCATCGGGCAGCTTGCCCTCGCCGTCCTTGATGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCGGGGCCGaa                                                                                       >  4:149474/1‑148 (MQ=255)                                                                 acaTCGGGCAGCTTGCCCTCGCCGTCCTTGATGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCGGGGCCGaa                                                                                       <  4:49993/146‑1 (MQ=255)                                                                 acaTCGGGCAGCTTGCCCTCGCCGTCCTTGATGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCGGGGCCGaa                                                                                       >  3:49993/1‑146 (MQ=255)                                                                 acaTCGGGCAGCTTGCCCTCGCCGTCCTTGATGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCGGGGCCGAacac                                                                                    >  3:73805/1‑149 (MQ=255)                                                                      gggCAGCTTGCCCTCGCCGTCCTTGATGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCGGGGCCGAACACCACg                                                                                >  3:372512/1‑148 (MQ=255)                                                                               gCCCTCGCCGTCCTTGATGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCGGGGCCGAACACCACGCAACTGGCCa                                                                      >  5:123829/1‑149 (MQ=255)                                                                                 ccTCGCCGTCCTTGATGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCGGGGCCGAACACCACGCAACTGGCCACg                                                                    >  8:156267/1‑149 (MQ=255)                                                                                      ccGTCCTTGATGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCGGGGCCGAACACCACGCAACTGGCCACGGTAGc                                                               <  4:33739/149‑1 (MQ=255)                                                                                         tCCTTGATGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCGGGGCCGAACACCACGCAACTGGCCACGGTAGCtt                                                             >  8:184795/1‑148 (MQ=255)                                                                                                     ccATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCGGGGCCGATCACCACGCAACTGGCCACGGTAGCTTCACCACTGGAGTg                                                >  7:171317/1‑149 (MQ=255)                                                                                                     ccATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCGGGGCCGAACACCACGCAACTGGCCACGGTAGCTTCACCACTGGAGTg                                                >  2:98238/1‑149 (MQ=255)                                                                                                      cATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCGGGGCCGAACACCACGCAACTGGCCACGGTAGCTTCACCACTGGAGTgg                                               <  5:194926/149‑1 (MQ=255)                                                                                                      cATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCGGGGCCGAACACCACGCAACTGGCCACGGTAGCTTCACCACTGGAGTgg                                               >  1:100642/1‑149 (MQ=255)                                                                                                       aTAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCGGGGCCGAACACCACGCAACTGGCCACGGTAGCTTCACCACTGGAGTGGc                                              <  4:73805/149‑1 (MQ=255)                                                                                                           cggcggGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCGGGGCCGAACACCACGCAACTGGCCACGGTAGCTTCACCACTGGAGTGGCTGAt                                          >  6:323993/1‑149 (MQ=255)                                                                                                                 gcaaGGGCCTGATGCATGGCCGCATAATCGTCACCTGCTGTGACGGCTGCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCGGGGCCGAACACCACGCAACTGGCCACGGTAGCTTCACCACTGGAGTGGCTGATGTCGTa                                    <  8:262339/146‑1 (MQ=255)                                                                                                                      ggCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCGGGGCCGAACACCACGCAACTGGCCACGGTAGCTTCACCACTGGAGTGGCTGATGTCGTAGCAtt                               <  5:323993/149‑1 (MQ=255)                                                                                                                      ggCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCGGGGCCGAACACCACGCAACTGGCCACGGTAGCTTCACCACTGGAGTGGCTGATGTCGTAGCAtt                               <  2:36057/149‑1 (MQ=255)                                                                                                                         cTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCGGGGCCGAACACCACGCAACTGGCCACGGTAGCTTCACCACTGGAGTGGCTGATGTCGTAGCATTCCa                            <  3:149474/149‑1 (MQ=255)                                                                                                                          tGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCGGGGCCGAACACCACGCAACTGGCCACGGTAGCTTCACCACTGGAGTGGCTGATGTCGTAGCATTCCAg                           >  5:85577/1‑149 (MQ=255)                                                                                                                          tGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCGGGGCCGAACACCACGCAACTGGCCACGGTAGCTTCACCACTGGAGTGGCTGATGTCGTAGCATTCCAg                           >  6:159593/1‑149 (MQ=255)                                                                                                                           gATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCGGGGCCGAACACCACGCAACTGGCCACGGTAg                                                                <  2:226360/111‑1 (MQ=255)                                                                                                                           gATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCGGGGCCGAACACCACGCAACTGGCCACGGTAg                                                                >  1:226360/1‑111 (MQ=255)                                                                                                                            aTGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCGGGGCCg                                                                                         <  3:133543/85‑1 (MQ=255)                                                                                                                            aTGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCGGGGCCg                                                                                         >  4:133543/1‑85 (MQ=255)                                                                                                                               cATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCGGGGCCGAACACCACGCAACTGGCCACGGTAGCTTCACCACTGGAGTGGCTGATGTCGTAGCATTCCAGGCGtt                      <  5:375807/149‑1 (MQ=255)                                                                                                                                 tGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCGGGGCCGAACACCACGCAACTGGCCACGGTAGCTTCACCACTGGAGTGGCTGATGTCGTAGCATTCCAGGCGTTGc                    <  7:184795/149‑1 (MQ=255)                                                                                                                                      gCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCGGGGCCGAACACCACGCAACTGGCCACGGTAGCTTCACCACTGGAGTGGCTGATGTCGTAGCATTCCAGGCGTTGCGGCAc               <  5:159593/149‑1 (MQ=255)                                                                                                                                       cgaagacGTCACCGCAGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGCCTTGATCGGCCCTTCGGGGCCGAACACCACGCAACTGGCCACGGTAGCTTCACCACTGGAGTGGCTGATGTCGTAGCATTCCAGGCGTTGCGGCACt              >  7:169196/7‑149 (MQ=255)                                                                                                                                       cgaagacGTCACCGCAGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCGGGGCCGAACACCACGCAACTGGCCACGGTAGCTTCACCACTGGAGTGGCTGATGTCGTAGCATTCCAGGCGTTGCGGCACt              >  8:312673/7‑149 (MQ=255)                                                                                                                                        gaagacGTCACCGCAGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCGGGGCCGAACACCACGCAACTGGCCACGGTAGCTTCACCACTGGAGTGGCTGATGTCGTAGCATTCCAGGCGTTGCGGCAc               <  3:226966/142‑1 (MQ=255)                                                                                                                                        gaagacGTCACCGCAGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCGGGGCCGAACACCACGCAACTGGCCACGGTAGCTTCACCACTGGAGTGGCTGATGTCGTAGCATTCCAGGCGTTGCGGCAc               >  4:226966/6‑147 (MQ=255)                                                                                                                                           aTCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCGGGGCCGAACACCACGCAACTGGCCACAGTAGCTTCACCACTGGAGTGGCTGATGTCGTAGCATTCCAGGCGTTGCGGCACTTCAt          >  1:324387/1‑149 (MQ=255)                                                                                                                                                   cTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCGGGGCCGAACACCACGCAACTGGCCACGGTAGCTTCACCACTGGAGTGGCTGATGTCGTAGCATTCCAGGCGTTGCGGCACTTCATCCAGGccc  <  4:372512/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                   GCCCAGCTCCTGCATCACATCGCGGGCCATGTTCAGCTGGCCCTTGCCGCCGTCCACCAGCAGCACATCGGGCAGCTTGCCCTCGCCGTCCTTGATGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCGGGGCCGAACACCACGCAACTGGCCACGGTAGCTTCACCACTGGAGTGGCTGATGTCGTAGCATTCCAGGCGTTGCGGCACTTCATCCAGGCCC  >  NC_002947/4633662‑4633956

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