Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,103,967	G→T	5.5%	R267L (CGC→CTC)	PP_5353 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	6,103,967	0	G	T	5.5%	295.1 / 5.7	91	R267L (CGC→CTC)	PP_5353	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (43/43);  new base T (3/2);  total (46/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.91e-01

TACATCGATAACCTGGAAGTGCAACTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAAT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATC  >  NC_002947/6103825‑6104115                                                                                                                                                 |                                                                                                                                                     tacATCGATAACCTGGAAGTGCAACTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAAT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGccc                                                                                                                                                 >  3:241805/1‑148 (MQ=255) tacATCGATAACCTGGAAGTGCAACTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAAT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGccc                                                                                                                                                 >  6:310025/1‑148 (MQ=255)  acATCGATAACCTGGAAGTGCAACTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAAT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGccc                                                                                                                                                 >  3:102833/1‑147 (MQ=255)    aTCGATAACCTGGAAGTGCAACTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAAT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCaatg                                                                                                                                             <  2:250075/149‑3 (MQ=255)    aTCGATAACCTGGAAGTGCAACTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAAT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAgg                                                                                                                                             >  5:249035/1‑149 (MQ=255)     tCGATAACCTGGAAGTGCAACTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAATAACGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAgg                                                                                                                                             <  1:298302/149‑1 (MQ=255)       gATAACCTGGAAGTGCAACTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAAT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCga                                                                                                                                          <  4:59204/149‑1 (MQ=255)       gATAACCTGGAAGTGCAACTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAAT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGCTCGCGCCCAAGGCga                                                                                                                                          <  7:4900/149‑1 (MQ=255)            ccTGGAAGTGCAACTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAAT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAAc                                                                                                                                      >  4:59887/1‑148 (MQ=255)             cTGGAAGTGCAACTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAAT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAAc                                                                                                                                      >  5:243608/1‑147 (MQ=255)                    tGCAACTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAAT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCTCGCCTTGGGCGCGAACCatgtgcca                                                                                                                             <  8:336343/149‑8 (MQ=255)                    tGCAACTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAAT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTgg                                                                                                                             <  4:241805/149‑1 (MQ=255)                     gCAACTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAAT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCTCGCCTTGGGCGCGAACCatgtgcca                                                                                                                             <  4:307558/148‑8 (MQ=255)                      cAACTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAAT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGc                                                                                                                            <  4:354260/148‑1 (MQ=255)                      cAAATGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTTTCGAAGATT‑TCGCGCTTCGACGAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATTGTTCGCGGCCAAGGCGAGAACCAGCCTTGGc                                                                                                                            <  3:122498/148‑1 (MQ=255)                       aaCTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAAT‑TCGCGCTGCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTg                                                                                                                          <  3:34012/149‑1 (MQ=255)                        aCTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAAT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTgg                                                                                                                         <  6:249035/149‑1 (MQ=255)                        aCTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAAT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTgg                                                                                                                         >  5:247633/1‑149 (MQ=255)                              cGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAAT‑TCTCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCGCAAGTCGAGAACCAGCCTTGGCTGGTAGCCg                                                                                                                   <  8:149030/149‑1 (MQ=255)                              cGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAAT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCg                                                                                                                   >  3:181409/1‑149 (MQ=255)                               ggTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAAT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCgg                                                                                                                  <  5:251697/149‑1 (MQ=255)                                   tCGATGACCGCGCCGAGCGGACCGACGCTACCCTGACCTTCCGGGGCGTGTCGAAGAAT‑TCGCGCTTCGACCAGGGCTAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCTCCCAAGGCGAGAACCAGCGTTGGCTGGTAGCCGGTAt                                                                                                               <  4:324506/148‑1 (MQ=255)                                     gATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAAT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTgg                                                                                                                         >  1:301422/1‑136 (MQ=255)                                     gATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAAT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTgg                                                                                                                         <  2:301422/136‑1 (MQ=255)                                     gATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAAT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCg                                                                                                            >  4:185334/1‑149 (MQ=255)                                          ccGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAAT‑TCGCGCTTCGACCACGGCGAA‑GTGTTCAGCGAAACCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCAGTCaa                                                                                                        >  6:174235/1‑148 (MQ=255)                                              gcCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAAT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAg                                                                                                   <  3:132764/149‑1 (MQ=255)                                              gcCCACCGCACCTAACGTACCCTGTCCTTACGTGGGGTTTCGAATAAT‑TGGGGGTTTGCTCATGGGGAATGTCTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACa                                                                                                    <  1:117811/149‑1 (MQ=255)                                                      cACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAAT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGc                                                                                           >  7:127943/1‑149 (MQ=255)                                                      cACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAAT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGc                                                                                           >  3:296732/1‑149 (MQ=255)                                                         cGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAAT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGc                                                                                        >  5:133575/1‑149 (MQ=255)                                                         cGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAAT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGc                                                                                        >  5:328896/1‑149 (MQ=255)                                                            cGCTACCCTGACCTTCCGTGGCGTGTCGAAGAAT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTg                                                                                     >  5:70200/1‑149 (MQ=255)                                                                 cccTGACCTTCCGTGGCGTGTCGAAGAAT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTaa                                                                                 <  4:287399/148‑1 (MQ=255)                                                                 cccTGACCTTCCGTGGCGTGTCGAAGAAT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTaa                                                                                 <  2:346478/148‑1 (MQ=255)                                                                  ccTGACCTTCCGTGGCGTGTCGAAGAAT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTaaa                                                                                >  2:166267/1‑148 (MQ=255)                                                                  ccTGACCTTCCGTGGCGTGTCGAAGAAT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTaaa                                                                                >  7:345296/1‑148 (MQ=255)                                                                    tGACCTTCCGTGGCGTGTCGAAGAAT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCGTTGTCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGc                                                                                        <  3:351565/138‑1 (MQ=255)                                                                    tGACCTTCCGTGGCGTGTCGAAGAAT‑TCGCGCTTCGACAAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTCGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCgacaacgg                                                                                >  4:351565/1‑139 (MQ=255)                                                                       ccTTCCGTGGCGTGTCGAAGAAT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAg                                                                           >  5:342965/1‑148 (MQ=255)                                                                       ccTTCCGTGGCGTGTCGAAGAAT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAg                                                                           >  7:346079/1‑148 (MQ=255)                                                                         ttCCGTGGCGTGTCGAAGAAT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAAc                                                                        >  3:63542/1‑149 (MQ=255)                                                                         ttCCGTGGCGTGTCGAAGAAT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAAc                                                                        >  4:276051/1‑149 (MQ=255)                                                                            cGTGGCGTGTCGAAGAAT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAAcccc                                                                     <  6:247633/149‑1 (MQ=255)                                                                                gCGTGTCGAAGAAT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAATCAGCTAATCGCTGCGTGCTAAACACAGAACCCCggg                                                                  <  6:243608/148‑1 (MQ=255)                                                                                gCGTGTCGAAGAAT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCggg                                                                  <  8:132991/148‑1 (MQ=255)                                                                                gCGTGTCGAAGAAT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCggg                                                                  <  2:254878/148‑1 (MQ=255)                                                                                gCGTGTCGAAGAAT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTCGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCACGGGc                                                                 >  6:264426/1‑149 (MQ=255)                                                                                     tCGAAGAAT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCg                                                                                                                   <  4:95013/94‑1 (MQ=255)                                                                                     tCGAAGAAT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCg                                                                                                                   >  3:95013/1‑94 (MQ=255)                                                                                     tCGAAGAAT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGcc                                                            >  1:259659/1‑149 (MQ=255)                                                                                       gaagaaT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCg                                                          <  4:102833/149‑1 (MQ=255)                                                                                        aagaaT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCgg                                                         >  6:219508/1‑149 (MQ=255)                                                                                        aagaaT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCgg                                                         >  6:190046/1‑149 (MQ=255)                                                                                        aagaaT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCAGGTATCCGTCAAAACAGCTAATCGCAGCGTGCTAAACACAGAACCCCGGGCTTGCCCgg                                                         >  6:346944/1‑149 (MQ=255)                                                                                        aagaaT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCACCGAAAGCTGGCACATGGTAAGCGCCCAAGGGGAGAACCAGCCTTGGCTGGTAGCCGCTATCCGTCAAAACAGCTAATCGCTGCCTGCTAAACACAGACCCACGGGCTTGCCCgg                                                         >  4:31648/1‑149 (MQ=255)                                                                                         agaaT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCggg                                                        <  6:342965/149‑1 (MQ=255)                                                                                         agaaT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCggg                                                        <  8:346079/149‑1 (MQ=255)                                                                                         agaaT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCggg                                                        >  1:144031/1‑149 (MQ=255)                                                                                         agaaT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCAAGCCCggg                                                        <  4:296732/149‑1 (MQ=255)                                                                                          gaaT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCgggg                                                       >  2:95508/1‑149 (MQ=255)                                                                                           aaT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGt                                                      <  6:70200/149‑1 (MQ=255)                                                                                                 gcgcTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAAAACCCCGGGCAAGCCCGGGGttctgt                                                 <  1:145511/149‑5 (MQ=255)                                                                                                    cTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTttttgctt                                              >  5:161083/1‑149 (MQ=255)                                                                                                       cGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTttttgcttt                                             <  5:190046/147‑1 (MQ=255)                                                                                                        gACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTttttgcttttg                                           <  4:181409/148‑1 (MQ=255)                                                                                                                cGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTATTCTTTTGCCAAAGTgg                                  >  7:7162/1‑149 (MQ=255)                                                                                                                  aa‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAac                                                                               >  5:94509/1‑102 (MQ=255)                                                                                                                  aa‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAac                                                                               <  6:94509/102‑1 (MQ=255)                                                                                                                       gtTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGctac                             >  7:17414/1‑148 (MQ=255)                                                                                                                       gtTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGctac                             >  5:213418/1‑148 (MQ=255)                                                                                                                        ttCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGctac                             >  6:14685/1‑147 (MQ=255)                                                                                                                           aGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAggg                        <  3:276051/149‑1 (MQ=255)                                                                                                                           aGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAgg                         <  6:180329/148‑1 (MQ=255)                                                                                                                           aGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAgg                         <  2:259659/148‑1 (MQ=255)                                                                                                                           aGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAgg                         <  5:219508/148‑1 (MQ=255)                                                                                                                           aGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAgg                         <  8:345296/148‑1 (MQ=255)                                                                                                                               aaaGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAgg                         >  2:34280/1‑144 (MQ=255)                                                                                                                               aaaGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGCTTTTTGCTTTTGCCAAAGTGGGCTACTAgg                         <  1:34280/144‑1 (MQ=255)                                                                                                                                 aGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAg                  <  5:310025/149‑1 (MQ=255)                                                                                                                                  gCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAgg                 <  3:237052/149‑1 (MQ=255)                                                                                                                                   cTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGgc                >  8:160733/1‑149 (MQ=255)                                                                                                                                      aaggctGGTTCTCGCCTTGGGCGCGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGgcg               >  2:372894/6‑147 (MQ=255)                                                                                                                                       aggctGGTTCTCGCCTTGGGCGCGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACACCGGGCTTGGCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAg            >  5:242158/5‑149 (MQ=255)                                                                                                                                        acaTGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGc           >  2:127263/1‑149 (MQ=255)                                                                                                                                        acaTGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGc           <  8:17414/149‑1 (MQ=255)                                                                                                                                           tGGTTCTCGCCTTGGGCGCGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGtt        >  8:99709/1‑149 (MQ=255)                                                                                                                                             gTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGt      >  6:71724/1‑149 (MQ=255)                                                                                                                                             gTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGt      <  5:224881/149‑1 (MQ=255)                                                                                                                                              ttCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGt      >  7:128062/1‑148 (MQ=255)                                                                                                                                               tCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCTGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTc     <  6:328896/148‑1 (MQ=255)                                                                                                                                                 gcgcCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATc  <  6:213418/149‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                     TACATCGATAACCTGGAAGTGCAACTGGACGGTGTCGATGACCGCGCCGACCGCACCGACGCTACCCTGACCTTCCGTGGCGTGTCGAAGAAT‑TCGCGCTTCGACCAGGGCGAA‑GTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATC  >  NC_002947/6103825‑6104115

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑