Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,381,492	C→A	6.4%	Y81* (TAC→TAA) ‡	dinB →	DNA polymerase IV

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,381,492	0	C	A	6.4%	200.3 / 3.7	62	Y81* (TAC→TAA) ‡	dinB	DNA polymerase IV
Reads supporting (aligned to +/- strand):  ref base C (31/27);  new base A (2/2);  total (33/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.67e-01

GGTGGGGGGCTCGCCCGACCATCGCGGGGTAATCGCCACCTGCAACTATGAAGCGCGTGCCTATGGTGTGCGCTCGGCCATGTCATCGCGGCATGCGCTGAAGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTA‑CCGTG‑‑AGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCTTTGTCCCTGGACGAAGCTTATCTGGATGTCAGCGACAGCCAGTGGTATTCGGGCAGCGCCACGCGCATTGCCGAGGATATTCGCCG  >  NC_002947/1381345‑1381638                                                                                                                                                     |                                                                                                                                                     ggTGGGGGGCTCGCCCGACCATCGCGGGGTAATCGCCACCTGCAACTATGAAGCGCGTGCCTATGGTGTGCGCTCGGCCATGTCATCGCGGCATGCGCTGAAGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTA‑cc                                                                                                                                                     <  1:126871/149‑1 (MQ=255)     gggggCTCGCCCGACCATCGCGGGGTAATCGCCACCTGCAACTATGAAGCGCGTGCCTATGGTGTGCGCTCGGCCATGTCATCGCGGCATGCGCTGAAGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTA‑CCGTG‑‑a                                                                                                                                               <  2:319386/149‑1 (MQ=255)     gggggCTCGCCCGACCATCGCGGGGTAATCGCCACCTGCAACTATGAAGCGCGTGCCTATGGTGTGCGCTCGGCCATGTCATCGCGGCATGCGCTGAAGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTA‑CCGTG‑‑a                                                                                                                                               <  2:314238/149‑1 (MQ=255)      ggggCTCGCCCGACCATCGCGGGGTAATCGCCACCTGCAACTATGAAGCGCGTGCCTATGGTGTGCGCTCGGCCATGTCATCGCGGCATGCGCTGAAGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTA‑CCGTG‑‑Ag                                                                                                                                              >  8:162925/1‑149 (MQ=255)      ggggCTCGCCCGACCATCGCGGGGTAATCGCCACCTGCAACTATGAAGCGCGTGCCTATGGTGTGCGCTCGGCCATGTCATCGCGGCATGCGCTGAAGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTA‑CCGTG‑‑Ag                                                                                                                                              <  7:8920/149‑1 (MQ=255)         gCTCGCCCGACCATCGCGGGGTAATCGCCACCTGCAACTATGAAGCGCGTGCCTATGGTGTGCGCTCGGCCATGTCATCGCGGCATGCGCTGAAGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTA‑CCGTG‑‑AGGcc                                                                                                                                           >  7:156214/1‑149 (MQ=255)               ccGACCATCGCGGGGTAATCGCCACCTGCAACTATGAAGCGCGTGCCTATGGTGTGCGCTCGGCCATGTCATCGCGGCATGCGCTGAAGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTA‑CCGTG‑‑AGGCCTCGCgg                                                                                                                                     <  4:83477/149‑1 (MQ=255)               ccGACCATCGCGGGGTAATCGCCACCTGCAACTATGAAGCGCGTGCCTATGGTGTGCGCTCGGCCATGTCATCGCGGCATGCGCTGAAGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTA‑CCGTG‑‑AGGCCTCGCgg                                                                                                                                     >  8:14366/1‑149 (MQ=255)                 gACCATCGCGGGGTAATCGCCACCTGCAACTATGAAGCGCGTGCCTATGGTGTGCGCTCGGCCATGTCATCGCGGCATGCGCTGAAGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTA‑CCGTG‑‑AGGCCTcg                                                                                                                                        >  4:37214/1‑144 (MQ=255)                 gACCATCGCGGGGTAATCGCCACCTGCAACTATGAAGCGCGTGCCTATGGTGTGCGCTCGGCCATGTCATCGCGGCATGCGCTGAAGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTA‑CCGTG‑‑AGGCCTcg                                                                                                                                        <  3:37214/144‑1 (MQ=255)                                   gCCACCTGCAACTATGAAGCGCGTGCCTATGGTGTGCGCTCGGCCATGTCATCGCGGCATGCGCTGAAGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTA‑CCGTG‑‑AGGCCTCGCGGGAAATCCACACGATCTTTcg                                                                                                                 <  7:195921/149‑1 (MQ=255)                                    ccACCTGCAACTATGAAGCGCGTGCCTATGGTGTGCGCTCGGCCATGTCATCGCGGCATGCGCTGAAGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTA‑CCGTG‑‑AGGCCTCGCGGGAAATCCACACGATCTTTcgc                                                                                                                >  2:103098/1‑149 (MQ=255)                                      aCCTGCAACTATGAAGCGCGTGCCTATGGTGTGCGCTCGGCCATGTCATCGCGGCATGCGCTGAAGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTA‑CCGTG‑‑AGGCCTCGCGGGAAATCCACACGATCTTTCGCGa                                                                                                              <  7:162925/149‑1 (MQ=255)                                         tGCAACTATGAAGCGCGTGCCTATGGTGTGCGCTCGGCCATGTCATCGCGGCATGCGCTGAAGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTA‑CCGTG‑‑AGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTa                                                                                                           >  6:116124/1‑149 (MQ=255)                                            aaCTATGAAGCGCGTGCCTATGGTGTGCGCTCGGCCATGTCATCGCGGCATGCGCTGAAGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTA‑CCGTG‑‑AGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTacac                                                                                                        <  7:86643/149‑1 (MQ=255)                                            aaCTATGAAGCGCGTGCCTATGGTGTGCGCTCGGCCATGTCATCGCGGCATGCGCTGAAGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTA‑CCGTG‑‑AGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTacac                                                                                                        >  4:139663/1‑149 (MQ=255)                                                 tGAAGCGCGTGCCTATGGTGTGCGCTCGGCCATGTCATCGCGGCATGCGCTGAAGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTACCCGTG‑‑AGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGa                                                                                                     <  1:75780/148‑1 (MQ=255)                                                 tGAAGCGCGTGCCTATGGTGTGCGCTCGGCCATGTCATCGCGGCATGCGCTGAAGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTA‑CCGTG‑‑AGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGc                                                                                                   >  6:137433/1‑149 (MQ=255)                                                         gTGCCTATGGTGTGCGCTCGGCCATGTCATCGCGGCATGCGCTGAAGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTA‑CCGTG‑‑AGGCCTCGCGGGAAATCcaca                                                                                                                           <  6:202943/117‑1 (MQ=255)                                                         gTGCCTATGGTGTGCGCTCGGCCATGTCATCGCGGCATGCGCTGAAGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTA‑CCGTG‑‑AGGCCTCGCGGGAAATCcaca                                                                                                                           >  5:202943/1‑117 (MQ=255)                                                                       cgcTCGGCCATGTCATCGCGGCATGCGCTGAAGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTA‑CCGTGTCAGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCTTTGTCCCt                                                                                <  1:103098/148‑1 (MQ=255)                                                                        gcTCGGCCATGTCATCGCGGCATGCGCTGAAGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTA‑CCGTG‑‑AGGCCTCGCGGGAAATc                                                                                                                               >  6:192644/1‑98 (MQ=255)                                                                        gcTCGGCCATGTCATCGCGGCATGCGCTGAAGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTA‑CCGTG‑‑AGGCCTCGCGGGAAATc                                                                                                                               <  5:192644/98‑1 (MQ=255)                                                                        gcTCGGCCATGTCATCGCGGCATGCGATGAAGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTA‑CCGTG‑‑AGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCTTTGTCCCTGGAc                                                                            >  1:182920/1‑149 (MQ=255)                                                                          tCGGCCATGTCATCGCGGCATGCGCTGAAGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTA‑CCGTG‑‑AGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCTTTGTCCCTGGACg                                                                           >  7:26193/1‑148 (MQ=255)                                                                          tCGGCCATGTCATCGCGGCATGCGCTGAAGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTA‑CCGTG‑‑AGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCTTTGTCCCTGGACg                                                                           >  6:88350/1‑148 (MQ=255)                                                                            ggCCATGTCATCGCGGCATGCGCTGAAGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTA‑CCGTG‑‑AGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCTTTGTCCCTGGACGaa                                                                         <  6:207688/148‑1 (MQ=255)                                                                                           gCATGCGCTGAAGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTA‑CCGTG‑‑AGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCTTTGTCCCTGGACGAAGCTTATCTGGATGTCa                                                         <  5:88350/149‑1 (MQ=255)                                                                                                cgcTGAAGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTA‑CCGTG‑‑AGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCTTTGTCCCTGGACGAAGCTTATCTGGATGTCAGCGa                                                     >  7:233666/1‑148 (MQ=255)                                                                                                 gcTGAAGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTA‑CCGTG‑‑AGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCTTTGTCCCTGGACGAAGCTTATCTGGATGTCAGCGa                                                     >  3:228831/1‑147 (MQ=255)                                                                                                 gcTGAAGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTA‑CCGTG‑‑AGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCTTTGTCCCTGGACGAAGCTTATCTGGATGTCAGCGACa                                                   <  5:137433/149‑1 (MQ=255)                                                                                                     aaGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTA‑CCGTG‑‑AGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCTTTGTCCCTGGACGAAGCTTATCTGGATGTCAGCGACAGc                                                 >  1:125136/1‑147 (MQ=255)                                                                                                     aaGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTA‑CCGTG‑‑AGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCTTTGTCCCTGGACGAAGCTTATCTGGATGTCAGCGACAGc                                                 <  8:26193/147‑1 (MQ=255)                                                                                                      aGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTA‑CCGTG‑‑AGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCTTTGTCCCTGGACGAAGCTTATCTGGATGTCAGCGACAGCCa                                               <  4:228831/148‑1 (MQ=255)                                                                                                        cTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTA‑CCGTG‑‑AGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCTTTGTCCCTGGACGAAGCTTATCTgg                                                               <  7:108390/130‑1 (MQ=255)                                                                                                        cTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTA‑CCGTG‑‑AGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCTTTGTCCCTGGACGAAGCTTATCTgg                                                               >  8:108390/1‑130 (MQ=255)                                                                                                              cccGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTA‑CCGTG‑‑AGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCTTTGTCCCTGGACGAAGCTTATCTGGATGTCAGCGACAGCCAGTGGTATTc                                      <  5:116124/149‑1 (MQ=255)                                                                                                                       ttGATCGTCAAGCCGCGCTTCGAGGCCTA‑CCGTG‑‑AGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCTTTGTCCCTGGACGAAGCTTATCTGGATGTCAGCGACAGCCAGTGGTATTCggg                                   <  3:333737/143‑1 (MQ=255)                                                                                                                         gATCGTCAAGCCGCGCTTCGAGGCCTA‑CCGTG‑‑AGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCTTTGTCCCTGGACGAAGCTTATCTGGATGTCAGCGACAGCCAGTGGTATTCGGGCAGCGCCa                           >  2:286905/1‑149 (MQ=255)                                                                                                                          aTCGTCAAGCCGCGCTTCGAGGCCTA‑CCGTG‑‑AGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCTTTGTCCCTGGACGAAGCTTATCTGGATGTCAGCGACAGCCAGTGGTATTCGGGCAGCGCCAc                          >  1:133973/1‑149 (MQ=255)                                                                                                                             gTCAAGCCGCGCTTCGAGGCCTA‑CCGTG‑‑AGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCTTTGTCCCTGGACGAAGCTTATCTGGATGTCAGCGACAGCCAGTGGTATTCGGGCAGCGCCAcgcg                       >  5:219814/1‑149 (MQ=255)                                                                                                                             gTCAAGCCGCGCTTCGAGGCCTA‑CCGTG‑‑AGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCTTTGTCCCTGGACGAAGCTTATCTGGATGTCAGCGACAGCCAGTGGTATTCGGGCAGCGCCAcgcg                       >  4:102138/1‑149 (MQ=255)                                                                                                                             gTCAAGCCGCGCTTCGAGGCCTA‑CCGTG‑‑AGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCTTTGTCCCTGGACGAAGCTTATCTGGATGTCAGCGACAGCCAGTGGTATTCGGGCAGCGCCAcgcg                       >  4:211108/1‑149 (MQ=255)                                                                                                                             gTCAAGCCGCGCTTCGAGGCCTA‑CCGTG‑‑AGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCTTTGTCCCTGGACGAAGCTTATCTGGATGTCAGCGACAGCCAGCGGTATTCGGGCAGCGCCAcgcg                       >  7:94137/1‑149 (MQ=255)                                                                                                                                aaGCCGCGCTTCGAGGCCTA‑CCGTG‑‑AGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCTTTGTCCCTGGACGAAGCTTATCTGGATGTCAGCGACAGCCAGTGGTATTCggg                                   >  3:71420/1‑134 (MQ=255)                                                                                                                                aaGCCGCGCTTCGAGGCCTA‑CCGTG‑‑AGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCTTTGTCCCTGGACGAAGCTTATCTGGATGTCAGCGACAGCCAGTGGTATTCggg                                   <  4:71420/134‑1 (MQ=255)                                                                                                                                aaGCCGCGCTTCGAGGCCTA‑CCGTG‑‑AGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCTTTGTCCCTGGACGAAGCTTATCTGGATGTCAGCGACAGCCAGTGGTATTCGGGCAgcgc                             <  6:336703/140‑1 (MQ=255)                                                                                                                                aaGCCGCGCTTCGAGGCCTA‑CCGTG‑‑AGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCTTTGTCCCTGGACGAAGCTTATCTGGATGTCAGCGACAGCCAGTGGTATTCGGGCAGCGCCACGCGCa                     <  8:233666/148‑1 (MQ=255)                                                                                                                                aaGCCGCGCTTCGAGGCCTA‑CAGTG‑‑AGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCTTTGTCCCTGGACGAAGCTTATCTGGATGTCAGCGACAGCCAGTGGGATTCGGGCAgcgc                             >  5:336703/1‑140 (MQ=255)                                                                                                                                       gcTTCGAGGCCTA‑CCGTG‑‑AGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCTTTGTCCCTGGACGAAGCTTATCTGGATGTCAGCGACAGCCAGTGGTATTCGGGCAGCGCCACGCGCATTGCCGAg             <  3:139663/149‑1 (MQ=255)                                                                                                                                         ttCGAGGCCTA‑CCGTG‑‑AGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCTTTGTCCCTGGACGAAGCTTATCTGGATGTCAGc                                                       <  7:152558/105‑1 (MQ=255)                                                                                                                                         ttCGAGGCCTA‑CCGTG‑‑AGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCTTTGTCCCTGGACGAAGCTTATCTGGATGTCAGc                                                       >  8:152558/1‑105 (MQ=255)                                                                                                                                         cgcGAGGCCTC‑ACGGT‑‑AGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCTTTGTCCCTGGACGAAGCTTATCTGGa                                                              <  3:180556/96‑1 (MQ=255)                                                                                                                                         cgcGAGGCCTC‑ACGGT‑‑AGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCTTTGTCCCTGGACGAAGCTTATCTGGa                                                              >  4:180556/3‑98 (MQ=255)                                                                                                                                          gcGAGGCCTC‑ACGGT‑‑AGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCTTTGTCCCTGGACGAAGCTTATCTGGATGTCAGCGACAGCCAGTg                                            >  3:324213/2‑115 (MQ=255)                                                                                                                                          gcGAGGCCTC‑ACGGT‑‑AGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCTTTGTCCCTGGACGAAGCTTATCTGGATGTCAGCGACAGCCAGTg                                            <  4:324213/114‑1 (MQ=255)                                                                                                                                           cGAGGCCTA‑CCGTG‑‑AGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCTTTGTCCCTGGACGAAGCTTATCTGGATGTCAGCGACAGCCAGTGGTATTCGGGCAGCGCCACGCGCATTGCCGAGGa           <  7:14366/147‑1 (MQ=255)                                                                                                                                               gCCTA‑CCGTG‑‑AGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCTTTGTCCCTGGACGAAGCTTATCTGGATGTCAGCGACAGCCAGTGGTATTCGGGCAGCGCCACGCGCATTGCCGAGGATATTCg     >  7:170639/1‑149 (MQ=255)                                                                                                                                                 cTA‑CCGTG‑‑AGGCCTCGCGGGAAATCCACATGATCTTTCGCGATTACACGGAGCTCATCGAGCCTTTGTCCCTGGACGAAGCTTATCTGGATGTCAGCGACAGCCAGTGGTATTCGGGCAGCGCCACGCGCATTGCCGAGGATATTCGcc   >  2:13499/1‑149 (MQ=255)                                                                                                                                                 cTA‑CCGTG‑‑AGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCTTTGTCCCTGGACGAAGCTTATCTGGATGTCAGCGACAGCCAGTGGTATTCGGGCAGCGCCACGCGCATTGCCGAGGATATTCGcc   >  8:112923/1‑149 (MQ=255)                                                                                                                                                  tA‑CCGTG‑‑AGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCTTTGTCCCTGGACGAAGCTTATCTGGATGTCAGCGACAGCCAGTGGTATTCGGGCAg                                >  8:89797/1‑119 (MQ=255)                                                                                                                                                  tA‑CCGTG‑‑AGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCTTTGTCCCTGGACGAAGCTTATCTGGATGTCAGCGACAGCCAGTGGTATTCGGGCAg                                <  7:89797/119‑1 (MQ=255)                                                                                                                                                  tA‑CCGTG‑‑AGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCTTTGTCCCTGGACGAAGCTTATCTGGATGTCAGCGACAGCCAGTGGTATTCGGGCAGCGCCACGCGCATTGCCGAGGATATTCGCCg  >  1:124032/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                     GGTGGGGGGCTCGCCCGACCATCGCGGGGTAATCGCCACCTGCAACTATGAAGCGCGTGCCTATGGTGTGCGCTCGGCCATGTCATCGCGGCATGCGCTGAAGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTA‑CCGTG‑‑AGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCTTTGTCCCTGGACGAAGCTTATCTGGATGTCAGCGACAGCCAGTGGTATTCGGGCAGCGCCACGCGCATTGCCGAGGATATTCGCCG  >  NC_002947/1381345‑1381638

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑