Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	821,741	A→G	6.8%	N371S (AAC→AGC)	PP_0707 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	821,741	0	A	G	6.8%	173.6 / 2.3	59	N371S (AAC→AGC)	PP_0707	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (27/28);  new base G (2/2);  total (29/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 = 9.90e-01

TCCTGGCGGTTGCCGGCGATCCACGACGAAGTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGT‑GGCGCTCACCTTTGCAGGCGCGCTGCTGTGTTACCGCCGCACCCTGCGCAAGCATGA‑CCTGGAGCGCCCGACGGGCCTGGCCGGGCTGATCCCGTTCGTGATGGTCATCTGGCTGACGTTG  >  NC_002947/821593‑821883                                                                                                                                                     |                                                                                                                                                 tccTGGCGGTTGCCGGCGATCCACGACGAAGTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCaa                                                                                                                                                  >  6:167026/1‑149 (MQ=255)          ttGCCGGCGATCCACGACGAAGTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGAc                                                                                                                                         >  6:114552/1‑149 (MQ=255)              cGGCGATCCACGACGAAGTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCaa                                                                                                                                                  >  2:158900/1‑136 (MQ=255)              cGGCGATCCACGACGAAGTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCaa                                                                                                                                                  <  1:158900/136‑1 (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:218020/1‑124 (MQ=255)                            gaAGTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGTTGACGGCCACCGTCAGCg                                                                                                                                                >  6:221649/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                                                                                               <  5:167026/149‑1 (MQ=255)                                                   gTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGT‑GGCGCTCACCTTTGCAGGCGCGCTGCtgt                                                                                               >  6:199641/1‑149 (MQ=255)                                                    ttCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGT‑GGCGCTCACCTTTGCAGGCGCGCTGCtgtg                                                                                              >  4:252280/1‑149 (MQ=255)                                                            ggTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGT‑GGCGCTCACCTTTGCAGGCGCGCTGCtgt                                                                                               >  4:137818/1‑140 (MQ=255)                                                            ggTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGT‑GGCGCTCACCTTTGCAGGCGCGCTGCtgt                                                                                               <  3:137818/140‑1 (MQ=255)                                                               gCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGT‑GGCGCTCACCTTTGCAGGCGCGCTGCTGTGTTACCGCCGCa                                                                                   >  1:107784/1‑149 (MQ=255)                                                                 tGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGT‑GGCGCTCACCTTTGCAGGCGCGCTGCTGTGTTACCGCCGCAcc                                                                                 >  5:48635/1‑149 (MQ=255)                                                                   gCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGT‑GGCGCTCACCTTTGCAGGCGCGCTGCTGTGTTACCGCCGCAcc                                                                                 <  2:35330/147‑1 (MQ=255)                                                                   gCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGT‑GGCGCTCACCTTTGCAGGCGCGCTGCTGTGTTACCGCCGCAcc                                                                                 <  1:42498/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                                                                                          >  8:76627/6‑131 (MQ=255)                                                                          gctccaTGGTCATCATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGT‑GGCGCTCACCTTTGCAGGCGCGCTGCTGTGTTAc                                                                                          <  7:76627/126‑1 (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                                                            <  5:219481/149‑1 (MQ=255)                                                                                     atgatgACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGT‑GGCGCTCACCTTTGCAGGCGCGCTGCTGTGTTACCGCCGCACCCTGCGCAAGCATGA‑CCTGGa                                                            <  7:213685/149‑1 (MQ=255)                                                                                      tgatgaCCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGT‑GGCGCTCACCTTTGCAGGCGCGCTGCTGTGTTACCGCCGCACCCTGCGCAAGCATGA‑CCTGGAg                                                           >  5:35733/1‑149 (MQ=255)                                                                                           aCCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGT‑GGCGCTCACCTTTGCAGGCGCGCTGCTGTGTTACCGCCGCACCCTGCGCAAGCATGA‑CCTGGAGCGccc                                                      <  1:85566/149‑1 (MQ=255)                                                                                            ccATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGT‑GGCGCTCACCTTTGCAGGCGCGCTGCTGTGTTACCGCCGCa                                                                                   <  2:222418/120‑1 (MQ=255)                                                                                            ccATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGT‑GGCGCTCACCTTTGCAGGCGCGCTGCTGTGTTACCGCCGCa                                                                                   >  1:222418/1‑120 (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                    <  6:35733/149‑1 (MQ=255)                                                                                                                             gCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGT‑GGCGCTCACCTTTGCAGGCGCGCTGCTGTGTTACCGCCGCACCCTGCGCAAGCATGA‑CCTGGAGCGCCCGACGGGCCTGGCCGGGCTGATCCCGTTCGTGATg                    <  7:232338/149‑1 (MQ=255)                                                                                                                                tGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGT‑GGCGCTCACCTTTGCAGGCGCGCTGCTGTGTTACCGCCGCACCCTGCGCAAGCATGA‑CCTGGAGCGCCCGACGGGCCTGGCCGGGCTGATCCCGTTCGTGATGGt                  >  3:97931/1‑148 (MQ=255)                                                                                                                                        aCGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGT‑GGCGCTCACCTTTGCAGGCGCGCTGCTGTGTTACCGCCGCACCCTGCGCAAGCATGACCCTGGAGCGCCCGACGGGCCTGGCCGGGCTGATCCCGTTCGTGATGGTCATCTGGc          <  1:74095/149‑1 (MQ=255)                                                                                                                                               ccGTCAACGGGCTGACCGCGCTGGTGGT‑GGCGCTCACCTTTGCAGGCGCGCTGCTGTGTTACCGCCGCACCCTGCGCAAGCATGA‑CCTGGAGCGCCCGACGGGCCTGGCCGGGCTGATCCCGTTCGTGATGGTCATCTGGCTGACGTTg  <  2:250338/149‑1 (MQ=255)                                                                                                                                                   cAACGGGCTGACCGCGCTGGTGGT‑GGCGCTCACCTTTGCAGGCGCGCTGCTGTGTTACCGCCGCACCCTGCGCAAGCATGA‑CCTGGAGCGCCCGACGGGCCTGGCCGGGCTGATCCCGTTc                          >  1:186795/1‑121 (MQ=255)                                                                                                                                                   cAACGGGCTGACCGCGCTGGTGGT‑GGCGCTCACCTTTGCAGGCGCGCTGCTGTGTTACCGCCGCACCCTGCGCAAGCATGA‑CCTGGAGCGCCCGACGGGCCTGGCCGGGCTGATCCCGTTc                          <  2:186795/121‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                 TCCTGGCGGTTGCCGGCGATCCACGACGAAGTCGCCAGTGCGCTGGGCTGGTTCCCGAAGGTGCTGGCATTGGCACTGATGGTCATGATGACCATGGAGCGCATGAACAGCGTGATCGGTGCCAGCCTGGCGCTGACGGTGGCCGTCAACGGGCTGACCGCGCTGGTGGT‑GGCGCTCACCTTTGCAGGCGCGCTGCTGTGTTACCGCCGCACCCTGCGCAAGCATGA‑CCTGGAGCGCCCGACGGGCCTGGCCGGGCTGATCCCGTTCGTGATGGTCATCTGGCTGACGTTG  >  NC_002947/821593‑821883

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