Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,655,516	C→G	7.8%	G137G (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,516	0	C	G	7.8%	196.2 / 6.3	64	G137G (GGC→GGG)	PP_4118	membrane protein
Reads supporting (aligned to +/- strand):  ref base C (30/29);  new base G (2/3);  total (32/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.44e-01

TTCGTGTTCGGCCTGGCCGCCTTCCTTTGCGCCCACCTTGCTTACCTGCGCGGCTATTGCGGCATTACCCTGCGCCCCGCCCTGCCC‑GCACTGATCTTCAGCGCAATCACTGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGCCCTGGCCATCAGCGCCATGCTCTGGCGCGCACTGGCCTGCGGGGGAGTGGCAGCACTGGGTGCCGGCCTGTTCGTGTTCTCCGACAGCCTGATCGGCATCGACCGCTTCAT  >  NC_002947/4655370‑4655662                                                                                                                                                    |                                                                                                                                                   ttcgtgttcgGCCTGGCCGCCTTCCTTTGCGCCCACCTTGCTTACCTGCGCGGCTATTGCGGCATTACCCTGCGCCCCGCCCTGCCC‑GCACTGATCTTCAGCGCAATCACTGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCt                                                                                                                                                  >  4:58042/1‑149 (MQ=255)    gtgttcgGCCTGGCCGCCTTCCTTTGCGCCCACCTTGCTTACCTGCGCGGCTATTGCGGCATTACCCTGCGCCCCGCCCTGCCC‑GCACTGATCTTCAGCGCAATCACTGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCTCgg                                                                                                                                               <  4:81423/149‑1 (MQ=255)            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                                                                                                                                    <  5:134962/148‑1 (MQ=255)                gccgccTTCCTTTGCGCCCACCTTGCTTACCTGCGCGGCTATTGCGGCATTACCCTGCGCCCCGCCCTGCCC‑GCACTGATCTTCAGCGCAATCACTGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCTCGGCCCGTTGCTGa                                                                                                                                    <  4:162913/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                                                                                                                                     >  2:59605/1‑99 (MQ=38)                                                                aTTACCCTGCGCCCCGCCCTGCCC‑GCACTGATCTTCAGCGCAATCACTGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGGCCGAGGCCGTTGCTg                                                                                                                                     <  1:59605/99‑1 (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)                                                                                                                                                    |                                                                                                                                                   TTCGTGTTCGGCCTGGCCGCCTTCCTTTGCGCCCACCTTGCTTACCTGCGCGGCTATTGCGGCATTACCCTGCGCCCCGCCCTGCCC‑GCACTGATCTTCAGCGCAATCACTGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGCCCTGGCCATCAGCGCCATGCTCTGGCGCGCACTGGCCTGCGGGGGAGTGGCAGCACTGGGTGCCGGCCTGTTCGTGTTCTCCGACAGCCTGATCGGCATCGACCGCTTCAT  >  NC_002947/4655370‑4655662

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