Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,633,811	C→A	7.6%	A434S (GCA→TCA) ‡	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,811	0	C	A	7.6%	200.9 / 4.5	64	A434S (GCA→TCA) ‡	uvrC	UvrABC system protein C
Reads supporting (aligned to +/- strand):  ref base C (29/30);  new base A (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.74e-01

CCAGCTCCTGCATCACATCGCGGGCCATGTTCAGCTGGCCCTTGCCGCCGTCCACCAGCAGCACATCGGGCAGCTTGCCCTCGCCGTCCTTGATGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCGGGGCCGAACACCACGCAACTGGCCACGGTAGCTTCACCACTGGAGTGGCTGATGTCGTAGCATTCCAGGCGTTGCGGCACTTCATCCAGGCCC  >  NC_002947/4633664‑4633956                                                                                                                                                    |                                                                                                                                                  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                                                                                                                                       <  2:249873/149‑1 (MQ=255)            aTCACATCGCGGGCCATGTTCAGCTGGCCCTTGCCGCCGTCCACCAGCAGCACATCGGGCAGCTTGCCCTCGCCGTCCTTGATGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTc                                                                                                                                       >  4:229407/1‑149 (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                                                                                                             >  6:358461/1‑140 (MQ=255)                                               gccgTCCACCAGCAGCACATCGGGCAGCTTGCCCTCGCCGTCCTTGATGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACt                                                                                                             <  5:358461/140‑1 (MQ=255)                                               gccgTCCACCAGCAGCACATCGGGCAGCTTGCCCTCGCCGTCCTTGATGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGcc                                                                                                    <  7:149682/149‑1 (MQ=255)                                                      accaGGAGCACATCGGGCAGCTTGCCCGCGCCGTCCTTGATGCGCCCATAGCGGCGGGGCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGCTGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCggg                                                                                             <  8:241510/149‑1 (MQ=255)                                                      accaGCAGCACATCGGGCAGCTTGCCCTCGCCGTCCTTGATGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCggg                                                                                             >  5:190464/1‑149 (MQ=255)                                                      accaGCAGCACATCGGGCAGCTTGCCCTCGCCGTCCTTGATGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCggg                                                                                             >  6:375807/1‑149 (MQ=255)                                                      accaGCAGCACATCGGGCAGCTTGCCCTCGCCGTCCTTGATGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCggg                                                                                             <  6:141616/149‑1 (MQ=255)                                                      accaGCAGCACATCGGGCAGCTTGCCCTCGCCGTCCTTGATGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCggg                                                                                             <  4:227459/149‑1 (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                               <  2:36057/149‑1 (MQ=255)                                                                                                                    ggCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCGGGGCCGAACACCACGCAACTGGCCACGGTAGCTTCACCACTGGAGTGGCTGATGTCGTAGCAtt                               <  5:323993/149‑1 (MQ=255)                                                                                                                       cTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCGGGGCCGAACACCACGCAACTGGCCACGGTAGCTTCACCACTGGAGTGGCTGATGTCGTAGCATTCCa                            <  3:149474/149‑1 (MQ=255)                                                                                                                        tGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCGGGGCCGAACACCACGCAACTGGCCACGGTAGCTTCACCACTGGAGTGGCTGATGTCGTAGCATTCCAg                           >  6:159593/1‑149 (MQ=255)                                                                                                                        tGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCGGGGCCGAACACCACGCAACTGGCCACGGTAGCTTCACCACTGGAGTGGCTGATGTCGTAGCATTCCAg                           >  5:85577/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               >  4:226966/6‑147 (MQ=255)                                                                                                                                      gaagacGTCACCGCAGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCGGGGCCGAACACCACGCAACTGGCCACGGTAGCTTCACCACTGGAGTGGCTGATGTCGTAGCATTCCAGGCGTTGCGGCAc               <  3:226966/142‑1 (MQ=255)                                                                                                                                         aTCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCGGGGCCGAACACCACGCAACTGGCCACAGTAGCTTCACCACTGGAGTGGCTGATGTCGTAGCATTCCAGGCGTTGCGGCACTTCAt          >  1:324387/1‑149 (MQ=255)                                                                                                                                                 cTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCGGGGCCGAACACCACGCAACTGGCCACGGTAGCTTCACCACTGGAGTGGCTGATGTCGTAGCATTCCAGGCGTTGCGGCACTTCATCCAGGccc  <  4:372512/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                  CCAGCTCCTGCATCACATCGCGGGCCATGTTCAGCTGGCCCTTGCCGCCGTCCACCAGCAGCACATCGGGCAGCTTGCCCTCGCCGTCCTTGATGCGCCCATAGCGGCGGGTCAGGGCCTGATGCATGGCCGCATAATCGTCACCTGCGGTGACGTCTTCGATGTTGAAGCGGCGGTAGTCGGACTTGATCGGCCCTTCGGGGCCGAACACCACGCAACTGGCCACGGTAGCTTCACCACTGGAGTGGCTGATGTCGTAGCATTCCAGGCGTTGCGGCACTTCATCCAGGCCC  >  NC_002947/4633664‑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: ‑