Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,995,966	A→C	6.0%	F356V (TTC→GTC)	PP_5252 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,995,966	0	A	C	6.0%	206.4 / 2.4	65	F356V (TTC→GTC)	PP_5252	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (25/36);  new base C (2/2);  total (27/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

ACCTTCTCGAACACGTCGAGCATGCGCGTGATCGATTCGTTGTAGGTGGCGTGCAGGCGGAATGGCCAGCGCTGCTCGACCAGGTGGCGCACCACCGGCTCCAGTTCTTCTTCCATGGTCTGCGGCAGGTCCGGGCGCGGCTCGAGGAAGTCCTCGAAGTCGGCGGCAGAGAACA‑CCAGCATTTCGCCCGCGCCATTGT‑GGCGCAAGAAGTCGGTACCGCTGTGCAGCTTGACGCTGGAGGTCCACTTCTTGAAGTCGTCCAGTTCTTCCTTGGGCTTCTGGGTGAACAGGTTGTA  >  NC_002947/5995818‑5996113                                                                                                                                                     |                                                                                                                                                      aCCTTCTCGAACACGTCGAGCATGCGCGTGATCGATTCGTTGTAGGTGGCGTGCAGGCGGAATGGCCAGCGCTGCTCGACCAGGTGGCGCACCACCGGCTCCAGTTCTTCTTCCATGGTCTGCGGCAGGTCCGGGCGCGGCTCGAGGaa                                                                                                                                                       <  4:335083/149‑1 (MQ=255)  ccTTCTCGAACACGTCGAGCATGCGCGTGATCGATTCGTTGTAGGTGGCGTGCAGGCGGAATGGCCAGCGCTGCTCGACCAGGTGGCGCACCACCGGCTCCAGTTCTTCTTCCATGGTCTGCGGCAGGTCCGGGCGCGGCTCGAGGAAg                                                                                                                                                      <  4:24950/149‑1 (MQ=255)      ctcGAACACGTCGAGCATGCGCGTGATCGATTCGTTGTAGGTGGCGTGCAGGCGGAATGGCCAGCGCTGCTCGACCAGGTGGCGCACCACCGGCTCCAGTTCTTCTTCCATGGTCTGCGGCAGGTCCGGGCGCGGCTCGAGGAAGTcc                                                                                                                                                   <  6:158826/148‑1 (MQ=255)         gAACACGTCGAGCATGCGCGTGATCGATTCGTTGTAGGTGGCGTGCAGGCGGAATGGCCAGCGCTGCTCGACCAGGTGGCGCACCACCGGCTCCAGTTCTTCTTCCATGGTCTGCGGCAGGTCCGGGCGCGGCTCGAGGAAGTCCTc                                                                                                                                                 <  2:322171/147‑1 (MQ=255)         gAACACGTCGAGCATGCGCGTGATCGATTCGTTGTAGGTGGCGTGCAGGCGGAATGGCCAGCGCTGCTCGACCAGGTGGCGCACCAACGGCTCCAGTTCTTCTTCCATGGTCTGCGGCAGGTCCGGGCGCGGCTCGAGGAAGTCCTCg                                                                                                                                                >  8:86435/1‑148 (MQ=255)          aaCACGTCGAGCATGCGCGTGATCGATTCGTTGTAGGTGGCGTGCAGGCGGAATGGCCAGCGCTGCTCGACCAGGTGGCGCACCACCGGCTCCAGTTCTTCTTCCATGGTCTGCGGCAGGTCCGGGCGCGGCTCGAGGAAGTCCTCGaa                                                                                                                                              <  1:245876/149‑1 (MQ=255)           acacGTCGAGCATGCGCGTGATCGATTCGTTGTAGGTGGCGTGCAGGCGGAATGGCCAGCGCTGCTCGACCAGGTGGCGCACCACCGGCTCCAGTTCTTCTTCCATGGTCTGCGGCAGGTCCGGGCGCGGCTCGAGGAAGTCCTCGaat                                                                                                                                             >  3:297902/1‑148 (MQ=255)            cacGTCGAGCATGCGCGTGATCGATTCGTTGTAGGTGGCGTGCAGGCGGAATGGCCAGCGCTGCTCGACCAGGTGGCGCACCACCGGCTCCAGTTCTTCTTCCATGGTCTGCGGCAGGTCCGGGCGCGGCTCGAGGACTTCCTCGagcc                                                                                                                                            <  1:124228/149‑4 (MQ=255)            cacGTCGAGCATGCGCGTGATCGATTCGTTGTAGGTGGCGTGCAGGCGGAATGGCCAGCGCTGCTCGACCAGGTGGCGCACCACCGGCTCCAGTTCTTCTTCCATGGTCTGCGGCAGGTCCGGGCGCGGCTCGAGGAAGTCCTCGAAGt                                                                                                                                            <  6:303213/149‑1 (MQ=255)             acGTCGAGCATGCGCGTGATCGATTCGTTGTAGGTGGCGTGCAGGCGGAATGGCCAGCGCTGCTCGACCAGGTGGCGCACCACCGGCTCCAGTTCTTCTTCCATGGTCTGCGGCAGGTCCGGGCGCGGCTCGAGGAAGTCCTCGAAGTc                                                                                                                                           <  2:305836/149‑1 (MQ=255)                 cGAGCATGCGCGTGATCGATTCGTTGTAGGTGGCGTGCAGGCGGAATGGCCAGCGCTGCTCGACCAGGTGGCGCACCACCGGCTCCAGTTCTTCTTCCATGGTCTGCGGCAGGTCCGGGCGCGGCTCGAGGAAGTCCTCGAAGTcggcg                                                                                                                                       <  5:56918/149‑1 (MQ=255)                        gcgcGTGATCGATTCGTTGTAGGTGGCGTGCAGGCGGAATGGCCAGCGCTGCTCGACCAGGTGGCGCACCACCGGCTCCAGTTCTTCTTCCATGGTCTGCGGCAGGTCCGGGCGCGGCTCGAGGAAGTCCTCGAAGTCGGCGGCagag                                                                                                                                 >  3:262061/1‑148 (MQ=255)                         cgcgTGATCGATTCGTTGTAGGTGGCGTGCAGGCGGAATGGCCAGCGCTGCTCGACCAGGTGGCGCACCACCGGCTCCAGTTCTTCTTCCATGGTCTGCGGCAGGTCCGGGCGCGGCTCGAGGAAGTCCTCGAAGTCGGCGGCAGAGaa                                                                                                                               <  8:71418/149‑1 (MQ=255)                            gTGATCGATTCGTTGTAGGTGGCGTGCAGGCGGAATGGCCAGCGCTGCTCGACCAGGTGGCGCACCACCGGCTCCAGTTCTTCTTCCATGGTCTGCGGCAGGTCCGGGCGCGGCTCGAGGAAGTCCTCGAAGTCGGCGGCAGAGAaca‑c                                                                                                                           <  2:21893/149‑1 (MQ=255)                            gTGATCGATTCGTTGTAGGTGGCGTGCAGGCGGAAGGGCCAGCGCTGCTCGACCAGGTGGCGCACCACCGGCTCCAGTTCTTCTTCCATGGTCTGCGGCAGGTCCGGGCGCGGCTCGAGGAAGTCCTCGAAGTCGGCGGCAGAGAac                                                                                                                              <  6:161868/147‑1 (MQ=255)                               aTCGATTCGTTGTAGGTGGCGTGCAGGCGGAATGGCCAGCGCTGCTCGACCAGGTGGCGCACCACCGGCTCCAGTTCTTCTTCCATGGTCTGCGGCAGGTCCGGGCGCGGCTCGAGGAAGTCCTCGAAGTCGGCGGCAGAGAACA‑CCa                                                                                                                         <  1:14341/148‑1 (MQ=255)                                  gATTCGTTGTAGGTGGCGTGCAGGCGGAATGGCCAGCGCTGCTCGACCAGGTGGCGCACCACCGGCTCCAGTTCTTCTTCCATGGTCTGCGGCAGGTCCGGGCGCGGCTCGAGGAAGTCCTCGAAGTCGGCGGCAGAGAACA‑CCAGCAt                                                                                                                     >  3:190839/1‑149 (MQ=255)                                    ttCGTTGTAGGTGGCGTGCAGGCGGAATGGCCAGCGCTGCTCGACCAGGTGGCGCACCACCGGCTCCAGTTCTTCTTCCATGGTCTGCGGCAGGTCCGGGCGCGGCTCGAGGACTTCCTCGagc                                                                                                                                             <  6:121052/124‑3 (MQ=255)                                      cGTTGTAGGTGGCGTGCAGGCGGAATGGCCAGCGCTGCTCGACCAGGTGGCGCACCACCGGCTCCAGTTCTTCTTCCATGGTCTGCGGCAGGTCCGGGCGCGGCTCGAGGAAGTCCTCGAAGTCGGCGGCAGAGAACA‑CCAGCATTTc                                                                                                                  <  2:269364/148‑1 (MQ=255)                                      cGTTGTAGGTGGCGTGCAGGCGGAATGGCCAGCGCTGCTCGACCAGGTGGCGCACCACCGGCTCCAGTTCTTCTTCCATGGTCTGCGGCAGGTCCGGGCGCGGCTCGAGGAAGTCCTCGAAGTCGGCGGCAGAGAACA‑CCAGCATTTc                                                                                                                  <  1:221603/148‑1 (MQ=255)                                        ttGTAGGTGGCGTGCAGGCGGAATGGCCAGCGCTGCTCGACCAGGTGGCGCACCACCGGCTCCAGTTCTTCTTCCATGGTCTGCGGCAGGTCCGGGCGCGGCTCGAGGAAGTCCTCGAAGTCGGCGGCAGAGAACA‑CCAGCATTTCGcc                                                                                                               >  1:47314/1‑149 (MQ=255)                                          gTAGGTGGCGTGCAGGCGGAATGGCCAGCGCTGCTCGACCAGGTGGCGCACCACCGGCTCCAGTTCTTCTTCCATGGTCTGCGGCAGGTCCGGGCGCGGCTCGAGGAAGTCCTCGAAGTCGGCGGCAGAGAACA‑CCAGCATTTCGccc                                                                                                              <  4:190839/148‑1 (MQ=255)                                           tAGGTGGCGTGCAGGCGGAATGGCCAGCGCTGCTCGACCAGGTGGCGCACCACCGGCTCCAGTTCTTCTTCCATGGTCTGCGGCAGGTCCGGGCGCGGCTCGAGGAAGTCCTCGAAGTCGGCGGCAGAGAACA‑CCAGCATTTCGCCcgc                                                                                                            <  6:319036/149‑1 (MQ=255)                                              gTGGCGTGCAGGCGGAATGGCCAGCGCTGCTCGACCAGGTGGCGCACCACCGGCTCCAGTTCTTCTTCCATGGTCTGCGGCAGGTCCGGGCGCGGCTCGAGGAAGTCCTCGAAGTCGGCGGCAGAGAACA‑CCAGCATTTCGCCCgcgc                                                                                                          >  3:198794/1‑148 (MQ=255)                                                    tGCAGGCGGAATGGCCAGCGCTGCTCGACCAGGTGGCGCACCACCTGCTCCAGTTCTTCTTCCATGGTCTGCGACAGGTCCGGGCGCGGCTCGAGGAAGTCCTCGAAGTCGGCGGCAGAGAACA‑CCAGCATTTCGCCCGCGCCATtgt‑g                                                                                                  <  6:167994/149‑1 (MQ=255)                                                     gCAGGCGGAATGGCCAGCGCTGCTCGACCAGGTGGCGCACCACCGGCTCCAGTTCTTCTTCCATGGTCTGCGGCAGGTCCGGGCGCGGCTCGAGGAAGTCCTCGAAGTCGGCGGCAGAGAACA‑CCAGCATTTCGCCCGCGCCATTGT‑gg                                                                                                 <  7:138501/149‑1 (MQ=255)                                                       aGGCGGAATGGCCAGCGCTGCTCGACCAGGTGGCGCACCACCGGCTCCAGTTCTTCTTCCATGGTCTGCGGCAGGTCCGGGCGCGGCTCGAGGAAGTCCTCGAAGTCGGCGGCAGAGAACA‑CCAGCATTTCGCCCGCGCCATTGT‑Ggcg                                                                                               >  2:182603/1‑149 (MQ=255)                                                                         tgctCGACCAGGTGGCGCACCACCGGCTCCAGTTCTTCTTCCATGGTCTGCGGCAGGTACGGGCGCGGCTCGAGGAAGTCCTCGAAGTCGGCGGCAGAGAACA‑CCAGCATTTCGCCCGCGCCATTGT‑GGCGCAAGAAGTCGGTACCGCt                                                                             >  8:290244/1‑149 (MQ=255)                                                                           ctCGACCAGGTGGCGCACCACCGGCTCCAGTTCTTCTTCCATGGTCTGCGGCAGGTCCGGGCGCGGCTCGAGGAAGTCCTCGAAGTCGGCGGCAGAGAACA‑CCAGCATTTCGCCCGCGCCATTGT‑GGCGCAAGAAGTCGGTACCGCtg                                                                            <  6:239351/148‑1 (MQ=255)                                                                              gACCAGGTGGCGCACCACCGGCTCCAGTTCTTCTTCCATGGTCTGCGGCAGGTCCGGGCGCGGCTCGAGGAAGTCCTCGAAGTCGGCGGCAGAGAACA‑CCAGCATTTCGCCCGCGCCATTGTCGGCGCAAGAAGTCGGTACCGCTGTGc                                                                         <  1:4563/149‑1 (MQ=255)                                                                              gACCAGGTGGCGCACCACCGGCTCCAGTTCTTCTTCCATGGTCTGCGGCAGGTCCGGGCGCGGCTCGAGGAAGTCCTCGAAGTCGGCGGCAGAGAACA‑CCAGCATTTCGCCCGCGCCATTGT‑GGCGCAAGAAGTCGGTACCGCtgtg                                                                          >  8:39511/1‑147 (MQ=255)                                                                                ccAGGTGGCGCACCACCGGCTCCAGTTCTTCTTCCATGGTCTGCGGCAGGTCCGGGCGCGGCTCGAGGAAGTCCTCGAAGTCGGCGGCAGAGAACA‑CCAGCATTTCGCCCGCGCCATTGT‑GGCGCAAGAAGTCGGTACCGCTGTGCAGc                                                                      >  2:337610/1‑149 (MQ=255)                                                                                         gcACCACCGGCTCCAGTTCTTCTTCCATGGTCTGCGGCAGGTCCGGGCGCGGCTCGAGGACGTCCTCGAAGTCGGCGGCAGAGAACA‑CCAGCATTTCGCCCGCGCCATTGT‑GGCGCAAGCAGGCGGTACCGCTGTGCAGCTTGACGGTg                                                             >  4:40213/1‑149 (MQ=255)                                                                                         gcACCACCGGCTCCAGTTCTTCTTCCATGGTCTGCGGCAGGTCCGGGCGCGGCTCGAGGAAGTCCTCGAAGTCGGCGGCAGAGAACA‑CCAGCATTTCGCCCGCGCCATTGT‑GGCGCAAGAAGTCGGTACCGCTGTGCAGCTTGACGCTg                                                             >  3:302913/1‑149 (MQ=255)                                                                                         gcACCACCGGCTCCAGTTCTTCTTCCATGGTCTGCGGCAGGTCCGGGCGCGGCTCGAGGAAGTCCTCGAAGTCGGCGGCAGAGAACA‑CCAGCATTTCGCCCGCGCCATTGT‑GGCGCAAGAAGTCGGTACCGCTGTGCAGCTTGACGCTg                                                             >  4:35590/1‑149 (MQ=255)                                                                                          caccacCGGCTCCAGTTCTTCTTCCATGGTCTGCGGCAGGTCCGGGCGCGGCTCGAGGAAGTCCTCGAAGTCGGCGGCAGAGAACA‑CCAGCATTTCGCCCGCGCCATTGT‑GGCGCAAGAAGTCGGTACCGCTGTGCAGCTTGACGCTgg                                                            <  4:198794/149‑1 (MQ=255)                                                                                              accGGCTCCAGTTCTTCTTCCATGGTCTGCGGCAGGTCCGGGCGCGGCTCGAGGAAGTCCTCGAAGTCGGCGGCAGAGAACA‑CCAGCATTTCGCCCGCGCCATTGG‑GGCGCAAGAAGTCGGTACCGCTGTGCAGCTTGACGCTGGAGGt                                                        <  8:102324/149‑1 (MQ=255)                                                                                                      cAGTTCTTCTTCCATGGTCTGCGGCAGGTCCGGGCGCGGCTCGAGGAAGTCCTCGAAGTCGGCGGCAGAGAACA‑CCAGCATTTCGCCCGCGCCATTGT‑GGCGCAAGAAGTCGGTACCGCTGTGCAGCTTGACGCTGGAGGTCCActtct                                                >  6:276022/1‑149 (MQ=255)                                                                                                       aGTTCTTCTTCCATGGTCTGCGGCAGGTCCGGGCGCGGCTCGAGGAAGTCCTCGAAGTCGGCGGCAGAGAACA‑CCAGCATTTCGCCCGCGCCATTGT‑GGCGCAAGAAGTCGGTACCGCTGTGCAGCTTGACGCTGGAGGTCCActtctt                                               <  7:261739/149‑1 (MQ=255)                                                                                                        gTTCTTCTTCCATGGTCTGCGGCAGGTCCGGGCGCGGCTCGAGGAAGTCCTCGAAGTCGGCGGCAGAGAACA‑CCAGCATTTCGCCCGCGCCATTGT‑GGCGCAAGAAGTCGGTACCGCTGTGCAGCTTGACGCTGGAGGTCCAct                                                   >  3:91010/1‑144 (MQ=255)                                                                                                        gTTCTTCTTCCATGGTCTGCGGCAGGTCCGGGCGCGGCTCGAGGAAGTCCTCGAAGTCGGCGGCAGAGAACA‑CCAGCATTTCGCCCGCGCCATTGT‑GGCGCAAGAAGTCGGTACCGCTGTGCAGCTTGACGCTGGAGGTCCAct                                                   <  4:91010/144‑1 (MQ=255)                                                                                                        gTTCTTCTTCCATGGTCTGCGGCAGGTCCGGGCGCGGCTCGAGGAAGTCCTCGAACTCGGCCGCAGAGAACA‑ACAGCATTTAGCCCGCGCCATAGT‑GGCGCAAGAAGACGGAACCGCAGTGCAGCTTGAAGCTGGAGGTCCActtctt                                               >  4:128096/1‑148 (MQ=255)                                                                                                          tcttcttcCATGGTCTGCGGCAGGTCCGGGCGCGGCTCGAGGAAGTCCTCGAAGTCGGCGGCAGAGAACA‑CGAGCATTTCGCCCGCGCCATTGT‑GGCGCAAGAAGTCGGTACCGCTGTGCAGCTTGACGCTGGAGGTCCACTTCTTg                                              >  6:106257/1‑147 (MQ=255)                                                                                                           cttcttcCATGGTCTGCGGCAGGTCCGGGCGCGGCTCGAGGAAGTCCTCGAAGTCGGCTGCAGAGAACA‑CCAGCATTTCGCCAGCGCCATTGT‑GGCTCAAGAAGTCGGTACCGCTGTGCAGCTTGACGCGGGAGGTCCACTTCTTGaa                                            <  7:86435/148‑1 (MQ=255)                                                                                                               ttcCATGGTCTGCGGCAGGTCCGGGCGCGGCTCGAGGAAGTCCTCGAAGTCGGCGGCAGAGAACA‑CCAGCATTTCGCCCGCGCCATTGT‑GGCGCAAGAAGTCGGTACCGCTGTGCAGCTTGACGCTGGAGGTCCACTTCTTGAAgtcgt                                       >  2:137768/1‑149 (MQ=255)                                                                                                                tcCATGGTCTGCGGCAGGTCCGGGCGCGGCTCGAGGAAGTCCTCGAAGTCGGTGGCAGAGAACA‑CCAGCATTTCGCCCGCGCCATTGT‑GGCGCAAGAAGTCGGTACCGCTGTGCAGCTTGACGCTGGAGGTCCACTTCTTGAAgtcgtc                                      >  4:311597/1‑149 (MQ=255)                                                                                                                 ccATGGTCTGCGGCAGGTCCGGGCGCGGCTCGAGGAAGTCCTCGAAGTCGGCGGCAGAGAACA‑CCAGCATTTCGCCCGCGCCATTGT‑GGCGCAAGAAGTCGGTACCGCTGTGCAGCTTGACGCTGGAGGTCCACTTCTTGAAgtcgtc                                      <  4:262061/148‑1 (MQ=255)                                                                                                                 ccATGGTCTGCGGCAGGTCCGGGCGCGGCTCGAGGAAGTCCTCGAAGTCGGCGGCAGAGAACA‑CCAGCATTTCGCCCGCGCCATTGT‑GGCGCAAGAAGTCGGTACCGCTGTGCAGCTTGACGCTGGAGGTCCACTTCTTGAAgtcgt                                       <  7:290244/147‑1 (MQ=255)                                                                                                                 acATGTTCAGCGGCAGGTCCGGGCGCGGCTCGAGGAAGTCCTCGAAGTCGGCGGCAGAGAACA‑CCAGCATTTCGCCCGCGCCATTGT‑GGCGCAAGAAGTCGGTACCGCTGTGCAGCTTGACGCTGGAGGTCCACTTCTTGAAgtcgtc                                      >  8:108219/2‑148 (MQ=255)                                                                                                                  cATGGTCTGCGGCAGGTCCGGGCGCGGCTCGAGGAAGTCCTCGAAGTCGGCGGCAGAGAACA‑CCAGCATTTCGCCCGCGCCATTGT‑GGCGCAAGAAGTCGGTACCGCTGTGCAGCTTGACGCTGGAGGTCCACTTCTTGAAGTCGTCCa                                    <  4:217233/149‑1 (MQ=255)                                                                                                                     ggTCTGCGGCAGGTCCGGGCGCGGCTCGAGGAAGTCCTCGAAGTCGGCGGCAGAGAACA‑CCAGCATTTCGCCCGCGCCATTGT‑GGCGCAAGAAGTCGGTACCGCTGTGCAGCTTGACGCTGGAGGTCCACTTCTTGAAGTCGTCCAGt                                  >  6:224598/1‑148 (MQ=255)                                                                                                                      gTCTGCGGCAGGTCCGGGCGCGGCTCGAGGAAGTCCTCGAAGTCGGCGGCAGAGAACA‑CCAGCATTTCGCCCGCGCCATTGT‑GGCGCAAGAAGTCGGTACCGCTGTGCAGCTTGACGCTGGAGGTCCACTTCTTGAAGTCGTCCAGtt                                 <  1:337610/148‑1 (MQ=255)                                                                                                                         tGCGGCAGGTCCGGGCGCGGCTCGAGGAAGTCCTCGAAGTCGGCGGCAGAGAACA‑CCAGCATTTCGCCCGCGCCATTGT‑GGCGCAAGAAGTCGGTACCGCTGTGCAGCTTGACGCTGGAGGTCCACTTCTTGAAGTCGTCCAGttcttc                             <  5:224598/149‑1 (MQ=255)                                                                                                                         tGCGGCAGGTCCGGGCGCGGCTCGAGGAAGTCCTCGAAGTCGGCGGCAGAGAACA‑CCAGCATTTCGCCCGCGCCATTGT‑GGCGCAAGAAGTCGGTACCGCTGTGCAGCATGACGCTGGAGGTCCACTTCTTGa                                             >  3:76591/1‑133 (MQ=255)                                                                                                                           cGGCAGGTCCGGGCGCGGCTCGAGGAAGTCCTCGAAGTCGGCGGCAGAGAACATCCAGCATTTCGCCCGCGCCATTGT‑GGCGCAAGAAGTCGGTACCGCTGTGCAGCTTGACGCTGGAGGTCCACTTCTTGAAGTCGTCCAGTTCTTcc                            >  1:288683/1‑149 (MQ=255)                                                                                                                            ggCAGGTCCGGGCGCGGCTCGATGAAGTCCTCGAAGTCGGCGGCAGAGAACA‑CCAGCATTTCGCCCGCGCCATTGT‑GGCGCAAGAAGTCGGTACCGCTGTGCAGCTTGACGCTGGAGGTCCACTTCTTGAAGTCGTCCAGTTCTTCCtt                          <  3:40213/149‑1 (MQ=255)                                                                                                                            ggCAGGTCCGGGCGCGGCTCGAGGAAGTCCTCGAAGTCGGCGGCAGAGAACA‑CCAGCATTTCGCCCGCGCCATTGT‑GGCGCAAGAAGTCGGTACCGCTGTGCAGCTTGACGCTGGAGGTCCACTTCTTGAAGTCGTCCAGTTCTTCCtt                          <  3:35590/149‑1 (MQ=255)                                                                                                                             gCAGGTCCGGGCGCGGCTCGAGGAAGTCCTCGAAGTCGGCGGCAGAGAACA‑CCAGCATTTCGCCCGCGCCATTGT‑GGCGCAAGAAGTCGGTACCGCTGTGCAGCTTGACGCTGGAGGTCCACTTCTTGAAGTCGTCCAGTTCTTCCTTg                         <  2:47314/149‑1 (MQ=255)                                                                                                                                  tCCGGGCGCGGCTCGAGGAAGTCCTCGAAGTCGGCGGCAGAGAACA‑CCAGCATTTCGCCCGCGCCATTGT‑GGCGCAAGAAGTCGGTACCGCTGTGCAGCTTGACGCTGGAGGTCCACTTCTTGAAGTCGTCCAGTTCTTCCTTGGGCtt                    >  4:176325/1‑149 (MQ=255)                                                                                                                                           ggCTCGAGGAAGTCCTCGAAGTCGGCGGCAGAGAACA‑CCAGCATTTCGCCCGCGCCATTGT‑GGCGCAAGAAGTCGGTACCGCTGTGCAGCTTGACGCTGGAGGTCCACTTCTTGAAGTCGTCCAGTTCTTCCTTGGGCTTCTGGGTGaa           <  6:72166/149‑1 (MQ=255)                                                                                                                                           acttCGAGGACTTCCTCGAAGTCGGCGGCAGAGAACA‑CCAGCATTTCGCCCGCGCCATTGT‑GGCGCAAGAAGTCGGTACCGCTGTGCAGCTTGACGCTGGAGGTCCACTTCTTGAAGTCGTCCAGTTCTTCCTTGGGCTTCTGGGTGaa           >  3:177796/4‑149 (MQ=255)                                                                                                                                                gAGGAAGTCCTCGAAGTCTGCGGCAGAGAACA‑CCATCATTTCGCCCTCGCCATTGT‑GGCGCAAGATGTCGGTACCGCTGTGCAGCTTGACGCTGGAGGTCCACTTCTTGAAGTCGTCGATTTCTTCCTTTGGCTTCTGGGTGAACAGGt      <  4:297902/149‑1 (MQ=255)                                                                                                                                                gAGGAAGTCCTCGAAGTCGGCGGCAGAGAACA‑CCAGCATTTCGCCCGCGCCATTGT‑GGCGCAAGAAGTCGGTACCGCTGTGCAGCTTGACGCTGGAGGTCCACTTCTTGAAGTCGTCCAGTTCTTCCTTGGGCTTCTGGGTGAACAGGt      >  5:165448/1‑149 (MQ=255)                                                                                                                                                 aGGAAGTCCTCGAAGTCGGCGGCAGAGAACA‑CCAGCATTTCGCCCGCGCCATTGT‑GGCGCAAGAAGTCGGTACCGCTGTGCAGCTTGACGCTGGAGGTCCACTTCTTGAAGTCGTCCAGTTCTTCCTTGGGCTTCTGGGTGAACAGGtt     <  6:165448/149‑1 (MQ=255)                                                                                                                                                  ggAAGTCCTCGAAGTCGGCGGCAGAGAACA‑CCAGCATTTCGCCCGCGCCATTGT‑GGCGCAAGAAGTCGGTACCGCTGTGCAGCTTGACGCTGGAGGTCCACTTCTTGAAGTCGTCCAGTTCTTCCTTGGGCTTCTGGGTGAACAGGTTg    >  1:86385/1‑149 (MQ=255)                                                                                                                                                   gAAGTCCTCGAAGTCGGCGGCAGAGAACA‑CCAGCATTTCGCCCGCGCCATTGT‑GGCGCAAGAAGTCGGTACCGCTGTGCAGCTTGACGCTGGAGGTCCACTTCTTGAAGTCGTCCAGTTCTTCCTTGGGCTTCTGGGTGAACAGGTTGt   <  2:288683/149‑1 (MQ=255)                                                                                                                                                    aaGTCCTCGAAGTCGGCGGCAGAGAACA‑CCAGCATTTCGCCCGCGCCATTGT‑GGCGCAAGAAGTCGGTACCGCTGTGCAGCTTGACGCTGGAGGTCCACTTCTTGAAGTCGTCCAGTTCTTCCTTGGGCTTCTGGGTGAACAGGTTGTa  <  7:166102/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                      ACCTTCTCGAACACGTCGAGCATGCGCGTGATCGATTCGTTGTAGGTGGCGTGCAGGCGGAATGGCCAGCGCTGCTCGACCAGGTGGCGCACCACCGGCTCCAGTTCTTCTTCCATGGTCTGCGGCAGGTCCGGGCGCGGCTCGAGGAAGTCCTCGAAGTCGGCGGCAGAGAACA‑CCAGCATTTCGCCCGCGCCATTGT‑GGCGCAAGAAGTCGGTACCGCTGTGCAGCTTGACGCTGGAGGTCCACTTCTTGAAGTCGTCCAGTTCTTCCTTGGGCTTCTGGGTGAACAGGTTGTA  >  NC_002947/5995818‑5996113

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