Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,738,440	T→A	7.4%	V251V (GTA→GTT)	PP_1532 ←	putative Phage integrase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,738,440	0	T	A	7.4%	201.4 / 4.5	67	V251V (GTA→GTT)	PP_1532	putative Phage integrase
Reads supporting (aligned to +/- strand):  ref base T (32/30);  new base A (2/3);  total (34/33)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

TTTGCTGATTGTGTCGAGCATGATGGGAGATCCCGCCTTCGCGCGGCTGATTGTCCGGTGGTGGTGCACCAGGTGCTTTGACAGCACCCGATCGCGGCAACGCCGT‑ATCACTTCACCCAAGCTCAGGCTGAT‑GCATTCCAGCGTTACGCTGGTGCTGATCCTTAGCCTGGCCCCCGTCTTCGATTGGACGACCTGCAGGAATCCGTCCACCTCATCCCGGAAAACGAATGAAGCAATGTCCTCGCGGCGCTGCCCGGTTATAAGAGCGAGCTCCATGGCGCGCTT  >  NC_002947/1738296‑1738580                                                                                                                                                   |                                                                                                                                             tttGCTGATTGTGTCGAGCATGATGGGAGATCCCGCCTTCGCGCGGCTGATTGTCCGGTGGTGGTGCACCAGGTGCTTTGACAGCACCCGATCGCGGCAACGCCGT‑ATCACTTCACCCAAGCTCAGGCTGAT‑GCATTCCAGCGTAACGc                                                                                                                                          >  7:227998/1‑149 (MQ=255)    gCTGATTGTGTCGAGCATGATGGGAGATCCCGCCTTCGCGCGGCTGATTGTCCGGTGGTGGTGCACCAGGTGCTTTGACAGCACCCGATCGCGGCAACGCCGT‑ATCACTTCACCCAAGCTCAGGCTGAT‑GCATTCCAGCGTAACGCTg                                                                                                                                        <  3:220457/148‑1 (MQ=255)     cTGATTGTGTCGAGCATGATGGGAGATCCCGCCTTCGCGCGGCTGATTGTCCGGTGGTGGTGCACCAGGTGCTTTGACAGCACCCGATCGCGGCAACGCCGT‑ATCACTTCACCCAAGCTCAGGCTGAT‑GCATTCCAGCGTTACGCTgg                                                                                                                                       <  7:170571/148‑1 (MQ=255)        aTTGTGTCGAGCATGATGGGAGATCCCGCCTTCGCGCGGCTGATTGTCCGGTGGTGGTGCACCAGGTGCTTTGACAGCACCCGATCGCGGCAACGCCGT‑ATCACTTCACCCAAGCTCAGGCTGAT‑GCATTCCAGCGTAACGCTggaat                                                                                                                                    <  8:227998/148‑4 (MQ=255)        aTTGTGTCGAGCATGATGGGAGATCCCGCCTTCGCGCGGCTGATTGTCCGGTGGTGGTGCACCAGGTGCTTTGACAGCACCCGATCGCGGCAACGCCGT‑ATCACTTCACCCAAGCTCAGGCTGAT‑GCATTCCAGCGTAACGCTggaat                                                                                                                                    <  1:99262/148‑4 (MQ=255)                gAGCATGATGGGAGATCCCGCCTTCGCGCGGCTGATTGTCCGGTGGTGGTGCACCAGGTGCTTTGACAGCACCCGATCGCGGCAACGCCGT‑ATCACTTCACCCAAGCTCAGGCTGAT‑GCATTCCAGCGTTACGCTGGTGCTGATCCt                                                                                                                             >  5:51327/1‑147 (MQ=255)                gAGCATGATGGGAGATCCCGCCTTCGCGCGGCTGATTGTCCGGTGGTGGTGCACCAGGTGCTTTGACAGCACCCGATCGCGGCAACGCCGT‑ATCACTTCACCCAAGCTCAGGCTGAT‑GCATTCCAGCGTTACGCTGGTGCTGATCCTTa                                                                                                                           <  8:98188/149‑1 (MQ=255)                gAGCATGATGGGAGATCCCGCCTTCGCGCGGCTGATTGTCCGGTGGTGGTGCACCAGGTGCTTTGACAGCACCCGATCGCGGCAACGCCGT‑ATCACTTCACCCAAGCTCAGGCTGAT‑GCATTCCAGCGTTACGCTGGTGCTGATCCTTa                                                                                                                           <  6:51327/149‑1 (MQ=255)                   cATGATGGGAGATCCCGCCTTCGCGCGGCTGATTGTCCGGTGGTGGTGCACCAGGTGCTTTGACAGCACCCGATCGCGGCAACGCCGT‑ATCACTTCACCCAAGCTCAGGCTGAT‑GCATTCCAGCGTTACGCTGGTGCTGATCCTTAGcc                                                                                                                        <  7:156190/149‑1 (MQ=255)                      gatgGGAGATCCCGCCTTCGCGCGGCTGATTGTCCGGTGGTGGTGCACCAGGTGCTTTGACAGCACCCGATCGCGGCAACGCCGT‑ATCACTTCACCCAAGCTCAGGCTGAT‑GCATTCCAGCGTTACGCTGGTGCTGATCCTTAGCCTgg                                                                                                                     >  1:78260/1‑149 (MQ=255)                      gatgGGAGATCCCGCCTTCGCGCGGCTGATTGTCCGGTGGTGGTGCACCAGGTGCTTTGACAGCACCCGATCGCGGCAACGCCGT‑ATCACTTCACCCAAGCTCAGGCTGAT‑GCATTCCAGCGTTACGCTGGTGCTGATCCTTAGCCTgg                                                                                                                     >  5:377058/1‑149 (MQ=255)                       atgGGAGATCCCGCCTTCGCGCAGCTGATTGTCCGGTGGTGGGGCACCAGGTGCTTTGACAGCACCCGATCGCGGCAACGCCGT‑ATCACTTCACCCAAGCTCAGGCTGAT‑GCATTCCAGCGTTACGCTGGTGCTGATCCTTAGCCTGGc                                                                                                                    <  5:57740/149‑1 (MQ=255)                        tgGGAGATCCCGCCTTCGCGCGGCTGATTGTCCGGTGGTGGTGCACCAGGTGCTTTGACAGCACCCGATCGCGGCAACGCCGT‑ATCACTTCACCCAAGCTCAGGCTGAT‑GCATTCCAGCGTTACGCTGGTGCTGATCCTTAGCCTGGcc                                                                                                                   >  4:187619/1‑149 (MQ=255)                           gagaTCCCGCCTTCGCGCGGCTGATTGTCCGGTGGTGGTGCACCAGGTGCTTTGACAGCACCCGATCGCGGCAACGCCGT‑ATCACTTCACCCAAGCTCAGGCTGAT‑GCATTCCAGCGTTACGCTGGTGCTGATCCTTAGCCTGGccccc                                                                                                                >  2:289735/1‑149 (MQ=255)                                 ccGCCTTCGCGCGGCTGATTGTCCGGTGGTGGTGCACCAGGTGCTTTGACAGCACCCGATCGCGGCAACGCCGT‑ATCACTTCACCCAAGCTCAGGCTGAT‑GCATTCCAGCGTTACGCTGGTGCTGATCCTTAGCCTGGCCCCCGTCTTc                                                                                                          <  6:377058/149‑1 (MQ=255)                                   gCCTTCGCGCGGCTGATTGTCCGGTGGTGGTGCACCAGGTGCTTTGACAGCACCCGATCGCGGCAACGCCGT‑ATCACTTCACCCAAGCTCAGGCTGAT‑GCATTCCAGCGTTACGCTGGTGCTGATCCTTAGCCTGGCCCCCGTCTTCGa                                                                                                        >  6:185295/1‑149 (MQ=255)                                   gCCTTCGCGCGGCTGATTGTCCGGTGGTGGTGCACCAGGTGCTTTGACAGCACCCGATCGCGGCAACGCCGT‑ATCACTTCACCCAAGCTCAGGCTGAT‑GCATTCCAGCGTTACGCTGGTGCTGATCCTTAGCCTGGCCCCCGTCTTCGa                                                                                                        >  4:289701/1‑149 (MQ=255)                                       tCGCGCGGCTGATTGTCCGGTGGTGGTGCACCAGGTGCTTTGACAGCACCCGATCGCGGCAACGCCGT‑ATCACTTCACCCAAGCTCAGGCTGAT‑GAATTCCATCGTTACGCTGGTGCTGATCCTTAGACTGGCACCCGTCTTCGATTg                                                                                                     >  1:378655/1‑148 (MQ=255)                                        cgcgcgGCTGATTGTCCGGTGGTGGTGCACCAGGTGCTTTGACAGCACCCGATCGCGGCAACGCCGT‑ATCACTTCACCCAAGCTCAGGCTGAT‑GCATTCCAGCGTTACGCTGGTGCTGATCCTTAGCCTGGCCCCCGTCTTCGATTgg                                                                                                    <  2:4585/148‑1 (MQ=255)                                        cgcgcgGCTGATTGTCCGGTGGTGGTGCACCAGGTGCTTTGACAGCACCCGATCGCGGCAACGCCGT‑ATCACTTCACCCAAGCTCAGGCTGAT‑GCATTCCAGCGTTACGCTGGTGCTGATCCTTAGCCTGGCCCCCGTCTTCGATTgg                                                                                                    <  5:230968/148‑1 (MQ=255)                                          cgcgGCTGATTGTCCGGTGGTGGTGCACCAGGTGCTTTGACAGCACCCGATCGCGGCAACGCCGT‑ATCACTTCACCCAAGCTCAGGCTGAT‑GCATTCCAGCGTTACGCTGGTGCTGATCCTTAGCCTGGCCCCCGTCt                                                                                                            >  2:240637/1‑138 (MQ=255)                                          cgcgGCTGATTGTCCGGTGGTGGTGCACCAGGTGCTTTGACAGCACCCGATCGCGGCAACGCCGT‑ATCACTTCACCCAAGCTCAGGCTGAT‑GCATTCCAGCGTTACGCTGGTGCTGATCCTTAGCCTGGCCCCCGTCt                                                                                                            <  1:240637/138‑1 (MQ=255)                                            cgGCTGATTGTCCGGTGGTGGTGCACCAGGTGCTTTGACAGCACCCGATCGCGGCAACGCCGT‑ATCACTTCACCCAAGCTCAGGCTGAT‑GCATTCCAGCGTTACGCTGGTGCTGATCCTTAGCCTGGCCCCCGTCTTCGATTGgacgac                                                                                               <  1:289735/149‑1 (MQ=255)                                               cTGATTGTCCGGTGGTGGTGCACCAGGTGCTTTGACAGCACCCGATCGCGGCAACGCCGT‑ATCACTTCACCCAAGCTCAGGCTGAT‑GCATTCCAGCGTTACGCTGGTGCTGATCCTTAGCCTGGCCCCCGTCTTCGATTGGACGACCTg                                                                                            <  2:126363/149‑1 (MQ=255)                                               cTGATTGTCCGGTGGTGGTGCACCAGGTGCTTTGACAGCACCCGATCGCGGCAACGCCGT‑ATCACTTCACCCAAGCTCAGGCTGAT‑GCATTCCAGCGTTACGCTGGTGCTGATCCTTAGCCTGGCCCCCGTCTTCGATTGGACGACCTg                                                                                            <  7:38395/149‑1 (MQ=255)                                                        cGGTGGTGGTGCACCAGGTGCTTTGACAGCACCCGATCGCGGCAACGCCGT‑ATAACTTCACCCAAGCTCAGGCTGAT‑GCATTCCAGCGTTACGCTGGTGCTGATCCTTAGCCTGGCCCCCGGCTTCGATTGGACGCCCTGCAGGAATcc                                                                                   >  1:220609/1‑149 (MQ=255)                                                          gtggtggtgCACCAGGTGCTTTGACAGCACCCGATCGCGGCAACGCCGT‑ATCACTTCACCCAAGCTCAGGCTGAT‑GCATTCCAGCGTTACGCTGGTGCTGATCCTTAGCCTGGCCCCCGTCTTCGATTGGACGACCTGCAGGAATCCGt                                                                                 >  3:193669/1‑149 (MQ=255)                                                               ggtgCACCAGGTGCTTTGACAGCACCCGATCGCGGCAACGCCGT‑ATCACTTCACCCAAGCTCAGGCTGAT‑GCATTCCAGCGTTACGCTGGTGCTGATCCTTAGCCTGGCCCCCGTCTTCGATTGGACGACCTGCAGGAATCCGTCCAcc                                                                            >  7:357586/1‑149 (MQ=255)                                                                gtgCACCAGGTGCTTTGACAGCACCCGATCGCGGCAACGCCGTAATCACTTCACCCAAGCTCAGGCTGAT‑GCATTCCAGCGTTACGCTGGTGCTGATCCTTAGCCTGGCCCCCGTCTTCGATTGGACGACCTGCAGGAATCCGTCCAcc                                                                            >  1:270637/1‑149 (MQ=255)                                                                gtgCACCAGGTGCTTTGACAGCACCCGATCGCGGCAACGCCGT‑ATCACTTCACCCAAGCTCAGGCTGAT‑GCATTCCAGCGTTACGCTGGTGCTGATCCTTAGCCTGGCCCCCGTCTTCGATTGgacgac                                                                                               >  7:335972/1‑129 (MQ=255)                                                                gtgCACCAGGTGCTTTGACAGCACCCGATCGCGGCAACGCCGT‑ATCACTTCACCCAAGCTCAGGCTGAT‑GCATTCCAGCGTTACGCTGGTGCTGATCCTTAGCCTGGCCCCCGTCTTCGATTGgacgac                                                                                               <  8:335972/129‑1 (MQ=255)                                                                   gatcgggTGCTGTCAAAGCACCCGATCGCGGCAACGCCGT‑ATCACTTCACCCAAGCTCAGGCTGAT‑GCATTCCAGCGTTACGCTGGTGCTGATCCTTAGCCTGGCCCCCGTCTTCGATTGGACGACCTGCAGGAATCCGTCCACCTCAt                                                                        >  3:62772/6‑149 (MQ=255)                                                                      cAGGTGCTTTGACAGCACCCGATCGCGGCAACGCCGT‑ATCACTTCACCCAAGCTCAGGCTGAT‑GCATTCCAGCGTTACGCTGGTGCTGATCCTTAGCCTGGCCCCCGTCTTCGATTGGACGACCTGCAGGAATCCGTCCACCTCATccc                                                                     >  7:179548/1‑149 (MQ=255)                                                                             tttGACAGCCCCCGATCGCGGCAACGCCGT‑ATCACTTCACCCAAGCTCAGGCTGAT‑GCATTCCAGCGTTACGCTGGTGCTGATCCTTAGCCTGGCCCCCGTCTTAGATTGGACGACCTGCAGGAATCCGTCCACCTCATCCCGGAAAAc                                                              >  5:105976/1‑149 (MQ=255)                                                                               tGACAGCACCCGATCGCGGCAACGCCGT‑ATCACTTCACCCAAGCTCAGGCTGAT‑GCATTCCAGCGTTACGCTGGTGCTGATCCTTAGCCTGGCCCCCGTCTTCGATTGGACGACCTGCAGGAATCCGTCCACCTCATCCCGGAAAAc                                                              <  7:27214/147‑1 (MQ=255)                                                                               tGACAGCACCCGAGCGCGGCAACGCCGT‑ATCACTTCACCCAAGCTCAGGCTGAT‑GCATTCCAGCGTTACGCTGGTGCTGATCCTTAGCCTGGCCCCCGTCTTCGATTGGACGACCTGCAGGAATCCGTCCACCTCATCCCGGAAAAc                                                              <  1:298022/147‑1 (MQ=255)                                                                                 aCAGCACCCGATCGCGGCAACGCCGT‑ATCACTTCACCCAAGCTCAGGCTGAT‑GCATTCCAGCGTTACGCTGGTGCTGATCCTTAGCCTGGCCCCCGTCTTCGATTGGACGACCTGCAGGAATCCGTCCACCTCATCCCGGAAAACGaa                                                           <  3:187619/148‑1 (MQ=255)                                                                                 aCAGCACCCGATCGCGGCAACGCCGT‑ATCACTTCACCCAAGCTCAGGCTGAT‑GCATTCCAGCGTTACGCTGGTGCTGATCCTTAGCCTGGCCCCCGTCTTCGATTGGACGACCTGCAGGAATCCGTCCACCTCATCCCGGAAAACGAAt                                                          <  5:185295/149‑1 (MQ=255)                                                                                    gCACCCGATCGCGGCAACGCCGT‑ATCACTTCACCCAAGCTCAGGCTGAT‑GCATTCCAGCGTTACGCTGGTGCTGATCCTTAGCCTGGCCCCCGTCTTCGATTGGACGACCTGCAGGAATCCGTCCACCTCATCCCGGAAAACGAATGaa                                                       <  8:179548/149‑1 (MQ=255)                                                                                    gCACCCGATCGCGGCAACGCCGT‑ATCACTTCACCCAAGCTCAGGCTGAT‑GCATTCCAGCGTTACGCTGGTGCTGATCCTTAGCCTGGCCCCCGTCTTCGATTGGACGACCTGCAGGAATCCGTCCACCTCATCCCGGAAAACGAATGaa                                                       <  7:214006/149‑1 (MQ=255)                                                                                      aCCCGATCGCGGCAACGCCGT‑ATCACTTCACCCAAGCTCAGGCTGAT‑GCATTCCAGCGTTACGCTGGTGCTGATCCTTAGCCTGGCCCCCGTCTTCGATTGGACGACCTGCAGGAATCCGTCCCCCTCATCCCGGAAAACGAATGAAg                                                      >  8:176661/1‑148 (MQ=255)                                                                                      aCCCGATCGCGGCAACGCCGT‑ATCACTTCACCCAAGCTCAGGCTGAT‑GCATTCCAGCGTTACGCTGGTGCTGATCCTTAGCCTGGCCCCCGTCTTCGATTGGACGACCTGCAGGAATCCGTCCACCTCATCCCGGAAAACGAATGAAg                                                      >  3:297536/1‑148 (MQ=255)                                                                                        ccGATCGCGGCAACGCGGT‑ATCACTTCACCCAAGCTCAGGCTGATGGCATTCCAGCGTTACGCTGGTGCTGATCCTTAGCCTGGCCCCCGTCTTCGATTGGACGACCTGCAGGAATCCGTCCACCTCATACCGGAAAACGAATGAAGc                                                     >  1:276067/1‑148 (MQ=255)                                                                                        ccGATCGCGGCAACGCCGT‑ATCACTTCACCCAAGCTCAGGCTGAT‑GCATTCCAGCGTTACGCTGGTGCTGATCCTTAGCCTGGCCCCCGTCTTCTATTTGACGACCTGCAGGAATCCGTCCACCTCATCCCGGAAAACGAATGAAGCaa                                                   >  5:286252/1‑149 (MQ=255)                                                                                        ccGATCGCGGCAACGCCGT‑ATCACTTCACCCAAGCTCAGGCTGAT‑GCATTCCAGCGTTACGCTGGTGCTGATCCTTAGCCTGGCCCCCGTCTTCGATTGGACGACCTGCAGGAATCCGTCCACCTCATCCCGGAAAACGAATGAAGCaa                                                   >  1:231767/1‑149 (MQ=255)                                                                                         cGATCGCGGCAACGCCGT‑ATCACTTCACCCAAGCTCAGGCTGAT‑GCATTCCAGCGTTACGCTGGTGCTGATCCTTAGCCTGGCCCCCGTCTTCGATTGGACGACCTGCAGGAATCCGTCCACCTCATCCCGGAAAACGAATGAAGCAAt                                                  <  4:193669/149‑1 (MQ=255)                                                                                               cgGCAACGCCGT‑ATCACTTCACCCAAGCTCAGGCTGAT‑GCATTCCAGCGTTACGCTGGTGCTGATCCTTAGCCTGGCCCCCGTCTTCGATTGGACGACCTGCAGGAATCCGTCCACCTCATCCCGGAAAACGAATGAAGCAATGTCCTc                                            >  7:262298/1‑149 (MQ=255)                                                                                               cgGCAACGCCGT‑ATCACTTCACCCAAGCTCAGGCTGAT‑GCATTCCAGCGTTACGCTGGTGCTGATCCTTAGCCTGGCCCCCGTCTTCGATTGGACGACCTGCAGGAATCCGTCCACCTCATCCCGGAAAACGAATGAAGCAATGTCCTc                                            >  8:20564/1‑149 (MQ=255)                                                                                                ggCAACGCCGT‑ATCACTTCACCCAAGCTCAGGCTGAT‑GCATTCCAGCGTTACGCTGGTGCTGATCCTTAGCCTGGCCCCCGTCTTCGATTGGACGACCTGCAGGAATCCGTCCACCTCATCCCGGAAAACGAATGAAGCAATGTCCTcg                                           <  4:297536/149‑1 (MQ=255)                                                                                                                   cACCCAAGCTCAGGCTGAT‑GCATTCCAGCGTTACGCTGGTGCTGATCCTTAGCCTGGCCCCCGTCTTCGATTGGACGACCTGCAGGAATCCGTCCACCTCATCCCGGAAAACGAATGAAGCAATGTCCTCGCGGCGCTGCCCGGTtat                          <  2:78260/148‑1 (MQ=255)                                                                                                                      cccaGCTCAGGCTGAT‑GCATTCCAGCGTTACGCTGGTGCTGATCCTTAGCCTGGCCCCCGTCTTCTATTGGACGACCTGCAGGAATCCGTCCACCTCATCCCGGAAAACGAATGAAGCAATGTCCTCGCGGCGCTTACCGGTTATaata                      <  2:378655/146‑3 (MQ=255)                                                                                                                       cAAGCTCAGGCTGAT‑GCATTCCAGCGTTACGCTGGTGCTGATCCTTAGCCTGGCCCCCGTCTTCGATTGGACGACCTGCAGGAATCCGTCCACCTCATCCCGGAAAACGAATGAAGCAATGTCCTCGCGGCGCTGCCCGGTTATAagag                     >  8:185184/1‑149 (MQ=255)                                                                                                                        aaGCTCAGGCTGAT‑GCATTCCAGCGTTACGCTGGTGCTGATCCTTAGCCTGGCCCCCGTCTTCGATTGGACGACCTGCAGGAATCCGTCCACCTCATCCCGGAAAACGAATGAAGCAATGTCCTCGCGGCGCTGCCCGGTTATAAgagc                    >  3:259777/1‑149 (MQ=255)                                                                                                                          gCTCAGGCTGAT‑GCATTCCAGCGTTACGCTGGTGCTGATCCTTAGCCTGGCCCCCGTCTTCGATTGGACGACCTGCAGGAATCCGTCCACCTCATCCCGGAAAACGAATGAAGCAATGTCCTCGCGGCGCTGCCCGGTTATAAgagcga                  <  4:62772/149‑1 (MQ=255)                                                                                                                              aGGCTGAT‑GCATTCCAGCGTTACGCTGGTGCTGATCCTTAGCCTGGCCCCCGTCTTCGATTGGACGACCTGCAGGAATCCGTCCACCTCATCCCGGAAAACGAATGAAGCAATGTCCTCGCGGCGCTGCCCGGTTATAAgagcgagc                <  3:289701/147‑1 (MQ=255)                                                                                                                                     t‑GCATTCCAGCGTTACGCTGGTGCTGATCCTTAGCCTGGCCCCCGTCTTCGATTGGACGACCTGCAGGAATCCGTCCACCTCATCCCGGAAAACGAATGAAGCAATGTCCTCGCGGCGCTGCCCGGTTATAAGAGCGAGCt               <  4:332928/141‑1 (MQ=255)                                                                                                                                     t‑GCATTCCAGCGTTACGCTGGTGCTGATCCTTAGCCTGGCCCCCGTCTTCGATTGGACGACCTGCAGGAATCCGTCCACCTCATCCCGGAAAACGAATGAAGCAATGTCCTCGCGGCGCTGCCCGGTTATAAGAGCGAGCt               >  3:332928/1‑141 (MQ=255)                                                                                                                                     t‑GCATTCCAGCGTTACGCTGGTGCTGATCCTTAGCCTGGCCCCCGTCTTCGATTGGACGACCTGCAGGAATCCGTCCACCTCATCCCGGAAAACGAATGAAGCAATGTCCTCGCGGCGCTGCCCGGTTATAAGAGCGAGCTCCATGgcg       >  2:373255/1‑149 (MQ=255)                                                                                                                                       gCATTCCAGCGTTACGCTGGTGCTGATCCTTAGCCTGGCCCCCGTCTTCGATTGGACGACCTGCAGGAATCCGTCCACCTCATCCCGGAAAACGAATGaa                                                       <  1:374718/100‑1 (MQ=255)                                                                                                                                       gCATTCCAGCGTTACGCTGGTGCTGATCCTTAGCCTGGCCCCCGTCTTCGATTGGACGACCTGCAGGAATCCGTCCACCTCATCCCGGAAAACGAATGaa                                                       >  2:374718/1‑100 (MQ=255)                                                                                                                                       gCATTCCAGCGTTACGCTGGTGCTGATCCTTAGCCTGGCCCCCGTCTTCGATTGGACGACCTGCAGGAATCCGTCCACCTCATCCCGGAAAACGAATGAAGCAATGTCCTCGCGGCGCTGCCCGGTTATAAGAGCGAGCTCCATGgcgc      <  8:357586/149‑1 (MQ=255)                                                                                                                                         aTTCCAGCGTTACGCTGGTGCTGATCCTTAGCCTGGCCCCCGTCTTCGATTGGACGACCTGCAGGAATCCGTCCACCTCATCCCGGAAAACGAATGAAGCAATGTCCTCGCGGCGCTGCCCGGTTATAAGAGCGAGCTCCATGgcgcgc    >  3:230503/1‑149 (MQ=255)                                                                                                                                          ttCCAGCGTTACGCTGGTGCTGATCCTTAGCCTGGCCCCCGTCTTCGATTGGACGCCCTGCAGGAATCCGTCCACCTCATCCCGGAAAACGAATGAAGCAATGTCCTCGCGGCGCTGCCCGGTTATAAGAGCGAGCTCCATGGCGCGCt   >  5:389251/1‑149 (MQ=255)                                                                                                                                          caccAGCGTAACGCTGGTGCTGATCCTTAGCCTGGCCCCCGTCTTCGATTTGACGACCTGCAGGAATCCGTCCACCTCATCCCGGAAAACGAATGAAGCAATGTCCTCGCGGCGCTGCCCGGTTATAAGAGCGAGCTCCATGGCGCGCt   >  6:264280/3‑149 (MQ=255)                                                                                                                                            ccAGCGTTACGCTGGTGCTGATCCTTAGCCTGGCCCCCGTCTTCGATTGGACGACCTGCAGGAATCCGTCCACCTCATCCCGGAAAACGAATGAAGCAATGTCCTCGCGGCGCTGCCCGGTTATAAGAGCGAGCTCCATGGCGCGCtt  <  2:270637/148‑1 (MQ=255)                                                                                                                                              aGCGTTACGCTGGTGCTGATCCTTAGCCTGGCCCCCGTCTTCGATTGGACGACCTGCAGGAATCCGTCCACCTCATCCCGGAAAACGAATGAAGCAATGTCCTCGCGGCGCTGCCCGGTTATAAGAGCGAGCTCCATgg         >  1:50757/1‑139 (MQ=255)                                                                                                                                              aGCGTTACGCTGGTGCTGATCCTTAGCCTGGCCCCCGTCTTCGATTGGACGACCTGCAGGAATCCGTCCACCTCATCCCGGAAAACGAATGAAGCAATGTCCTCGCGGCGCTGCCCGGTTATAAGAGCGAGCTCCATgg         <  2:50757/139‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                             TTTGCTGATTGTGTCGAGCATGATGGGAGATCCCGCCTTCGCGCGGCTGATTGTCCGGTGGTGGTGCACCAGGTGCTTTGACAGCACCCGATCGCGGCAACGCCGT‑ATCACTTCACCCAAGCTCAGGCTGAT‑GCATTCCAGCGTTACGCTGGTGCTGATCCTTAGCCTGGCCCCCGTCTTCGATTGGACGACCTGCAGGAATCCGTCCACCTCATCCCGGAAAACGAATGAAGCAATGTCCTCGCGGCGCTGCCCGGTTATAAGAGCGAGCTCCATGGCGCGCTT  >  NC_002947/1738296‑1738580

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