Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,655,520	G→A	7.8%	G139S (GGC→AGC) ‡	PP_4118 →	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,655,520	0	G	A	7.8%	179.6 / 4.2	64	G139S (GGC→AGC) ‡	PP_4118	membrane protein
Reads supporting (aligned to +/- strand):  ref base G (31/28);  new base A (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 = 1.91e-01

GTGTTCGGCCTGGCCGCCTTCCTTTGCGCCCACCTTGCTTACCTGCGCGGCTATTGCGGCATTACCCTGCGCCCCGCCCTGCCC‑GCACTGATCTTCAGCGCAATCACTGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGCCCTGGCCATCAGCGCCATGCTCTGGCGCGCACTGGCCTGCGGGGGAGTGGCAGCACTGGGTGCCGGCCTGTTCGTGTTCTCCGACAGCCTGATCGGCATCGACCGCTTCATCAG  >  NC_002947/4655373‑4655665                                                                                                                                                     |                                                                                                                                                  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                                                                                                                                       <  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)                                                                                                                                                     |                                                                                                                                                  GTGTTCGGCCTGGCCGCCTTCCTTTGCGCCCACCTTGCTTACCTGCGCGGCTATTGCGGCATTACCCTGCGCCCCGCCCTGCCC‑GCACTGATCTTCAGCGCAATCACTGGCATCGCGCTGCTCGGCGTACTGGCCAGCAACGGCCTCGGCCCGTTGCTGATCCCGGTAGCACTCTACGCCCTGGCCATCAGCGCCATGCTCTGGCGCGCACTGGCCTGCGGGGGAGTGGCAGCACTGGGTGCCGGCCTGTTCGTGTTCTCCGACAGCCTGATCGGCATCGACCGCTTCATCAG  >  NC_002947/4655373‑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: ‑