Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,065,695	G→C	7.3%	R212P (CGC→CCC) ‡	PP_2677 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,065,695	0	G	C	7.3%	218.4 / 5.9	68	R212P (CGC→CCC) ‡	PP_2677	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (31/32);  new base C (2/3);  total (33/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.31e-01

GAGGTGCTGGGCGGGCGCTACCCGCGTGGCGACAGCAGCCGCCTGCGCCTGCAGGTGGCCCACCCCATGGATAACGGCCTGGTCAGCGGCATTCCTGAGTTCTTTCTCAACCATGCCGAGTTACAAGGTGCAGATGGCCAGCGCCTGGCCAGCCTGGAGCTGTACCCGGCGGTCAGCGAAAACC‑CCAACCTGGGCTTTGACATTCAGGGCCCCGGGCCAACACGGCTGCTGTTGCGGGACAACAGCGGCAACCAGTTCGAGGCGGCGCTGCCCTGAGCGCGTTGTCTGATT  >  NC_002947/3065553‑3065843                                                                                                                                               |                                                                                                                                                      gAGGTGCTGGGCGGGCGCTACCCGCGTGGCGACAGCAGCCGCCTGCGCCTGCAGGTGGCCCACCCCATGGATAACGGCCTGGTCAGCGGCATTCCTGAGTTCTTTCTCAACCATGCCGAGTTACAAGGTGCAGATGGCCAGCGCCTGGc                                                                                                                                                 <  6:132052/149‑1 (MQ=255)    gTGCTGGGCGGGCGCTACCCGCGTGGCGACAGCAGCCGCCTGCGCCTGCAGGTGGCCCACCCCATGGATAACGGCCTGGTCAGCGGCATTCCTGAGTTCTTTCTCAACCATGCCGAGTTACAAGGTGCAGATGGCCAGGCGCTGGCCat                                                                                                                                              <  1:293402/149‑2 (MQ=255)     tGCTGGGCGGGCGCTACCCGCGTGGCGACAGCAGCCGCCTGCGCCTGCAGGTGGCCCACCCCATGGATAACGGCCTGGTCAGCGGCATTCCTGAGTTCTTTCTCAACCATGCCGAGTTACAAGGTGCAGATGGCCAGGCGCTGGCCatc                                                                                                                                             <  5:96445/149‑3 (MQ=255)      gCTGGGCGGGCGCTACCCGCGTGGCGACAGCAGCCGCCTGCGCCTGCAGGTGGCCCACCCCATGGATAACGGCCTGGTCAGCGGCATTCCTGAGTTCTTTCTCAACCATGCCGAGTTACAAGGTGCAGATGGCCAGCGCCTGGCCAGcc                                                                                                                                            >  4:123277/1‑149 (MQ=255)           gcgggcgCTACCCGCGTGGCGACAGCAGCCGCCTGCGCCTGCAGGTGGCCCACCCCATGGATAACGGCCTGGTCAGCGGCATTCCTGAGTTCTTTCTCAACCATGCCGAGTTACAAGGTGCAGATGGCCAGCGCCTGGCCAGCCTGGAg                                                                                                                                       <  3:4627/149‑1 (MQ=255)            cgggcgCTACCCGCGTGGCGACAGCAGCCGCCTGCGCCTGCAGGTGGCCCACCCCATGGATAACGGCCTGGTCAGCGGCATTCCTGAGTTCTTTCTCAACCATGCCGAGTTACAAGGTGCAGATGGCCAGCGCCTGGCCAGCCTGGAGc                                                                                                                                      >  7:46824/1‑149 (MQ=255)                 gcTACCCGCGTGGCGACAGCAGCCGCCTGCGCCTGCAGGTGGCCCACCCCATGGATAACGGCCTGGTCAGCGGCATTCCTGAGTTCTTTCTCAACCATGCCGAGTTACTAGGTGCTGATGGCCAGCGCCTGGCCAGCCTGGAGCGGtcg                                                                                                                                 >  7:144499/1‑147 (MQ=255)                 gcTACCCGCGTGGCGACAGCAGCCGCCTGCGCCTGCAGGTGGCCCACCCCATGGATAACGGCCTGGTCAGCGGCATTCCTGAGTTCTTTCTCAACCATGCCGAGTTACAAGGTGCAGATGGCCAGCGCCTGGCCAGCCTGGAGCTGTAc                                                                                                                                 >  3:128419/1‑149 (MQ=255)                 gcTACCCGCGTGGCGACAGCAGCCGCCTGCGCCTGCAGGTGGCCCACCCCATGGATAACGGCCTGGTCAGCGGCATTCCTGAGTTCTTTCTCAACCATGCCGAGTTACAAGGTGCAGATGGCCAGCGCCTGGCCAGCCTGGAGCTGTAc                                                                                                                                 <  3:123277/149‑1 (MQ=255)                    aCCCGCGTGGCGACAGCAGCCGCCTGCGCCTGCAGGTGGCCCACCCCATGGATAACGGCCTGGTCAGCGGCATTCCTGAGTTCTTTCTCAACCATGCCGAGTTACAAGGTGCAGATGGCCAGCGCCTGGCCAGCCTGGAGCTGTACCcg                                                                                                                              >  8:257011/1‑149 (MQ=255)                      ccGCGTGGCGACAGCAGCCGCCTGCGCCTGCAGGTGGCCCACCCCATGGATAACGGCCTGGTCAGCGGCATTCCTGAGTTCTTTCTCAACCATGCCGAGTTACAAGGTGCAGATGGCCAGCGCCTGGCCAGCCTg                                                                                                                                          <  5:236403/135‑1 (MQ=255)                      ccGCGTGGCGACAGCAGCCGCCTGCGCCTGCAGGTGGCCCACCCCATGGATAACGGCCTGGTCAGCGGCATTCCTGAGTTCTTTCTCAACCATGCCGAGTTACAAGGTGCAGATGGCCAGCGCCTGGCCAGCCTg                                                                                                                                          >  6:236403/1‑135 (MQ=255)                       cgcgTGGCGACAGCAGCCGCCTGCGCCTGCAGGTGGCCCACCCCATGGATAACGGCCTGGTCAGCGGCATTCCTGAGTTCTTTCTCAACCATGCCGAGTTACAAGGTGCAGATGGCCAGCGCCTGGCCAGCCTGGAGCTGTACCcggcg                                                                                                                           >  1:62611/1‑149 (MQ=255)                       cgcgTGGCGACAGCAGCCGCCTGCGCCTGCAGGTGGCCCACCCCATGGATAACGGCCTGGTCAGCGGCATTCCTGAGTTCTTTCTCAACCATGCCGAGTTACAAGGTGCAGATGGCCAGCGCCTGGCCAGCCTGGAGCTGTACCcggcg                                                                                                                           <  4:101187/149‑1 (MQ=255)                        gcgTGGCGACAGCAGCCGCCTGCGCCTGCAGGTGGCCCACCCCATGGATAACGGCCTGGTCAGCGGCATTCCTGAGTTCTTTCTCAACCATGCCGAGTTACAAGGTGCAGATGGCCAGCGCCTGGCCAGCCTGGAGCTGTACCcggcgg                                                                                                                          >  4:281887/1‑149 (MQ=255)                         cgTGGCGACAGCAGCCGCCTGCGCCTGCAGGTGGCCCACCCCATGGATAACGGCCTGGTCAGCGGCATTCCTGAGTTCTTTCTCAACCATGCCGAGTTACAAGGTGCAGATGGCCAGCGCCTGGCCAGCCTGGAGCTGTACCCGGCGGt                                                                                                                         <  4:285193/149‑1 (MQ=255)                             ccacctgcaggcgcaggcGGCTGCAGGTGGCCCACCCCATGGATAACGGCCTGGTCAGCGGCATTCCTGAGTTCTTTCTCAACCATGCCGAGTTACAAGGTGCAGATGGCCAGCGCCTGGCCAGCCTGGAGCTGTACCCGGCGGTCAGc                                                                                                                     <  2:249100/143‑1 (MQ=255)                                 cagcagCCGCCTGCGCCTGCAGGTGGCCCACCCCATGGATAACGGCCTGGTCAGCGGCATTCCTGAGTTCTTTCTCAACCATGCCGAGTTACAAGGTGCAGATGGCCAGCGCCTGGCCAGCCTGGAGCTGTACCCGGCGGTCAGCGaaa                                                                                                                 <  4:128419/149‑1 (MQ=255)                                   gcagcCGCCTGCGCCTGCAGGTGGCCCACCCCATGGATAACGGCCTGGTCAGCGGCATTCCTGAGTTCTTTCTCAACCATGCCGAGTTACAAGGTGCAGATGGCCAGCGCCTGGCCAGCCTGGAGCTGTACCCGGCGGTCAGCGAAAAc                                                                                                               >  1:352518/1‑149 (MQ=255)                                      gccgccTGCGCCTGCAGGTGGCCCACCCCATGGATAACGGCCTGGTCAGCGGCATTCCTGAGTTCTTTCTCAACCATGCCGAGTTACAAGGTGCAGATGGCCAGCGCCTGGCCAGCCTGGAGCTGTACCCGGCGGTCAGCGAAAAcc‑cc                                                                                                           <  3:370643/149‑1 (MQ=255)                                      gccgccTGCGCCTGCAGGTGGCCCACCCCATGGATAACGGCCTGGTCAGCGGCATTCCTGAGTTCTTTCTCAACCATGCCGAGTTACAAGGTGCAGATGGCCAGCGCCTGGCCAGCCTGGAGCTGTACCCGGCGGTCAGCGAAAAcc‑c                                                                                                            >  5:81545/1‑148 (MQ=255)                                              cgcctgcAGGTGGCCCACCCCATGGATAACGGCCTGGTCAGCGGCATTCCTGAGTTCTTTCTCAACCATGCCGAGTTACAAGGTGCAGATGGCCAGCGCCTGGCCAGCCTGGAGCTGTACCCGGCGGTCAGCGAAAACC‑CCAACCTggg                                                                                                   <  1:325268/149‑1 (MQ=255)                                              cgcctgcAGGTGGCCCACCCCATGGATAACGGCCTGGTCAGCGGCATTCCTGAGTTCTTTCTCAACCATGCCGAGTTACAAGGTGCAGATGGCCAGCGCCTGGCCAGCCTGGAGCTGTACCCGGCGGTCAGCGAAAACC‑CCAACCTgg                                                                                                    <  8:46824/148‑1 (MQ=255)                                               gcctgcAGGTGGCCCACCCCATGGATAACGGCCTGGTCAGCGGGATTCCTGAGTTCTTTCTCAACCATGCCGAGTTACAAGGTGCAGATGGCCAGCGCCTGGCCAGCCTGGAGCTGTACCCGGCGGTCAGCGAAAACC‑CCAACCTGGGc                                                                                                  >  8:106150/1‑149 (MQ=255)                                               gcctgcAGGTGGCCCACCCCATGGATAACGGCCTGGTCAGCGGCATTCCTGAGTTCTTTCTCAACCATGCCGAGTTACAAGGTGCAGATGGCCAGCGCCTGGCCAGCCTGGAGCTGTACCCGGCGGTCAGCGAAAACC‑CCAACCTGGGc                                                                                                  <  7:257011/149‑1 (MQ=255)                                                 ctgcAGGTGGCCCACCCCATGGATAACGGCCTGGTCAGCGGCATTCCTGAGTTCTTTCTCAACCATGCCGAGTTACAAGGTGCAGATGGCCAGCGCCTGGCCAGCCTGGAGCTGTACCCGGCGGTCAGCGAAAACC‑CCAACCTGGGc                                                                                                  <  5:66727/147‑1 (MQ=255)                                                 ctgcAGGTGGCCCACCCCATGGATAACGGCCTGGTCAGCGGCATTCCTGAGTTCTTTCTCAACCATGCCGAGTTACAAGGTGCAGATGGCCAGCGCCTGGCCAGCCTGGAGCTGTACCCGGCGGTCAGCGAAAACC‑CCAACCTGGGc                                                                                                  >  6:66727/1‑147 (MQ=255)                                                 ctgcAGGTGGCCCACCCCATGGATAACGGCCTGGTCAGCGGCATTCCTGAGTTCTTTCTCAACCATGCCGAGTTACAAGGTGCAGATGGCCAGCGCCTGGCCAGCCTGGAGCTGTACCCGGCGGTCAGCGAAAACC‑CCAACCTGGGCtt                                                                                                >  7:61119/1‑149 (MQ=255)                                                    cAGGTGGCCCACCCCATGGATAACGGCCTGGTCAGCGGCATTCCTGAGTTCTTTCTCAACCATGCCGAGTTACAAGGTGCAGATGGCCAGCGCCTGGCCAGCCTGGAGCTGTACCCGGCGGTCAGCGAAAACC‑CCAACCTGGGCTTTGa                                                                                             <  7:197268/149‑1 (MQ=255)                                                       gTGGCCCACCCCATGGATAACGGCCTGGTCAGCGGCATTCCTGAGTTCTTTCTCAACCATGCCGAGTTACAAGGTGCAGATGGCCAGCGCCTGGCCAGCCTGGAg                                                                                                                                       >  6:341261/1‑105 (MQ=255)                                                       gTGGCCCACCCCATGGATAACGGCCTGGTCAGCGGCATTCCTGAGTTCTTTCTCAACCATGCCGAGTTACAAGGTGCAGATGGCCAGCGCCTGGCCAGCCTGGAg                                                                                                                                       <  5:341261/105‑1 (MQ=255)                                                       gTGGCCCACCCCATGGATAACGGCCTGGTCAGCGGCATTCCTGAGTTCTTTCTCAACCATGCCGAGTTACAAGGTGCAGATGGCCAGCGCCTGGCCAGCCTGGAGCTGTACCCGGCGGTCAGCGAAAACC‑CCAACCTGGGCTTTGACAt                                                                                          >  4:29235/1‑149 (MQ=255)                                                          gCCCACCCCATGGATAACGGCCTGGTCAGCGGCATTCCTGAGTTCTTTCTCAACCATGCCGAGTTACAAGGTGCAGATGGCCAGCGCCTGGCCAGCCTGGAGCTGTACCCGGCGGTCAGCGAAAACC‑CCAACCTGGGCTTTGACAtt                                                                                         >  6:194418/1‑147 (MQ=255)                                                            ccACCCCATGGATAACGGCCTGGTCAGCGGCATTCCTGAGTTCTTTCTCAACCATGCCGAGTTACAAGGTGCAGATGGCCAGCGCCTGGCCAGCCTGGAGCTGTACCCGGCGGTCAGCGAAAACC‑CCAACCTGGGCTTTGACATTCAgg                                                                                     <  3:281887/149‑1 (MQ=255)                                                                     ggATAACGGCCTGGTCAGCGGCATTCCTGAGTTCTTTCTCAACCATGCCGAGTTACAAGGTGCAGATGGCCAGCGCCTGGCCAGCCTGGAGCTGTACCCGGCGGTCAGCGAAAACC‑CCAACCTGGGCTTTGACATTCAgg                                                                                     >  2:8926/1‑140 (MQ=255)                                                                     ggATAACGGCCTGGTCAGCGGCATTCCTGAGTTCTTTCTCAACCATGCCGAGTTACAAGGTGCAGATGGCCAGCGCCTGGCCAGCCTGGAGCTGTACCCGGCGGTCAGCGAAAACC‑CCAACCTGGGCTTTGACATTCAgg                                                                                     <  1:8926/140‑1 (MQ=255)                                                                       aTAACGGCCTCGTCAGCGGCAATACTGAGATCATTCTCAACAATGCCGAGTAACAAGGTGCAGATTGCCAGCGCCTGGCCAGCCTGGAGCAGTACCCGGCGGTCAGCGAAAACC‑CCAACCTGGGCTTTGACATTCAGCGCCCCCGGca                                                                           >  2:286945/1‑147 (MQ=255)                                                                        tAACGGCCTGGTCAGTGGCATTCCTGAGTTCTTTCTCAACCATGCCGAGTTACAAGGTGCAGATGGCCAGCGCCTGGCCAGCCTGGAGCTGTACCCGGCGGTCAGCGAAAACC‑CCAACCTg                                                                                                     <  7:317172/121‑1 (MQ=255)                                                                        tAACGGCCTGGTCAGTGGCATTCCTGAGTTCTTTCTCAACCATGCCGAGTTACAAGGTGCAGATGGCCAGCGCCTGGCCAGCCTGGAGCTGTACCCGGCGGTCAGCGAAAACC‑CCAACCTg                                                                                                     >  8:317172/1‑121 (MQ=255)                                                                          aCGGCCTGGTCAGCGGCATTCCTGAGTTCTTTCTCAACCATGCCGAGTTACAAGGTGCAGATGGCCAGCGCCTGGCCAGCCTGGAGCTGTACCCGGCGGTCAGCGAAAACC‑CCAACCTGGGCTTTGACATTCAGGGCCCCGGGCCAaca                                                                       <  6:81545/149‑1 (MQ=255)                                                                                      gCGGCATTCCTGAGTTCTTTCTCAACCATGCCGAGTTACAAGGTGCAGATGGCCAGCGCCTGGCCAGCCTGGAGCTGTACCCGGCGGTCAGCGAAAACC‑CCAACCTGGGCTTTGACATTCAGGGCCCCGGGCCAACACgg                                                                    >  6:204698/1‑140 (MQ=255)                                                                                      gCGGCATTCCTGAGTTCTTTCTCAACCATGCCGAGTTACAAGGTGCAGATGGCCAGCGCCTGGCCAGCCTGGAGCTGTACCCGGCGGTCAGCGAAAACC‑CCAACCTGGGCTTTGACATTCAGGGCCCCGGGCCAACACgg                                                                    <  5:204698/140‑1 (MQ=255)                                                                                       cGGCATTCCTGAGTTCTTTCTCAACCATGCCGAGTTACAAGGTGCTGATGGCCAGCGCCTGGCCAGCCTGGAGCTGTACCCGGCGGTCAGCGAAAACC‑CCAACCTGGGCTTTGACATTCAGGGCCCCGGGCCAACACGGCTGCTGTTGc                                                          >  1:270254/1‑149 (MQ=255)                                                                                            ttCCTGCGTTCTTTCTCCACCATGCCGAGGTACACGGTGCAGAGGGCCAGCGCCTGGCCAGCCTGGAGCTGTACCCGGCGGTCAGCGAAAACC‑CCAACCTGGGCGTTGACATTCAGGGCCCCGGGCCAACACGGCTGCTGTTGCGGGac                                                     <  2:62611/149‑1 (MQ=255)                                                                                               cTGAGTTCTTTCTCAACCATGCCGAGTTACAAGGTGCAGATGGCCAGCGCCTGGCCAGCCTGGAGCTGTACCCGGCGGTCAGCGAAAACC‑CCAACCTGGGCTTTGACATTCAGGGCCCCGGGCCAACACGGCTGCTGTTGCTGGacaa                                                   >  8:342453/1‑148 (MQ=255)                                                                                                tGAGTTCTTTCTCAACCATGCCGAGTTACAAGGTGCAGATGGCCAGCGCCTGGCCAGCCTGGAGCTGTACCCGGCGGTCAGCGAAAACC‑CCAACCTGGGCTTTGACATTCAGGGCCCCGGGCCAACACGGCTGCTGTTGCGGGacaa                                                   >  1:358365/1‑147 (MQ=255)                                                                                                  aGTTCTTTCTCAACCATGCCGAGTTACAAGGTGCAGATGGCCAGCGCCTGGCCAGCCTGGAGCTGTACCCGGCGGTCAGCGAAAACC‑CCAACCTGGGCTTTGACATTCAGGGCCCCGGGCCAACACGGCTGCTGTTGCGGGacaaca                                                 <  4:317600/147‑1 (MQ=255)                                                                                                   gTTCTTTCTCAACCATGCCGAGTTACAAGGTGCAGATGGCCAGCGCCTGGCCAGCCTGGAGCTGTACCCGGCGGTCAGCGAAAACC‑CCAACCTGGGCTTTGACATTCAGGGCCCCGGGCCAACACGGCTGCTGTTGCGGGACAACAGCg                                              <  3:29235/149‑1 (MQ=255)                                                                                                     tctttctCAACCATGCCGAGTTACAAGGTGCAGATGGCCAGCGCCTGGCCAGCCTGGAGCTGTACCCGGCGGTCAGCGAAAACC‑CCAACCTGGGCTTTGACATTCAGGGCCCCGGGCCAACACGGCTGCTGTTGCGGGACAACAGCgg                                             <  2:358365/148‑1 (MQ=255)                                                                                                         tctcAACCATGCCGAGTTACAAGGTGCAGATGGCCAGCGCCTGGCCAGCCTGGAGCTGTACCCGGCGGTCAGCGAAAACC‑CCAACCTGGGCTTTGACATTCAGGGCCCCGGGCCAACACGGCTGCTGTTGCGGGACAACAGCGGCa                                           <  1:286945/146‑1 (MQ=255)                                                                                                         tctcAACCATGCCGAGTTACAAGGTGCAGATGGCCAGCGCCTGGCCAGCCTGGAGCTGTACCCGGCGGTCAGCGAAAACC‑CCAACCTGGGCTTTGACATTCAGGGCCCCGGGCCAACACGGCTGCTGTTGCGGGACAACAGCGGCa                                           <  7:106150/146‑1 (MQ=255)                                                                                                             aaCCATGCCGAGTTACAAGGTGCAGATGGCCAGCGCCTGGCCAGCCTGGAGCTGTACCCGGCGGTCAGCGAAAACCACCAACCTGGGCTTTGACATTCAGGGCCCCGGGc                                                                            <  1:45296/110‑1 (MQ=255)                                                                                                             aaCCATGCCGAGTTACAAGGTGCAGATGGCCAGCGCCTGGCCAGCCTGGAGCTGTACCCGGCGGTCAGCGAAAACC‑CCAACCTGGGCTTTGACATTCAGGGCCCCGGGc                                                                            >  2:45296/1‑109 (MQ=255)                                                                                                              aCCATGCCGAGTTACAAGGTGCAGATGGCCAGCGCCTGGCCAGCCTGGAGCTGTACCCGGCGGTCAGCGAAAACC‑CCAACCTGGGCTTTGACa                                                                                           <  1:184623/93‑1 (MQ=255)                                                                                                              aCCATGCCGAGTTACAAGGTGCAGATGGCCAGCGCCTGGCCAGCCTGGAGCTGTACCCGGCGGTCAGCGAAAACC‑CCAACCTGGGCTTTGACa                                                                                           >  2:184623/1‑93 (MQ=255)                                                                                                                  tGCCGAGTTACAAGGTGCAGATGGCCAGCGCCTGGCCAGCCTGGAGCTGTACCCGGCGGTCAGCGAAAACC‑CCAACCTGGGCTTTGACATTCAGGGCCCCGGGCCAACACGGCTGCTGTTGCGGGACAACAGCGGCAACCAGTTCGAgg                               >  2:366328/1‑149 (MQ=255)                                                                                                                           aCAAGGTGCAGATGGCCAGCGCCTGGCCAGCCTGGAGCTGTACCCGGCGGTCAGCGAAAACC‑CCAACCTGGGCttt                                                                                               >  1:258558/1‑76 (MQ=255)                                                                                                                           aCAAGGTGCAGATGGCCAGCGCCTGGCCAGCCTGGAGCTGTACCCGGCGGTCAGCGAAAACC‑CCAACCTGGGCttt                                                                                               <  2:258558/76‑1 (MQ=255)                                                                                                                            cAAGGTGCAGATGGCCAGCGCCTGGCCAGCCTGGAGCTGTACCCGGCGGTCAGCGAAAACC‑CCAACCTGGGCTTTGACATTCAGGGCCCCGGGCCAACACGGCTGCTGTTGCGGGACAACAGCGGCAACCAGTTCGAGGCGGCGCTGcc                     <  2:352518/149‑1 (MQ=255)                                                                                                                                gTGCAGATGGCCAGCGCCTGGCCAGCCTGGAGCTGTACCCGGCGGTCAGCGAAAACC‑CCAACCTGGGCTTTGACATTCAGGGCCCCGGGCCAACACGGCTGCTGTTGCGTGACAACAGCGGCAACCAGTTCGAGGCGGCGCTGCCCTGa                 <  8:56425/149‑1 (MQ=255)                                                                                                                                gTGCAGATGGCCAGCGCCTGGCCAGCCTGGAGCTGTACCCGGCGGTCAGCGAAAACC‑CCAACCTGGGCTTTGACATTCAGGGCCCCGGGCCAACACGGCTGCTGTTGCGGGACAACAGCGGCAACCAGTTCGAGGCGGCGCTGCCCt                   >  6:199237/1‑147 (MQ=255)                                                                                                                                 ccaggctGGCCAGGCGCTGGCCATCCTGGAGCTGTACCCGGCGGTCAGCGAAAACC‑CCAACCTGGGCTTTGACATTCAGGGCCCCGGGCCAACCCGGCTGCTGTTGCGGGACAACAGCGGCAACCAGTTCGAGGCGGCGCTGCCCTGAg                >  6:274758/7‑149 (MQ=255)                                                                                                                                   cAGATGGCCAGCGCCTGGCCAGCCTGGAGCTGTACCCGGCGGTCAGCGAAAACC‑CCAACCTGGGCTTTGACATTCAGGGCCCCGGGCCAACACGGCTGCTGTTGCGGGACAACAGCGGCAACCAGTTCGAGGCGGCGCTGCCCTGAgc               <  5:194418/148‑1 (MQ=255)                                                                                                                                   aggctGGCCAGGCGCTGGCCAGCCTGGAGCTGTACCCGGCGGTCAGCGAAAACC‑CCAACCTGGGCTTTGACATTCAGGGCCCCGGGCCAACACGGCTGCTGTTGCGGGACAACAGCGGCAACCAGTTCGAGCCGGCGCTGCCCTGa                 <  7:324535/142‑1 (MQ=255)                                                                                                                                   aggctGGCCAGGCGCTGGCCAGCCTGGAGCTGTACCCGGCGGTCAGCGAAAACC‑CCAACCTGGGCTTTGACATTCAGGGCCCCGGGCCAACACGGCTGCTGTTGCGGGACAACAGCGGCAACCAGTTCGAGCCGGCGCTGCCCTGa                 >  8:324535/5‑146 (MQ=255)                                                                                                                                     gATGGCCAGCGCCTGGCCAGCCTGGAGCTGTACCCGGCGGTCAGCGAAAACC‑CCAACCTGGGCTTTGACATTCAGGGCCCCGGGCCAACACGGCTGCTGTTGCGGGACAACAGCGGCAACCAGTTCGAGGCGGCGCTGCCCTGAGcgcg            >  2:313982/1‑149 (MQ=255)                                                                                                                                           cAGCGCCTGGCCAGCCTGGAGCTGTACCCGGCGGTCAGCGAAAACC‑CCAACCTGGGCTTTGACATTCAGGGCCCCGGGCCAACACGGCTGCTGTTGCGGGACAACAGCGGCAACCAGTTCGAGGCGGCGCTGCCCTGAGCGCGTTGTCt      >  7:153040/1‑149 (MQ=255)                                                                                                                                           cAGCGCCTGGCCAGCCTGGAGCTGTACCCGGCGGTCAGCGAAAACC‑CCAACCTGGGCTTTGACATTCAGGGCCCCGGGCCAACACGGCTGCTGTTGCGGGACAACAGCGGCAACCAGTTCGAGGCGGCGCTGCCCTGAGCGCGTTGTCt      >  7:126316/1‑149 (MQ=255)                                                                                                                                              cgcCTGGCCAGCCTGGAGCTGTACCCGGCGGTCAGCGAAAACC‑CCAACCTGGGCTTTGACATTCAGGGCCCCGGGCCAACACGGCTGCTGTTGCGGGACAACAGCGGCAACCAGTTCGAGGCGGCGCTGCCCTGAGCGCGTTGTCTGAt   <  1:313982/149‑1 (MQ=255)                                                                                                                                               gcCTGGCCAGCCTGGAGCTGTACCCGGCGGTCAGCGAAAACC‑CCAACCTGGGCTTTGACATTCAGGGCCCCGGGCCAACACGGCTGCTGTTGCGGGACAACAGCGGCAACCAGTTCGAGGCGGCGCTGCCCTGAGCGCGTTGTCTGAtt  >  3:107062/1‑149 (MQ=255)                                                                                                                                               |                                                                                                                                                      GAGGTGCTGGGCGGGCGCTACCCGCGTGGCGACAGCAGCCGCCTGCGCCTGCAGGTGGCCCACCCCATGGATAACGGCCTGGTCAGCGGCATTCCTGAGTTCTTTCTCAACCATGCCGAGTTACAAGGTGCAGATGGCCAGCGCCTGGCCAGCCTGGAGCTGTACCCGGCGGTCAGCGAAAACC‑CCAACCTGGGCTTTGACATTCAGGGCCCCGGGCCAACACGGCTGCTGTTGCGGGACAACAGCGGCAACCAGTTCGAGGCGGCGCTGCCCTGAGCGCGTTGTCTGATT  >  NC_002947/3065553‑3065843

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