Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,381,495	T→G	6.7%	R82R (CGT→CGG)	dinB →	DNA polymerase IV

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,381,495	0	T	G	6.7%	218.5 / 4.5	71	R82R (CGT→CGG)	dinB	DNA polymerase IV
Reads supporting (aligned to +/- strand):  ref base T (32/34);  new base G (3/2);  total (35/36)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGGGCTCGCCCGACCATCGCGGGGTAATCGCCACCTGCAACTATGAAGCGCGTGCCTATGGTGTGCGCTCGGCCATGTCATCGCGGCATGCGCTGAAGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTACCGTGAGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCTTTGTCCCTGGAC‑‑GAAGCTTATCTGGATGTCAGCGACAGCCAGTGGTATTCGGGCAGCGCCACGCGCATTGCCGAGGATATTCGCCGGCG  >  NC_002947/1381349‑1381641                                                                                                                                                   |                                                                                                                                                     gggggCTCGCCCGACCATCGCGGGGTAATCGCCACCTGCAACTATGAAGCGCGTGCCTATGGTGTGCGCTCGGCCATGTCATCGCGGCATGCGCTGAAGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTACCGt                                                                                                                                                      >  6:166439/1‑147 (MQ=255) gggggCTCGCCCGACCATCGCGGGGTAATCGCCACCTGCAACTATGAAGCGCGTGCCTATGGTGTGCGCTCGGCCATGTCATCGCGGCATGCGCTGAAGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTACCGt                                                                                                                                                      >  4:109031/1‑147 (MQ=255)  ggggCTCGCCCGACCATCGCGGGGTAATCGCCACCTGCAACTATGAAGCGCGTGCCTATGGTGTGCGCTCGGCCATGTCATCGCGGCATGCGCTGAAGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTACCGTGAg                                                                                                                                                   <  3:367256/149‑1 (MQ=255)     gCTCGCCCGACCATCGCGGGGTAATCGCCACCTGCAACTATGAAGCGCGTGCCTATGGTGTGCGCTCGGCCATGTCATCGCGGCATGCGCTGAAGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTACCGTGAGGcc                                                                                                                                                <  2:360128/149‑1 (MQ=255)      cTCGCCCGACCATCGCGGGGTAATCGCCACCTGCAACTATGAAGCGCGTGCCTATGGTGTGCGCTCGGCCATGTCATCGCGGCATGCGCTGAAGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTCACGGTAGGCCt                                                                                                                                               >  5:114592/1‑149 (MQ=255)       tCGCCCGACCATCGCGGGGTAATCGCCACCTGCAACTATGAAGCACGTGCCTATGGTGTGCGCTCGGCCATGTCATCGCGGCATGCGCTGAAGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTACCGTGAGGCCTc                                                                                                                                              >  2:109738/1‑149 (MQ=255)           ccGACCATCGCGGGGTAATCGCCACCTGCAACTATGAAGCGCGTGCCTATGGTGTGCGCTCGGCCATGTCATCGCGGCATGCGCTGAAGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTCACGGTAGGCCTcgaag                                                                                                                                          <  6:114592/149‑4 (MQ=255)           ccGACCATCGCGGGGTAATCGCCACCTGCAACTATGAAGCGCGTGCCTATGGTGTGCGCTCGGCCATGTCATCGCGGCATGCGCTGAAGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTACCGTGAGGCCTCGCgg                                                                                                                                          >  2:65207/1‑149 (MQ=255)           ccGACCATCGCGGGGTAATCGCCACCTGCAACTATGAAGCGCGTGCCTATGGTGTGCGCTCGGCCATGTCATCGCGGCATGCGCTGAAGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTACCGTGAGGCCTCGCgg                                                                                                                                          >  2:314672/1‑149 (MQ=255)            cGACCATCGCGGGGTAATCGCCACCTGCAACTATGAAGCGCGTGCCTATGGTGTGCGCTCGGCCATGTCATCGCGGCATGCGCTGAAGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTCACGGTAGGCCTcgaagc                                                                                                                                         <  1:327421/149‑5 (MQ=255)               ccATCGCGGGGTAATCGCCACCTGCAACTATGAAGCGCGTGCCTATGGTGTGCGCTCGGCCATGTCATCGCGGCATGCGCTGAAGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTACCGTGAGGCCTCGCGGGaaa                                                                                                                                      >  5:261927/1‑149 (MQ=255)               ccATCGCGGGGTAATCGCCACCTGCAACTATGAAGCGCGTGCCTATGGTGTGCGCTCGGCCATGTCATCGCGGCATGCGCTGAAGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTACCGTGAGGCCTCGCGGGaaa                                                                                                                                      >  8:203809/1‑149 (MQ=255)                    gcgGGGTAATCGCCACCTGCAACTATGAAGCGCGTGCCTATGGTGTGCGCTCGGCCATGTCATCGCGGCATGCGCTGAAGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTACCGTGAGGCCTCGCGGGAAATCcac                                                                                                                                 <  3:207533/149‑1 (MQ=255)                          tAATCGCCACCTGCAACTATGAAGCGCGTGCCTATGGTGTGCGCTCGGCCATGTCATCGCGGCATGCGCTGAAGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTACCGTGAGGCCTCGCGGGAAATCCACACGATc                                                                                                                           <  8:273359/149‑1 (MQ=255)                             tCGCCACCTGCAACTATGAAGCGCGTGCCTATGGTGTGCGCTCGGCCATGTCATCGCGGCATGCGCTGAAGCTTTGCCCCGAGCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTACCGTGAGGCCTCGCGGGAAATCCACACGATCttt                                                                                                                        >  6:314970/1‑149 (MQ=255)                              cGCCACCTGCAACTATGAAGCGCGTGCCTATGGTGTGCGCTCGGCCATGTCATCGCGGCATGCGCTGAAGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTACCGTGAGGCCTCGCGGGAAATCCACACGATCTTTc                                                                                                                       >  8:201441/1‑149 (MQ=255)                               gCCACCTGCAACTATGAAGCGCGTGCCTATGGTGTGCGCTCGGCCATGTCATCGCGGCATGCGCTGAAGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTACCGTGAGGCCTCGCGGGAAATCCACACGATCTTTcg                                                                                                                      <  5:200774/149‑1 (MQ=255)                               gCCACCTGCAACTATGAAGCGCGTGCCTATGGTGTGCGCTCGGCCATGTCATCGCGGCATGCGCTGAAGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTACCGTGAGGCCTCGCGGGAAATCCACACGATCTTTc                                                                                                                       <  5:341127/148‑1 (MQ=255)                                     tGCAACTATGAAGCGCGTGCCTATGGTGTGCGCTCGGCCATGTCATCGCGGCATGCGCTGAAGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTACCGTGAGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTa                                                                                                                >  6:384257/1‑149 (MQ=255)                                          cTATGAAGCGCGTGCCTATGGTGTGCGCTCGGCCATGTCATCGCGGCATGCGCTGAAGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTACCGTGAGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACg                                                                                                            >  8:232449/1‑148 (MQ=255)                                           tATGAAGCGCGTGCCTATGGTGTGCGCTCGGCCATGTCATCGCGGCATGCGCTGAAGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTACCGTGAGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACg                                                                                                            >  4:72722/1‑147 (MQ=255)                                           tATGAAGCGCGTGCCTATGGTGTGCGCTCGGCCATGTCATCGCGGCATGCGCTGAAGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTACCGTGAGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACg                                                                                                            >  4:168331/1‑147 (MQ=255)                                             tGAAGCGCGTGCCTATGGTGTGCGCTCGGCCATGTCATCGCGGCATGCGCTGAAGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTACCGTGAGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGc                                                                                                        >  5:165743/1‑149 (MQ=255)                                                  cgcgTGCCTATGGTGTGCGCTCGGCCATGTCATCGCGGCATGCGCTGAAGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTACCGTGAGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATc                                                                                                   <  7:203809/149‑1 (MQ=255)                                                     gTGCCTATGGTGTGCGCTCGGCCATGTCATCGCGGCATGCGCTGAAGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTACCGTGAGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAg                                                                                                >  8:266743/1‑149 (MQ=255)                                                         cTATGGTGTGCGCTCGGCCATGTCATCGCGGCATGCGCTGAAGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTACCGTGAGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCtt                                                                                            >  1:69090/1‑149 (MQ=255)                                                          tATGGTGTGCGCTCGGCCATGTCATCGCGGCATGCGCTGAAGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTACCGTGAGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCttt                                                                                           >  5:217732/1‑149 (MQ=255)                                                          tATGGTGTGCGCTCGGCCATGTCATCGCGGCATGCGCTGAAGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTACCGTGAGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCttt                                                                                           >  7:79416/1‑149 (MQ=255)                                                           aTGGTGTGCGCTCGGCCATGTCATCGCGGCATGCGCTGAAGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTACCGTGAGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCTTTg                                                                                          >  7:261097/1‑149 (MQ=255)                                                            tGGTGTGCGCTCGGCCATGTCATCGCGGCATGCGCTGAAGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTACCGTGAGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCttt                                                                                           <  1:65207/147‑1 (MQ=255)                                                                 tgCGCTCGGCCATGTCATCGCGGCATGCGCTGAAGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTACCGTGAGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCTTTGTCCCt                                                                                     <  6:217732/148‑1 (MQ=255)                                                                   cgcTCGGCCATGTCATCGCGGCATGCGCTGAAGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTACCGTGAGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCTTTGTCCCTGGa                                                                                  >  6:319599/1‑149 (MQ=255)                                                                   cgcTCGGCCATGTCATCGCGGCATGCGCTGAAGCTTTGCCCCGAGCTGTTGATCGTCAAGCCGCTCTTCGAGGCCTACCGTGAGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCTTTGTCCCTGGa                                                                                  <  3:59418/149‑1 (MQ=255)                                                                    gcTCGGCCATGTCATCGCGGCATGCGCTGAAGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTACCGTGAGGCCTCGCGGGAAATCCACACGATCt                                                                                                                          <  4:125584/108‑1 (MQ=255)                                                                    gcTCGGCCATGTCATCGCGGCATGCGCTGAAGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTACCGTGAGGCCTCGCGGGAAATCCACACGATCt                                                                                                                          >  3:125584/1‑108 (MQ=255)                                                                        ggCCATGTCATCGCGGCATGCGCTGAAGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTACCGTGAGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCTTTGTCCCTGGAC‑‑GAAg                                                                           <  7:327464/149‑1 (MQ=255)                                                                        ggCCATGTCATCGCGGCATGCGCTGAAGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTACCGTGAGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCTTTGTCCCTGGAC‑‑GAAg                                                                           <  1:109738/149‑1 (MQ=255)                                                                           cATGTCATCGCGGCATGCGCTGAAGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTACCGTGAGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCTTTGTCCCTGGAC‑‑GAAGCtt                                                                        <  3:168331/149‑1 (MQ=255)                                                                                 aTCGCGGCATGCGCTGAAGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTACCGTGAGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCTTTGTCCCTGGAC‑‑GAAGCTTATCTg                                                                   <  3:153382/148‑1 (MQ=255)                                                                                 aTCGCGGCATGCGCTGAAGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTACCGTGAGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCTTTGTCCCTGGAC‑‑GAAGCTTATCTg                                                                   >  5:33036/1‑148 (MQ=255)                                                                                 aTCGCGGCATGCGCTGAAGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTACCGTGAGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCTTTGTCCCTGGAC‑‑GAAGCTTATCTg                                                                   >  7:41759/1‑148 (MQ=255)                                                                                  tCGCGGCATGCGCTGAAGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTACCGTGAGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCTTTGTCCCTGGAC‑‑GAAGCTTATCTGGa                                                                 <  5:296307/149‑1 (MQ=255)                                                                                  tCGCGGCATGCGCTGAAGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTACCGTGAGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCTTTGTCCCTGGAC‑‑GAAGCTTATCTGGa                                                                 >  4:63198/1‑149 (MQ=255)                                                                                   cgcgGCATGCGCTGAAGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTACCGTGAGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCTTTGTCCCTGGAC‑‑GAAGCTTATCTGGAt                                                                <  4:359464/149‑1 (MQ=255)                                                                                    gcgGCATGCGCTGAAGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTACCGTGAGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCTTTGTCCCTGGAC‑‑GAAGCTTATCTGGAt                                                                <  3:109031/148‑1 (MQ=255)                                                                                         aTGCGCTGAAGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTACCGTGAGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCTTTGTCCCTGGAC‑‑GAAGCTTATCTGGATGTCAGc                                                          <  4:77723/149‑1 (MQ=255)                                                                                          tGCGCTGAAGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTACCGTGAGGCCTCGCGGGGAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCTTTGTCCCTGGAC‑‑GAAGCTTATCTGGATGTCAGc                                                          >  4:10460/1‑148 (MQ=255)                                                                                             gcTGAAGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTACCGTGAGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCTTTGTCCCTGGAC‑‑GAAGCTTATCTGGATGTCAGCGACa                                                      <  4:316582/149‑1 (MQ=255)                                                                                             gcTGAAGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTACCGTGAGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCTTTGTCCCTGGAC‑‑GAAGCTTATCTGGATGTCAGCGACa                                                      <  1:314672/149‑1 (MQ=255)                                                                                                 aaGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTACCGTGAGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCTTTGTCCCTGGAC‑‑GAAGCTTATCTGGATGTCAGCGACAGcc                                                   <  7:266743/148‑1 (MQ=255)                                                                                                 aaGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTACCGTGAGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCTTTGTCCCTGGAC‑‑GAAGCTTATCTGGATGTCAGCGACAGc                                                    >  3:200031/1‑147 (MQ=255)                                                                                                    cTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTACCGTGAGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCTTTGTCCCTGGAC‑‑GAAGCTTATCTGGATGTCAGCGACAGCCAGTg                                               <  3:72722/149‑1 (MQ=255)                                                                                                    cTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTACCGTGAGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCTTTGTCCCTGGAC‑‑GAAGCTTATCTGGATGTCAGCGACAGCCAGTg                                               <  6:261927/149‑1 (MQ=255)                                                                                                     tttGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTACCGTGAGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCttt                                                                                           <  5:344967/106‑1 (MQ=255)                                                                                                     tttGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTACCGTGAGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCttt                                                                                           >  6:344967/1‑106 (MQ=255)                                                                                                        gCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTACCGTGAGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCTTTGTCCCTGGAC‑‑GAAGCTTATCTGGATGTCAGCGACAGCCAGTGGTa                                            <  8:79416/148‑1 (MQ=255)                                                                                                           ccGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTACCGTGAGGCCTCGCGGGAAATCCCCACGATCTTTCGCGATTCCACGCAGCTCATCGAGCCTTTGTCCCTGGAC‑‑GAAGCTTATCTGGATGTCAGCGACAGCCAGTGGTATTCg                                        <  6:165743/149‑1 (MQ=255)                                                                                                              aTCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTACCGTGAGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCTTTGTCCCTGGAC‑‑GAAGCTTATCTGGATGTCAGCGACAGCCAGTGGTATTCggg                                      >  1:348804/1‑148 (MQ=255)                                                                                                                  gTTGATCGTCAAGCCGCGCTTCGAGGCCTACCGTGAGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCTTTGTCCCTGGAC‑‑GAAGCTTATCTGGATGTCAGCGACAGCCAGTGGTATTCGGGCAgcg                                 >  2:191298/1‑149 (MQ=255)                                                                                                                  gTTGATCGTCAAGCCGCGCTTCGAGGCCTACCGTGAGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCTTTGTCCCTGAACCTTGAGCTTATCTGGATGTCAGCGACAGCCAGTGGTATTCGGGCAg                                   >  1:126017/1‑149 (MQ=255)                                                                                                                      aTCGTCAAGCCGCGCTTCGAGGCCTACCGTGAGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCTTTGTCCCTGGAC‑‑GAAGCTTATCTGGATGTCAGCGACAGCCAGTGGTATTCGGGCAGCGCCAc                             <  3:63198/149‑1 (MQ=255)                                                                                                                             aGCCGCGCTTCGAGGCCTACCGTGAGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCTTTGTCCCTGGAC‑‑GAAGCTTATCTGGATGTCAGCGACAGCCAGTGGTATTCGGGCAGCGCCACGCGCAtt                      <  5:384257/149‑1 (MQ=255)                                                                                                                              gCCGCGCTTCGAGGCCTACCGTGAGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCAGCGAGCCTTTGTCCCTGGAC‑‑GAAGCTTATCTGGATGTCAGCGACAGCCAGTGGTATTCGGGCAGCGCCACGCGCATTg                     <  5:319599/149‑1 (MQ=255)                                                                                                                                cgcgCTTCGAGGCCTACCGTGAGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCTTTGTCCCTGGAC‑‑GAAGCTTATCTGGATGTCAGCGACAGCCAGTGGTATTCGGGCAGCGCCACGCGCATTGcc                   >  1:62799/1‑149 (MQ=255)                                                                                                                                  cgcTTCGAGGCCTACCGTGAGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCTTTGTCCCTGGAC‑‑GAAGCTTATCTGGATGTCAGCGACAGCCAGTGGTATTCGGGCAGCGCCACGCGCATTGCCGa                 <  5:166439/149‑1 (MQ=255)                                                                                                                                   gcTTCGAGGCCTACCGTGAGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCTTTGTCCCTGGAC‑‑GAAGCTTATCTGGATGTCAGCGACAGCCAGTGGTATTCGGGCAGCGCCACGCGCATTGCCGAg                <  5:314970/149‑1 (MQ=255)                                                                                                                                    ccgcGAGGCCTCACGGTAGGCCTCGCGGGAAATCCCCACGATCTTTCGCGATTACACGGAGCTCATCGAGCCTTTGTCCCTGGAC‑‑GAAGCTTATCTGGATGTCAGCGACAGCCAGTGGTATTCGGGCAGCGCCACGCGCATTGCCGAg                >  6:268517/4‑148 (MQ=255)                                                                                                                                    ccgcGAGGCCTCACGGTAGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCTTTGTCCCTGGAC‑‑GAAGCTTATCTGGATGTCAGCGACAGCCAGTGGTATTCGGGCAGCGCCACGCGCATTGCCGAg                >  2:20382/4‑148 (MQ=255)                                                                                                                                    cTTCGAGGCCTACCGTGAGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCTTTGTCCCTGGAC‑‑GAAGCTTATCTGGATGTCAGCGACAGCCAGTGGTATTCGGGCAGCGCCACGCGCATTGCCGAg                >  8:355035/1‑148 (MQ=255)                                                                                                                                         aGGCCTACCGTGAGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCTTTGTCCCTGGAC‑‑GAAGCTTATCTGGATGTCAGCGACAGCCAGTGGTATTCgg                                       >  4:243199/1‑120 (MQ=255)                                                                                                                                         aGGCCTACCGTGAGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCTTTGTCCCTGGAC‑‑GAAGCTTATCTGGATGTCAGCGACAGCCAGTGGTATTCgg                                       <  3:243199/120‑1 (MQ=255)                                                                                                                                         aGGCCTACCGTGAGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCTTTGTCCCTGGAC‑‑GAAGCTTATCTGGATGTCAGCGACAGCCAGTGGTATTCGGGCAGCGCCACGCGCATTGCCGAGGATAtt          >  1:166669/1‑149 (MQ=255)                                                                                                                                                 cGTGAGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCTTTGTCCCTGGAC‑‑GAAGCTTATCTGGATGTCAGCGACAGCCAGTGGTATTCGGGCAGCGCCACGCGCATTGCCGAGGATATTCGCCGGCg  <  2:166669/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                     GGGGGCTCGCCCGACCATCGCGGGGTAATCGCCACCTGCAACTATGAAGCGCGTGCCTATGGTGTGCGCTCGGCCATGTCATCGCGGCATGCGCTGAAGCTTTGCCCCGATCTGTTGATCGTCAAGCCGCGCTTCGAGGCCTACCGTGAGGCCTCGCGGGAAATCCACACGATCTTTCGCGATTACACGGAGCTCATCGAGCCTTTGTCCCTGGAC‑‑GAAGCTTATCTGGATGTCAGCGACAGCCAGTGGTATTCGGGCAGCGCCACGCGCATTGCCGAGGATATTCGCCGGCG  >  NC_002947/1381349‑1381641

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