Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,381,496	G→T	6.2%	E83* (GAG→TAG)	dinB →	DNA polymerase IV

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,381,496	0	G	T	6.2%	206.2 / 3.6	64	E83* (GAG→TAG)	dinB	DNA polymerase IV
Reads supporting (aligned to +/- strand):  ref base G (32/28);  new base T (2/2);  total (34/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.28e-01

GGGGGCTCGCCCGACCATCGCGGGGTAATCGCCACCTGCAACTATGAAGCGCGTGCCTATGGTGTGCGCTCGGCCATGTCATCGCGGCATGCGCTGAAGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTA‑CCGTG‑‑AGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCTTTGTCCCTGGACGAAGCTTATCTGGATGTCAGCGACAGCCAGTGGTATTCGGGCAGCGCCACGCGCATTGCCGAGGATATTCGCCGGCGTG  >  NC_002947/1381349‑1381643                                                                                                                                                     |                                                                                                                                                      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                                                                                                                                          >  8:14366/1‑149 (MQ=255)           ccGACCATCGCGGGGTAATCGCCACCTGCAACTATGAAGCGCGTGCCTATGGTGTGCGCTCGGCCATGTCATCGCGGCATGCGCTGAAGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTA‑CCGTG‑‑AGGCCTCGCgg                                                                                                                                          <  4:83477/149‑1 (MQ=255)             gACCATCGCGGGGTAATCGCCACCTGCAACTATGAAGCGCGTGCCTATGGTGTGCGCTCGGCCATGTCATCGCGGCATGCGCTGAAGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTA‑CCGTG‑‑AGGCCTcg                                                                                                                                             <  3:37214/144‑1 (MQ=255)             gACCATCGCGGGGTAATCGCCACCTGCAACTATGAAGCGCGTGCCTATGGTGTGCGCTCGGCCATGTCATCGCGGCATGCGCTGAAGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTA‑CCGTG‑‑AGGCCTcg                                                                                                                                             >  4:37214/1‑144 (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                                                                                                                                    <  5:192644/98‑1 (MQ=255)                                                                    gcTCGGCCATGTCATCGCGGCATGCGCTGAAGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTA‑CCGTG‑‑AGGCCTCGCGGGAAATc                                                                                                                                    >  6:192644/1‑98 (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                                                                    >  8:108390/1‑130 (MQ=255)                                                                                                    cTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTA‑CCGTG‑‑AGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCTTTGTCCCTGGACGAAGCTTATCTgg                                                                    <  7:108390/130‑1 (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                                                                   >  4:180556/3‑98 (MQ=255)                                                                                                                                     cgcGAGGCCTC‑ACGGT‑‑AGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCTTTGTCCCTGGACGAAGCTTATCTGGa                                                                   <  3:180556/96‑1 (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)                                                                                                                                                   gTG‑‑ATGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCTTTGTCCCTGGACGAAGCTTATCTGGATGTCAGCGACAGCCAGTGGTATTCGGGCAGCGCCACGCGCATTGCCGAGGATATTCGCCGGCgt   <  6:133120/149‑1 (MQ=255)                                                                                                                                                    tG‑‑AGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCTTTGTCCCTGGACGAAGCTTATCTGGATGTCAGCGACAGCCAGTGGTATTCGGGCAGCGCCACGCGCATTGCCGAGGATATTCGCCGGCgtg  >  7:200854/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                      GGGGGCTCGCCCGACCATCGCGGGGTAATCGCCACCTGCAACTATGAAGCGCGTGCCTATGGTGTGCGCTCGGCCATGTCATCGCGGCATGCGCTGAAGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTA‑CCGTG‑‑AGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCTTTGTCCCTGGACGAAGCTTATCTGGATGTCAGCGACAGCCAGTGGTATTCGGGCAGCGCCACGCGCATTGCCGAGGATATTCGCCGGCGTG  >  NC_002947/1381349‑1381643

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