Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	821,732	T→C	6.5%	V368A (GTG→GCG) ‡	PP_0707 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	821,732	0	T	C	6.5%	176.5 / 2.3	61	V368A (GTG→GCG) ‡	PP_0707	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (29/28);  new base C (2/2);  total (31/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.67e-01

AGCGTCCGTCCTGGCGGTTGCCGGCGATCCACGACGAAGTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGT‑GGCGCTCACCTTTGCAGGCGCGCTGCTGTGTTACCGCCGCACCCTGCGCAAGCATGA‑CCTGGAGCGCCCGACGGGCCTGGCCGGGCTGATCCCGTTCGTGATGGTCATCTGGC  >  NC_002947/821585‑821875                                                                                                                                                    |                                                                                                                                                  aGCGTCCGTCCTGGCGGTTGCCGGCGATCCACGACGAAGTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTg                                                                                                                                                  <  2:126857/149‑1 (MQ=255)       cgtacTGGCGGTTGCCGGCGATAAACGACGAAGTCGCAAGTGAGCTGGGCTGGTTCCCGAAGGTGATGGCATTGGAACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGACTCGCGCTGACGGTGGACGt                                                                                                                                             >  2:102610/1‑148 (MQ=255)         tccTGGCGGTTGCCGGCGATCCACGACGAAGTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCaa                                                                                                                                          >  6:167026/1‑149 (MQ=255)                  ttGCCGGCGATCCACGACGAAGTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGAc                                                                                                                                 >  6:114552/1‑149 (MQ=255)                      cGGCGATCCACGACGAAGTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCaa                                                                                                                                          <  1:158900/136‑1 (MQ=255)                      cGGCGATCCACGACGAAGTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCaa                                                                                                                                          >  2:158900/1‑136 (MQ=255)                           aTCCACGACGAAGTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCtggt                                                                                                                        <  3:25557/149‑1 (MQ=255)                              cACGACGAAGTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCtggtggt                                                                                                                     >  4:170301/1‑149 (MQ=255)                               acgacgAAGTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCtggtggt‑g                                                                                                                   >  2:42498/1‑149 (MQ=255)                                 gacgaAGTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGT‑Ggc                                                                                                                 >  8:213685/1‑149 (MQ=255)                                    gaAGTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGTTGACGGCCACCGTCAGCg                                                                                                                                        >  6:221649/1‑124 (MQ=255)                                    gaAGTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGTTGACGGCCACCGTCAGCg                                                                                                                                        >  6:218020/1‑124 (MQ=255)                                    gaAGTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGTTGACGGCCACCGTCAGCg                                                                                                                                        <  5:221649/124‑1 (MQ=255)                                    gaAGTCGACAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGTTGACGGCCACCGTCAGCg                                                                                                                                        <  5:218020/124‑1 (MQ=255)                                       gTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGT‑GGCGCTCAc                                                                                                           >  1:250338/1‑149 (MQ=255)                                             aGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGT‑GGCGCTCACCTTTGc                                                                                                     <  1:102610/149‑1 (MQ=255)                                                       gctggTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGT‑GGCGCTCACCTTTGCAGGCGCgctg                                                                                           <  3:170301/149‑1 (MQ=255)                                                           gTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGTCGGCGCTCACCTTTGCAGGCGCgctgct                                                                                         <  1:148775/148‑1 (MQ=255)                                                           gTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGT‑GGCGCTCACCTTTGCAGGCGCGCTGCtgt                                                                                       >  6:199641/1‑149 (MQ=255)                                                           gTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGT‑GGCGCTCACCTTTGCAGGCGCGCTGCtgt                                                                                       <  5:167026/149‑1 (MQ=255)                                                            ttCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGT‑GGCGCTCACCTTTGCAGGCGCGCTGCtgtg                                                                                      >  4:252280/1‑149 (MQ=255)                                                                    ggTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGT‑GGCGCTCACCTTTGCAGGCGCGCTGCtgt                                                                                       <  3:137818/140‑1 (MQ=255)                                                                    ggTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGT‑GGCGCTCACCTTTGCAGGCGCGCTGCtgt                                                                                       >  4:137818/1‑140 (MQ=255)                                                                       gCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGT‑GGCGCTCACCTTTGCAGGCGCGCTGCTGTGTTACCGCCGCa                                                                           >  1:107784/1‑149 (MQ=255)                                                                         tGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGT‑GGCGCTCACCTTTGCAGGCGCGCTGCTGTGTTACCGCCGCAcc                                                                         >  5:48635/1‑149 (MQ=255)                                                                           gCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGT‑GGCGCTCACCTTTGCAGGCGCGCTGCTGTGTTACCGCCGCAcc                                                                         <  1:42498/147‑1 (MQ=255)                                                                           gCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGT‑GGCGCTCACCTTTGCAGGCGCGCTGCTGTGTTACCGCCGCAcc                                                                         <  2:35330/147‑1 (MQ=255)                                                                           gCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGT‑GGCGCTCACCTTTGCAGGCGCGCTGCTGTGTTACCGCCGCACCCt                                                                       >  6:219481/1‑149 (MQ=255)                                                                                 gCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGT‑GGCGCTCACCTTTGCAGGCGAGCTGCTGTGTTACAGCAGCACCATGCGCaa                                                                 >  8:28002/1‑149 (MQ=255)                                                                                 gCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGAGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGT‑GGCGCTCACCTTTGCAGGCGCGCTGCTGTGTTACAGACGCACCCTgcgc                                                                   >  8:32191/1‑147 (MQ=255)                                                                                  gctccaTGGTCATCATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGT‑GGCGCTCACCTTTGCAGGCGCGCTGCTGTGTTAc                                                                                  <  7:76627/126‑1 (MQ=255)                                                                                  gctccaTGGTCATCATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGT‑GGCGCTCACCTTTGCAGGCGCGCTGCTGTGTTAc                                                                                  >  8:76627/6‑131 (MQ=255)                                                                                  cACTTATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGTTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGT‑GGCGCTCACCTTTGCAGGCGCGCTGCTGTGTTACCGCCGCACCCTGCGCAAg                                                                <  7:28002/149‑1 (MQ=255)                                                                                  cACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCCGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGT‑GGCGCTCACCTTTGCAGGCCCGCTGCTGTGTTACCGCCGCACCCTGCGCAAg                                                                <  7:32191/149‑1 (MQ=255)                                                                                   aCTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGT‑GGCGCTCACCTTTGCAGGCGCGCTGCTGTGTTACCGCCGCACCCTGCGCAAg                                                                >  5:105282/1‑148 (MQ=255)                                                                                   aCTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGT‑GGCGCTCACCTTTGCAGGCGCGCTGCTGTGTTACCGCCGCACCCTGCGCAAg                                                                >  8:232338/1‑148 (MQ=255)                                                                                     tGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGT‑GGCGCTCACCTTTGCAGGCGCGCTGCTGTGTTACCGCCGCACCCTGCGCAAGCAt                                                             >  2:85566/1‑149 (MQ=255)                                                                                        tGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGT‑GGCGCTCACCTTTGCAGGCGCGCTGCTGTGTTACCGCCGCACCCTGCGCAAGCATGA‑c                                                         <  5:199641/149‑1 (MQ=255)                                                                                           tCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGT‑GGCGCTCACCTTTGCAGGCGCGCTGCTGTGTTACCGCCGCACCCTGCGCAAGCATGA‑CCTg                                                      <  3:81785/149‑1 (MQ=255)                                                                                             atgatgACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGT‑GGCGCTCACCTTTGCAGGCGCGCTGCTGTGTTACCGCCGCACCCTGCGCAAGCATGA‑CCTGGa                                                    <  7:213685/149‑1 (MQ=255)                                                                                             atgatgACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGT‑GGCGCTCACCTTTGCAGGCGCGCTGCTGTGTTACCGCCGCACCCTGCGCAAGCATGA‑CCTGGa                                                    <  5:219481/149‑1 (MQ=255)                                                                                              tgatgaCCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGT‑GGCGCTCACCTTTGCAGGCGCGCTGCTGTGTTACCGCCGCACCCTGCGCAAGCATGA‑CCTGGAg                                                   >  5:35733/1‑149 (MQ=255)                                                                                                   aCCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGT‑GGCGCTCACCTTTGCAGGCGCGCTGCTGTGTTACCGCCGCACCCTGCGCAAGCATGA‑CCTGGAGCGccc                                              <  1:85566/149‑1 (MQ=255)                                                                                                    ccATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGT‑GGCGCTCACCTTTGCAGGCGCGCTGCTGTGTTACCGCCGCa                                                                           >  1:222418/1‑120 (MQ=255)                                                                                                    ccATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGT‑GGCGCTCACCTTTGCAGGCGCGCTGCTGTGTTACCGCCGCa                                                                           <  2:222418/120‑1 (MQ=255)                                                                                                     cATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGT‑GGCGCTCACCTTTGc                                                                                                     <  6:13672/93‑1 (MQ=255)                                                                                                     cATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGT‑GGCGCTCACCTTTGc                                                                                                     >  5:13672/1‑93 (MQ=255)                                                                                                          aGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGT‑GGCGCTCACCTTTGCAGGCGCGCTGCTGTGTTACCGCCGCACCCTGCGCAAGCATGA‑CCTGGAGCGCCCGACGGGc                                       <  6:105282/149‑1 (MQ=255)                                                                                                            cgcATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGT‑GGCGCTCACCTTTGCAGGCGCGCTGCTGTGTTACCGCCGCACCCTGCGCAAGCATGA‑CCTGGAGCGCCCGACGGGCCt                                     <  3:252280/149‑1 (MQ=255)                                                                                                                tGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGT‑GGCGCTCACCTTTGCAGGCGCGCTGCTGTGTTACCGCCGCACCCTGCGCAAGCATGA‑CCTGGAGCGCCCGACGGGCCTGGcc                                 >  8:99665/1‑149 (MQ=255)                                                                                                                 gAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGCTGGT‑GGCGCTCACCTTTGCAGGCGCGCTGCTGTGTTACCGCCGCACCCTGCGCAAGCATGA‑CCTGGAGAGCCCGACGGGCATGTCCg                                >  5:26836/1‑149 (MQ=255)                                                                                                                               gTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGT‑GGCGCTCACCTTTGCAGGCGCGCTGCTGTGTTACCGCCGCACCCTGCGCAAGCATGA‑CCTGGAGCGCCCGACGGGCCTGGCCGGGCTGATCCCGTTc                  <  7:99665/149‑1 (MQ=255)                                                                                                                               gTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGT‑GGCGCTCACCTTTGCAGGCGCGCTGCTGTGTTACCGCCGCACCCTGCGCAAGCATGA‑CCTGGAGCGCCCGACGGGCCTGGCCGGGCTGATCCCGTTc                  <  5:114552/149‑1 (MQ=255)                                                                                                                                tGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGT‑GGCGCTCACCTTTGCAGGCGCGCTGCTGTGTTACCGCCGCACCCTGCGCAAGCATGA‑CCTGGAGCGCCCGACGGGCCTGGCCGGGCTGATCCCGTTCg                 >  2:74095/1‑149 (MQ=255)                                                                                                                                  ccAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGT‑GGCGCTCACCTTTGCAGGCGCGCTGCTGTGTTACCGCCGCCCCCTGCGCAAGCATGC‑CCTGGAGCGCCGGACGGGCCTGGCCGGGCTGATCCCGTTCGt                >  1:228825/1‑148 (MQ=255)                                                                                                                                  ccAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGT‑GGCGCTCACCTTTGCAGGCGCGCTGCTGTGTTACCGCCGCACCCTGCGCAAGCATGA‑CCTGGAGCGCCCGACGGGCCTGGCCGGGCTGATCCCGTTCGt                >  8:34232/1‑148 (MQ=255)                                                                                                                                  ccAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGT‑GGCGCTCACCTTTGCAGGCGCGCTGCTGTGTTACCGCCGCACCCTGCGCAAGCATGA‑CCTGGAGCGCCCGACGGGCCTGGCCGGGCTGATCCCGTTCGTg               <  6:48635/149‑1 (MQ=255)                                                                                                                                     gCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGT‑GGCGCTCACCTTTGCAGGCGCGCTGCTGTGTTACCGCCGCACCCTGCGCAAGCATGA‑CCTGGAGCGCCCGACGGGCCTGGCCGGGCTGATCCCGTTCGTGATg            <  7:232338/149‑1 (MQ=255)                                                                                                                                     gCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGT‑GGCGCTCACCTTTGCAGGCGCGCTGCTGTGTTACCGCCGCACCCTGCGCAAGCATGA‑CCTGGAGCGCCCGACGGGCCTGGCCGGGCTGATCCCGTTCGTGATg            <  6:35733/149‑1 (MQ=255)                                                                                                                                        tGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGT‑GGCGCTCACCTTTGCAGGCGCGCTGCTGTGTTACCGCCGCACCCTGCGCAAGCATGA‑CCTGGAGCGCCCGACGGGCCTGGCCGGGCTGATCCCGTTCGTGATGGt          >  3:97931/1‑148 (MQ=255)                                                                                                                                                aCGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGT‑GGCGCTCACCTTTGCAGGCGCGCTGCTGTGTTACCGCCGCACCCTGCGCAAGCATGACCCTGGAGCGCCCGACGGGCCTGGCCGGGCTGATCCCGTTCGTGATGGTCATCTGGc  <  1:74095/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                  AGCGTCCGTCCTGGCGGTTGCCGGCGATCCACGACGAAGTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGT‑GGCGCTCACCTTTGCAGGCGCGCTGCTGTGTTACCGCCGCACCCTGCGCAAGCATGA‑CCTGGAGCGCCCGACGGGCCTGGCCGGGCTGATCCCGTTCGTGATGGTCATCTGGC  >  NC_002947/821585‑821875

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