Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,655,522	C→G	7.8%	G139G (GGC→GGG) ‡	PP_4118 →	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,655,522	0	C	G	7.8%	197.0 / 6.3	64	G139G (GGC→GGG) ‡	PP_4118	membrane protein
Reads supporting (aligned to +/- strand):  ref base C (31/28);  new base G (2/3);  total (33/31)
Fisher's exact test for biased strand distribution p-value = 6.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.67e-01

CCTGGCCGCCTTCCTTTGCGCCCACCTTGCTTACCTGCGCGGCTATTGCGGCATTACCCTGCGCCCCGCCCTGCCC‑GCACTGATCTTCAGCGCAATCACTGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGCCCTGGCCATCAGCGCCATGCTCTGGCGCGCACTGGCCTGCGGGGGAGTGGCAGCACTGGGTGCCGGCCTGTTCGTGTTCTCCGACAGCCTGATCGGCATCGACCGCTTCATCAG  >  NC_002947/4655381‑4655665                                                                                                                                               |                                                                                                                                                ccTGGCCGCCTTCCTTTGCGCCCACCTTGCTTACCTGCGCGGCTATTGCGGCATTACCCTGCGCCCCGCCCTGCCC‑GCACTGATCTTCAGCGCAATCACTGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCTCGGCCCGTTGc                                                                                                                                          <  2:345226/149‑1 (MQ=255)  cTGGCCGCCTTCCTTTGCGCCCACCTTGCTTACCTGCGCGGCTATTGCGGCATTACCCTGCGCCCCGCCCTGCCC‑GCACTGATCTTCAGCGCAATCACTGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCTCGGCCCGTTGCt                                                                                                                                         >  1:167486/1‑149 (MQ=255)     gccgccTTCCTTTGCGCCCACCTTGCTTACCTGCGCGGCTATTGCGGCATTACCCTGCGCCCCGCCCTGCCC‑GCACTGATCTTCAGCGCAATCACTGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGGCCGAGGCCGTTGCTg                                                                                                                                        <  6:131826/147‑1 (MQ=255)     gccgccTTCCTTTGCGCCCACCTTGCTTACCTGCGCGGCTATTGCGGCATTACCCTGCGCCCCGCCCTGCCC‑GCACTGATCTTCAGCGCAATCACTGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCTCGGCCCGTTGCTGa                                                                                                                                       <  4:162913/148‑1 (MQ=255)     gccgccTTCCTTTGCGCCCACCTTGCTTACCTGCGCGGCTATTGCGGCATTACCCTGCGCCCCGCCCTGCCC‑GCACTGATCTTCAGCGCAATCACTGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCTCGGCCCGTTGCTGa                                                                                                                                       <  5:134962/148‑1 (MQ=255)      ccgccTTCCTTTGCGCCCACCTTGCTTACCTGCGCGGCTATTGCGGCATTACCCTGCGCCCCGCCCTGCCC‑GCACTGATCTTCAGCGCAATCACTGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCTCGGCCCGTTGCTGATc                                                                                                                                     >  3:20977/1‑149 (MQ=255)      ccgccTTCCTTTGCGCCCACCTTGCTTACCTGCGCGGCTATTGCGGCATTACCCTGCGCCCCGCCCTGCCC‑GCACTGATCTTCAGCGCAATCACTGGCATCGCGCTGCTCGGCGGACTGGCCAGCAACGGCCTCGGCCCGTTGCTGATc                                                                                                                                     <  3:125226/149‑1 (MQ=255)           ttccttTGCGCCCACCTTGCTTACCTGCGCGGCTATTGCGGCATTACCCTGCGCCCCGCCCTGCCC‑GCACTGATCTTCAGCGCAATCACTGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCTCGGCCCGTTGCTGATCCCgg                                                                                                                                 >  2:226865/1‑148 (MQ=255)              cttTGCGCCCACCTTGCTTACCTGCGCGGCTATTGCGGCATTACCCTGCGCCCCGCCCTGCCC‑GCACTGATCTTCAGCGCAATCACTGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCTCGGCCCCTTGCTGATCCCGGTAg                                                                                                                              >  7:143031/1‑148 (MQ=255)                    gcCCACCTTGCTTACCTGCGCGGCTATTGCGGCATTACCCTGCGCCCCGCCCTGCCC‑GCACTGATCTTCAGCGCAATCACTGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTa                                                                                                                       >  8:160268/1‑149 (MQ=255)                      ccACCTTGCTTACCTGCGCGGCTATTGCGGCATTACCCTGCGCCCCGCCCTGCCC‑GCACTGATCTTCAGCGCAATCACTGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACg                                                                                                                     >  3:269256/1‑149 (MQ=255)                        aCCTTGCTTACCTGCGCGGCTATTGCGGCATTACCCTGCGCCCCGCCCTGCCC‑GCACTGATCTTCAGCGCAATCACTGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGcc                                                                                                                   >  5:142634/1‑149 (MQ=255)                             gCTTACCTGCGCGGCTATTGCGGCATTACCCTGCGCCCCGCCCTGCCC‑GCACTGATCTTCAGCGCAATCACTGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGCCCTg                                                                                                                <  4:12490/147‑1 (MQ=255)                               ttACCTGCGCGGCTATTGCGGCATTACCCTGCGCCCCGCCCTGCCCGGCACTGATCATCAGCGCAATCACTGGCATCGCGCTGCTCGGCTTACTGGCCAGCAACGGCCTCGGCCCGTTGGTGATCCCGGTAGCACTCAACGCCCAGGc                                                                                                              >  1:249324/1‑148 (MQ=255)                                  ccTGCGCGGCTATTGCGGCATTACCCTGCGCCCCGCCCTGCCC‑GCACTGATCTTCAGCGCAATCACTGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGCCCTGGCCATc                                                                                                          <  3:58042/148‑1 (MQ=255)                                           cTATTGCGGCATTACCCTGCGCCCCGCCCTGCCC‑GCACTGATCTTCAGCGCAATCACTGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGCCCTGGCCATCAGCGCCATGc                                                                                                >  8:167437/1‑149 (MQ=255)                                                gCGGCATTACCCTGCGCCCCGCCCTGCCC‑GCACTGATCTTCAGCGCAATCACTGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGCCCTGGCCAt                                                                                                           >  1:66797/1‑133 (MQ=255)                                                gCGGCATTACCCTGCGCCCCGCCCTGCCC‑GCACTGATCTTCAGCGCAATCACTGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGCCCTGGCCAt                                                                                                           <  2:66797/133‑1 (MQ=255)                                                    cATTACCCTGCGCCCCGCCCTGCCC‑GCACTGATCTTCAGCGCAATCACTGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGCCCTGGCCATCAGCGCCATGCTCTGgcgcg                                                                                       >  5:244671/1‑149 (MQ=255)                                                     aTTACCCTGCGCCCCGCCCTGCCC‑GCACTGATCTTCAGCGCAATCACTGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGGCCGAGGCCGTTGCTg                                                                                                                                        <  1:59605/99‑1 (MQ=38)                                                     aTTACCCTGCGCCCCGCCCTGCCC‑GCACTGATCTTCAGCGCAATCACTGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGGCCGAGGCCGTTGCTg                                                                                                                                        >  2:59605/1‑99 (MQ=38)                                                          ccTGCGCCCCGCCCTGCCC‑GCACTGATCTTCAGCGCAATCACTGGCATCGCGCTGCTCGGCGTACTGGCCAGCAATGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGCCCTGGCCATCAGCGCCATGCTCTGGCGCGCACTgg                                                                                 >  6:65090/1‑149 (MQ=255)                                                          ccTGCGCCCCGCCCTGCCC‑GCACTGATCTTCAGCGCAATCACTGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTAc                                                                                                                      <  2:116016/112‑1 (MQ=255)                                                          ccTGCGCCCCGCCCTGCCC‑GCACTGATCTTCAGCGCAATCACTGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGCCCTGGCCATCAGCGCCATGCTCTGGCGCGCACTgg                                                                                 >  2:269441/1‑149 (MQ=255)                                                           cTGCGCCCCGCCCTGCCC‑GCACTGATCTTCAGCGCAATCACTTGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGCCCTGGCCATCAGCGCCATGCTCTGGCGCGCACTGGc                                                                                <  4:20977/149‑1 (MQ=255)                                                           cTGCGCCCCGCCCTGCCC‑GCACTGATCTTCAGCGCAATCACTGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGCCCTGGCCATCAGCGCCATGCTCTGGCGCGCACTGGc                                                                                <  4:269256/149‑1 (MQ=255)                                                            tGCGCCCCGCCCTGCCC‑GCACTGATCTTCAGCGCAATCACTGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGccc                                                                                                                  <  1:141965/114‑1 (MQ=255)                                                            tGCGCCCCGCCCTGCCC‑GCACTGATCTTCAGCGCAATCACTGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGccc                                                                                                                  >  2:141965/1‑114 (MQ=255)                                                                 cccgcccTGCCC‑GCACTGATCTTCAGCGCAATCACTGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGCCCTGGCCATCAGCGCCATGCTCTGGCGCGCACTGGCCTGCgg                                                                          >  2:221056/1‑149 (MQ=255)                                                                      ccTGCCC‑GCACTGATCTTCAGCGCAATCACTGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGCCCTGGCCATCAGCGCCATGCTCTGGCGCGCACTGGCCTGCgggg                                                                        >  1:216855/1‑146 (MQ=255)                                                                      ccTGCCC‑GCACTGATCTTCAGCGCAATCACTGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGCCCTGGCCATCAGCGCCATGCTCTGGCGCGCACTGGCCTGCGGGGGa                                                                      <  2:216855/148‑1 (MQ=255)                                                                           cc‑GCACTGATCTTCAGCGCAATCACTGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCTCGGCCCGTTGCTGATGCCGGTAGCACTCTACGCCCTGGCCATCAGCGCCATGCTCTGGCGAGCACTGGCCTGCGGGGGAGTGgca                                                                <  8:143031/149‑1 (MQ=255)                                                                                  tGATCTTCAGCGCAATCACTGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGCCCTGGCCATCAGCGCCATGCTCTGGCGCGCACTGGCCTGCGGGGGAGTGGCAGCACTg                                                          >  3:273901/1‑149 (MQ=255)                                                                                    aTCTTCAGCGCAATCACTGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGCCCTGGCCATCAGCGCCATGCTCTGGCGCGCACTGGCCTGCGGGGGAGTGGCAGCACTggg                                                        >  3:38549/1‑149 (MQ=255)                                                                                       ttCAGCGCAATCACTGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCACGGCCCGTTGAGGATCCCGGTAACACACAACGCCCAGGCCAACAGCGCAATGAAAGGGCGCGCACGGGCAAGCGGGGGAGTGGCCGCAAAGGGTGc                                                     >  2:309104/1‑149 (MQ=255)                                                                                        tCAGCGCAATCACTGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGCCCTGGCCATCAGCGCCATGCTCTGGCGCGCACTGGCCTGCGGGGGAGTGGCAGCACTGGGTGcc                                                    >  3:65921/1‑149 (MQ=255)                                                                                             gcAATCACTGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGCCCTGGCCATCAGCGCCATGCTCTGGCGCGCACTGGCCTGCGGGGGAGTGGCAGCACTGGGTGCCGGCCt                                               <  7:167437/149‑1 (MQ=255)                                                                                              cAATCACTGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGCCCTGGCCATCAGCGCCATGCTCTGGCGCGCACTGGCCTGCGGGGGAGTGGCAGCACTGGGTGCCg                                                   >  8:304006/1‑144 (MQ=255)                                                                                              cAATCACTGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGCCCTGGCCATCAGCGCCATGCTCTGGCGCGCACTGGCCTGCGGGGGAGTGGCAGCACTGGGTGCCg                                                   <  7:304006/144‑1 (MQ=255)                                                                                                 tCACTGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGCCCTGGCCATCAGCGCCATGCTCTGGCGCGCACTGGCCTGCGGGGGAGTGGCAGCACTGGGTGCCGGCCTGTTc                                           <  6:142634/149‑1 (MQ=255)                                                                                                    cTGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGCCCTGGCCATCAGCGCCATGCTCTGGCGCGCACTGGCCTGCGGGGGAGTGGCAGCACTGGGTGCCGGCCTGTTCgtg                                        <  4:65921/149‑1 (MQ=255)                                                                                                     tGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGCCCTGGCCATCAGCGCCATGCTCTGGCGCGCACTGGCCTGCGGGGGAGTGgcag                                                               <  1:339310/125‑1 (MQ=255)                                                                                                     tGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGCCCTGGCCATCAGCGCCATGCTCTGGCGCGCACTGGCCTGCGGGGGAGTGgcag                                                               >  2:339310/1‑125 (MQ=255)                                                                                                      ggCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGCCCTGGCCATCAGCGCCATGCTCTGGCGCGCACTGGCCTGCGGGGGAGTGGCAGCACTGGGTGCCGGCCTGTTCGTGtt                                      <  4:273901/149‑1 (MQ=255)                                                                                                         aTCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGCCCTGGCCATCAGCGCCATGCTCTGGCGCGCACTGGCCTGCGGGGGAGTGGCAGCACTGGGTGCCGGCCTGTTCGTGTtctc                                   >  7:72269/1‑149 (MQ=255)                                                                                                             cgcTGCTCGGCGTACTGGCCAGCAATGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGCCCTGGCCATCAGCGCCATGCTCTGGCGCGCACTGGCCTGCGGGGGAGTGGCAGCACTGGGTGCCGGCCTGTTCGTGTTCTCCGAc                               <  5:65090/149‑1 (MQ=255)                                                                                                                  ctCGGCGTACTGGCCAGCAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGCCCTTGCCATCAGCGCCATGCTCTGGCGCGCACTTGCCTGCGGGTGAGTGGCAGCACTGGGTGCCGGCCTGTTCGTGTTCTCCGACAGc                            <  2:167486/147‑1 (MQ=255)                                                                                                                     ggCGTACTGGCCAGCAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGCCCTGGCCATCAGCGCCATGCTCTGGCGCGCACTGGCCTGCGGGGGAGTGGCAGCACTGGGTGCCGGCCTGTTCGTGTTCTCCGACAGCCTGAt                       <  4:38549/149‑1 (MQ=255)                                                                                                                     ggCGTACTGGCCAGCAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGCCCTGGCCATCAGCGCCATGCTCTGGCGCGCACTGGCCTGCGGGGGAGTGGCAGCACTGGGTGCCGGCCTGTTCGTGTTCTCCGACAGCCTGAt                       >  6:150424/1‑149 (MQ=255)                                                                                                                      gCGTACTGGCCAGCAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGCCCTGGCCATCAGCGCCATGCTCTGGCGCGCACTGGCCTGCGGGGGAGTGGCAGCACTGGGTGCCGGCCTGTTCGTGTTCTCCGACAGCCTGATc                      <  8:72269/149‑1 (MQ=255)                                                                                                                       cGTACTGGCCAGCAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGCCCTGGCCATCAGCGCCATGCTCTGGCGCGCACTTGCCTGCGGGGGAGTGGCAGCACTGGGTGCCGGCCTGTTCGTGTTCTCCGACAGCCTGATCg                     <  1:269441/149‑1 (MQ=255)                                                                                                                             gatcAGCAACGGGCCGAGGCCGTTGCTGGTCCCGGTAGCACTCTACGCCCTGGCCATCAGCGCCATGCTCTGGCGCGCACTGGCCTGCGGGGGAGTGGCAGCACTGGGTGCCGGCCTGTTCGTGTTCTcc                                  >  7:39433/4‑130 (MQ=255)                                                                                                                             gatcAGCAACGGGCCGAGGCCGTTGCTGGTCCCGGTAGCACTCTACGCCCTGGCCATCAGCGCCATGCTCTGGCGCGCACTGGCCTGCGGGGGAGTGGCAGCACTGGGTGCCGGCCTGTTCGTGTTCTcc                                  <  8:39433/127‑1 (MQ=255)                                                                                                                               ccAGCAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGCCCTGGCCATCAGCGCCATGCTCTGGCGCGCACTGGCCTGCGGGGGAGTGGCAGCACTGGGTGCCGGCCTGTTc                                           <  7:253394/119‑1 (MQ=255)                                                                                                                               ccAGCAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGCCCTGGCCATCAGCGCCATGCTCTGGCGCGCACTGGCCTGCGGGGGAGTGGCAGCACTGGGTGCCGGCCTGTTc                                           >  8:253394/1‑119 (MQ=255)                                                                                                                                cAGCAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGCCCTGGCCATCAGCGCCATGCTCTGGCGCGCACTGGCCTGCGGGGGAGTGGCAGCACTGGGTGCCGGCCTGTTCgt                                         >  1:335731/1‑120 (MQ=255)                                                                                                                                cAGCAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGCCCTGGCCATCAGCGCCATGCTCTGGCGCGCACTGGCCTGCGGGGGAGTGGCAGCACTGGGTGCCGGCCTGTTCgt                                         <  2:335731/120‑1 (MQ=255)                                                                                                                                 aGCAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGCCCTGGCCATCAGCGCCATGCTCTGGCGCGCACTGGCCTGCGGGGGAGTGGCAGCACTGGGTGCCGGCCTGTTCGTGTTCTCCGACAGCCTGATCGGCATCGACCg           >  4:304316/1‑149 (MQ=255)                                                                                                                                 aGCAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGCCCTGGCCATCAGCGCCATGCTCTGGCGCGCACTGGCCAGCGGACGAGTGGCAGCACTGGGTGCCGGCCTGTTCGTGTTCTCCGACAGCCTGATCGGCATCGACCg           >  1:110691/1‑149 (MQ=255)                                                                                                                                   cAACGGCCTCGGCCCGTTGGTGATCCCGGTAGCACTCTACGCCCTGGCCATCAGCGCCATGCTCTGGCGCGCACTGGCCTGCGGGGGAGTGGCAGCACTGGGTGCCGGCCTGTTCGTGTTCTCCGACAGCCTGATCg                     >  6:251445/1‑137 (MQ=255)                                                                                                                                   cAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGCCCTGGCCATCAGCGCCATGCTCTGGCGCGCACTGGCCTGCGGGGGAGTGGCAGCACTGGGTGCCGGCCTGTTCGTGTTCTCCGACAGCCTGATCg                     <  5:251445/137‑1 (MQ=255)                                                                                                                                   cAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGCCCTGGCCATCAGCGCCATGCTCTGGCGCGCACTGGCCTGCGGGGGAGTGGCAGCACTGGGTGCCGACCTGTTCGTGTTCTCCGACAGCCTGATCGGCATCGACCGCt         <  1:309104/149‑1 (MQ=255)                                                                                                                                         ccTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGCCCTGGCCATCAGCGCCATGCTCTGGCGCGCACTGGCCTGCGGGGGAGTGGCAGCACTGGGTGCCGGCCTGTTCGTGTTCTCCGACAGCCTGATCGGCATCGACCGCTtcat     >  6:90405/1‑147 (MQ=255)                                                                                                                                          cTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGCCCTGGCCATCAGCGCCATGCTCTGGCGCGCACTGGCCTGCGGGGGAGTGGCAGCACTGGGTGCCGGCCTGTTCGTGTTCTCCGACAGCCTGATCGGCATCGACCGCTTCATCAg  >  4:333118/1‑149 (MQ=255)                                                                                                                                               |                                                                                                                                                CCTGGCCGCCTTCCTTTGCGCCCACCTTGCTTACCTGCGCGGCTATTGCGGCATTACCCTGCGCCCCGCCCTGCCC‑GCACTGATCTTCAGCGCAATCACTGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGCCCTGGCCATCAGCGCCATGCTCTGGCGCGCACTGGCCTGCGGGGGAGTGGCAGCACTGGGTGCCGGCCTGTTCGTGTTCTCCGACAGCCTGATCGGCATCGACCGCTTCATCAG  >  NC_002947/4655381‑4655665

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