Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	806,328	G→A	7.0%	L334L (CTC→CTT)	PP_0693 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	806,328	0	G	A	7.0%	168.3 / 2.6	57	L334L (CTC→CTT)	PP_0693	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (27/26);  new base A (2/2);  total (29/28)
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

AGGTCGCGTACGCGCTCGTCGAGGTCTTCGGCACTGAACTGCCGTGACTGGGCCAGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGTTCTGCAGGCGGTCCTGCTCGCCACGGTAGTGCTCGGACTGGATGACCAGCTCTTCCTT  >  NC_002947/806182‑806474                                                                                                                                                   |                                                                                                                                                           aGGTCGCGTACGCGCTCGTCGAGGTCTTCGGCACTGAACTGCCGTGACTGGGCCAGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAg                                                                                                                                                          <  1:127045/149‑1 (MQ=255) aGGTCGCGTACGCGCTCGTCGACGTCTTCGGCACTGAACTGCCGTGACTGGGCCAGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAg                                                                                                                                                          <  2:242558/149‑1 (MQ=255)     cgcgTACGCGCTCGTCGAGGTCTTCGGCACTGAACTGCCGTGACTGGGCCAGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAG‑‑‑CAc‑‑‑‑c                                                                                                                                               <  5:153863/149‑1 (MQ=255)      gcgTACGCGCTCGTCGAGGTCTTCGGCACTGAACTGCCGTGACTGGGCCAGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAG‑‑‑CAc‑‑‑‑c                                                                                                                                               >  7:221325/1‑148 (MQ=255)      gcgTACGCGCTCGTCGAGGTCTTCGGCACTGAACTGCCGTGACTGGGCCAGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAG‑‑‑CAc‑‑‑‑c                                                                                                                                               <  4:40996/148‑1 (MQ=255)        gTACGCGCTCGTCGAGGTCTTCGGCACTGAACTGCCGTGACTGGGCCAGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCAGCTCCAGTACCTTGCCGTCATCGACCAGCGCAAAGCGGTGCTTGAG‑‑‑CAC‑‑‑‑CGa                                                                                                                                             >  8:90576/1‑148 (MQ=255)           cgcgCTCGTCGAGGTCTTCGGCACTGAACTGCCGTGACTGGGCCAGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTCAAG‑‑‑CAC‑‑‑‑Cgctttg                                                                                                                                         <  4:49326/149‑6 (MQ=255)              gcTCGTCGAGGTCTTCGGCACTGAACTGCCGTGACTGGGCCAGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCg                                                                                                                                     >  5:268609/1‑149 (MQ=255)                  gtcgAGGTCTTCGGCACTGAACTGCCGTGACTGGGCCAGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGcc                                                                                                                                   <  3:103815/147‑1 (MQ=255)                     gAGGTCTTCGGCACTGAACTGCCGTGACTGGGCCAGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGc                                                                                                                               <  6:268609/148‑1 (MQ=255)                          cTTCGGCACTGAACTGCCGTGACTGGGCCAGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAgcgc                                                                                                                         >  2:234577/1‑149 (MQ=255)                             cGGCACTGAACTGCCGTGACTGGGCCAGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCt                                                                                                                        <  1:234577/147‑1 (MQ=255)                                          ccGTGACTGGGCCAGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCTGGATACCGGTAAtt                                                                                                          >  4:32252/1‑148 (MQ=255)                                            gTGACTGGGCCAGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTcc                                                                                                       >  6:69336/1‑148 (MQ=255)                                                 tGGGCCAGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCTg                                                                                                                       <  5:267136/128‑1 (MQ=255)                                                 tGGGCCAGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCTg                                                                                                                       >  6:267136/1‑128 (MQ=255)                                                   ggCAAGCGCGCCGGCCAGCGCGTCGTTGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCa                                                                                                <  4:159522/149‑1 (MQ=255)                                                       aGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGc                                                                                                                               >  4:252083/1‑114 (MQ=255)                                                       aGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGgcg                                                                                            <  6:51464/149‑1 (MQ=255)                                                          gcgcCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCAt                                                                                         <  1:5985/149‑1 (MQ=255)                                                                 ccAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTCGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTcc                                                                                  >  3:211085/1‑149 (MQ=255)                                                                    gCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTgg                                                                               <  8:221325/149‑1 (MQ=255)                                                                         tcgtGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAGGTAAACGGGAGGATGAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTcc                                                                                  >  1:86821/1‑149 (MQ=255)                                                                         tcgtGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTgcc                                                                          <  1:31468/149‑1 (MQ=255)                                                                            tGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTgccgc                                                                        >  4:26310/1‑148 (MQ=255)                                                                            tGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTgccgc                                                                        >  4:67684/1‑148 (MQ=255)                                                                              ggCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCa                                                                      <  3:32252/148‑1 (MQ=255)                                                                                   cTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTgg                                                                >  5:42660/1‑149 (MQ=255)                                                                                    tttGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGt                                                               >  7:25608/1‑149 (MQ=255)                                                                                     ttGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTc                                                              >  1:134998/1‑149 (MQ=255)                                                                                               gcagTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGTTCTGCAgg                                                    <  3:67684/149‑1 (MQ=255)                                                                                                     gctgctCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGTTCTGCAGGCGGTcc                                              >  1:34052/1‑149 (MQ=255)                                                                                                     gctgctCCAGTACCTTGACGTCATCGACCAGCGCAAAACGGTGCTTGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTAGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGTTCTGCAGGCGGTcc                                              >  1:47374/1‑149 (MQ=255)                                                                                                        gctCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGTTCTGCAGGCGGTCCTGc                                           >  2:3960/1‑149 (MQ=255)                                                                                                         ctCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGTTCTGCAGGCGGTCCt                                             <  6:195766/146‑1 (MQ=255)                                                                                                         ctCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGTTCTGCAGGCGGTCCt                                             >  5:195766/1‑146 (MQ=255)                                                                                                                 ccTTGACGGCATCGACCAGCGCAAAGCGGTGCTTGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGTTCTGCAGGCGGTCCTGCTCGCCACgg                                  <  7:116739/149‑1 (MQ=255)                                                                                                                  cTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGTTCTGCAGGCGGTCCTGCTCGCCACgg                                  >  2:148075/1‑148 (MQ=255)                                                                                                                     gACGTCATCGACCAGCGCAAAGCGGTGCTTGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGTTCTGCAGGCGGTCCTGCTCGCCACGGTAg                               <  1:148075/148‑1 (MQ=255)                                                                                                                           aTCGACCAGCGCAAAGCGGTGCTTGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGTTCTGCAGGCGGTCCTGCTCGCCACGGTAGTGCTCg                         >  4:230765/1‑148 (MQ=255)                                                                                                                            tCGACCAGCGCAAAGCGGTGCTCGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGTTCTGCAGGCGGTCCTGCTCGCCACGGTAGTGCTCGGa                       <  4:211085/149‑1 (MQ=255)                                                                                                                                ccAGCGCAAAGCGGTGCTTGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGTTCTGCAGGCGGTCCTGCTCGCCACGGTAGTGCTCGGACTgg                   <  1:3960/149‑1 (MQ=255)                                                                                                                                      cAAAGCGGTGCTTGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGTTCTGCAGGCGGTc                                               >  8:81489/1‑115 (MQ=255)                                                                                                                                      cAAAGCGGTGCTTGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGTTCTGCAGGCGGTc                                               <  7:81489/115‑1 (MQ=255)                                                                                                                                       aaaGCGGTGCTTGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGTTCTGCAGGCGGTCCTGCTCGCCACGGTAGTGCTCGGACTGGATGAcc             <  1:251865/148‑1 (MQ=255)                                                                                                                                       aaaGCGGTGCTTGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGTTCTGCAGGCGGTCCTGCTCGCCACGGTAGTGCTCGGACTGGATGAcc             >  2:251865/1‑148 (MQ=255)                                                                                                                                       aaaGCGGTGCTTGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGTTCTGCAGGCGGTCCTGCTCGCCACGGTAGTGCTCGGACTGGATGAc              >  1:64876/1‑147 (MQ=255)                                                                                                                                        tgtcGGTGCTCAAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGTTCTGCAGGCGGTCCTGCTCGCCACGGTAGTGCTCGGACTGGATGACCAg           >  2:234879/4‑149 (MQ=255)                                                                                                                                         aGCGGTGCTTGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGTTCTGCAGGCGGTCCTGCTCGCCACGGTAGTGCTCGGACTGGATGACCAGc          <  2:134998/149‑1 (MQ=255)                                                                                                                                          tcGGTGCTCAAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGTTCTGCAGGCGGTCCTGCTCGCCACGGTAGTGCTCGGACTGGATGACCAGct         >  5:102781/2‑149 (MQ=255)                                                                                                                                             gTGCTTGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGTTCTGCAGGCGGTCCTGCTCGCCACGGTAGTGCTCGGACTGGATGACCAGCTCtt      >  5:171612/1‑149 (MQ=255)                                                                                                                                             gTGCTCAAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGTTCTGCAGGCGGTCCTGCTCGCCACGGTAGTGCTCGGACTGGATGACCAGCTCtt      <  6:102781/149‑1 (MQ=255)                                                                                                                                              tGCTTGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGTTCTGCAGGCGGTCCTGCTCGCCACGGTAGTGCTCGGACTGGATGACCAGCTCTTc     <  8:25608/149‑1 (MQ=255)                                                                                                                                                cTTGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCGTGGGTGCCGCCACGCTGGTCGTTCTGCAGGCGGTCCTGCTCGCCACGGTAGTGCTCGGACTGGATGACCAGCTCTTCCt   >  5:194173/1‑149 (MQ=255)                                                                                                                                                cTTGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGTTCTGCAGGCGGTCCTGCTCGCCACGGTAGTGCTCGGACTGGATGACCAGCTCTTcc    <  2:34052/148‑1 (MQ=255)                                                                                                                                                cTTGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGTTCTGCAGGCGGTCCTGCTCGCCACGGTAGTGCTCGGACTGGATGACCAGCTCTTcc    <  3:230765/148‑1 (MQ=255)                                                                                                                                                 ttGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGTTCTGCAGGCGGTCCTGCTCGCCACGGTAGTGCTCGGACTGGATGACCAGCTCTTCCtt  >  7:247248/1‑149 (MQ=255)                                                                                                                                                   gAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGTTCTGCAGGCGGTCCTGCTCGCCACGGTAGTGCTCGGACTGGATGAc              >  7:58653/1‑135 (MQ=255)                                                                                                                                                   gAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGTTCTGCAGGCGGTCCTGCTCGCCACGGTAGTGCTCGGACTGGATGAc              <  8:58653/135‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                           AGGTCGCGTACGCGCTCGTCGAGGTCTTCGGCACTGAACTGCCGTGACTGGGCCAGCGCGCCGGCCAGCTCGTCGTGGGCGTCTTTGGCCGCCAGCAGTTGCTGCTCCAGTACCTTGACGTCATCGACCAGCGCAAAGCGGTGCTTGAG‑‑‑CAC‑‑‑‑CGA‑CAGCTCGCCTAGCCAGCGCTGGATACCGGTAATTTCACGCTCCAGGCGCATCAGCTCCTGGGTGCCGCCACGCTGGTCGTTCTGCAGGCGGTCCTGCTCGCCACGGTAGTGCTCGGACTGGATGACCAGCTCTTCCTT  >  NC_002947/806182‑806474

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