Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,639,337	T→C	6.2%	V117A (GTG→GCG)	PP_3206 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,639,337	0	T	C	6.2%	188.7 / 2.1	65	V117A (GTG→GCG)	PP_3206	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (30/30);  new base C (2/2);  total (33/32)
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.00e+00

ACGGTGTCGATGTGGTCTTCCATCTCGTCAGCGTGCCTGGCGGCGCTGCCGAGACGCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTA‑CGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAACGCATGGTGGAGATCGCGTTGCAGGATCTTGCCCGCCGTGGCGAGGTT  >  NC_002947/3639193‑3639482                                                                                                                                                 |                                                                                                                                                   aCGGTGTCGATGTGGTCTTCCATCTCGTCAGCGTGCCTGGCGGCGCTGCCGAGACGCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGAGCTGgtgt                                                                                                                                                  <  8:188975/147‑1 (MQ=255)   ggTGTCGATGTGGTCTTCCATCTCGTCAGCGTGCCTGGCGGCGCTGCCGAGACGCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTA‑CGc                                                                                                                                             >  8:372922/1‑149 (MQ=255)           tgtgGTCTTCCATCTCGTCAGCGTGCCTGGCGGCGCTGCCGAGACGCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTA‑CGCCAGCAGCg                                                                                                                                     >  5:383119/1‑149 (MQ=255)               gTCTTCCATCTCGTCAGCGTGCCTGGCGGCGCTGCCGAGACGCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGCGCTGGTGTA‑CGCCAGCACCggtg                                                                                                                                  <  2:14241/148‑4 (MQ=255)                  ttCCATCTCGTCAGCGTGCCTGGCGGCGCTGCCGAGACGCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTA‑CGCCAGCAGCGTCGCGgt                                                                                                                              >  5:84637/1‑149 (MQ=255)                   tCCATCTCGTCAGCGTGCCTGGCGGCGCTGCCGAGACGCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTA‑CGCCAGCAGCGTCGCGgtg                                                                                                                             >  8:208630/1‑149 (MQ=255)                    ccATCTCGTCAGCGTGCCTGGCGGCGCTGCCGAGACGCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTA‑CGCCAGCAGCGTCGCGgtg                                                                                                                             >  4:293773/1‑148 (MQ=255)                    ccATCTCGTCAGCGTGCCTGGCGGCGCTGCCGAGACGCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTA‑CGCCAGCAGCGTCGCGgt                                                                                                                              <  4:327964/147‑1 (MQ=255)                       tctcGTCAGCGTGCCTGGCGGCGCTGCCGAGACGCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTA‑CGCCAGCAGCGTCGCGGTGTAcg                                                                                                                         >  4:81377/1‑149 (MQ=255)                       tctcGTCAGCGTGCCTGGCGGCGCTGCCGAGACGCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTA‑CGCCAGCAGCGTCGCGGTGTAcg                                                                                                                         >  4:86976/1‑149 (MQ=255)                            tCAGCGTGCCTGGCGGCGCTGCCGAGACGCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTA‑CGCCAGCAGCGTCGCGGTGTACggcggc                                                                                                                    <  8:262564/149‑1 (MQ=255)                              aGCGTGCCTGGCGGCGCTGCCGAGACGCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTA‑CGCCAGCAGCGTCGCGGTGTACGGCGGCaa                                                                                                                  <  6:383119/149‑1 (MQ=255)                              aGCGTGCCTGGCGGCGCTGCCGAGACGCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTA‑CGCCAACAGCGTCGCGGTGTACGGCGGCaa                                                                                                                  <  3:156823/149‑1 (MQ=255)                                 gTGCCTGGCGGCGCTGCCGAGACGCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTA‑CGCCAGCAGCGTCGCGGTGTACGGCGGCAAGtt                                                                                                               <  7:208630/149‑1 (MQ=255)                                    ccTGGCGGCGCTGCCGAGACGCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTA‑CGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGc                                                                                                             <  7:372922/148‑1 (MQ=255)                                       ggcggcGCTGCCGAGACGCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACAGGTGCTGGAGGAGCCCCAGCAGCGTAGCGGTGTACGGCGGCAAGTTGCCGa                                                                                                          >  1:342076/1‑149 (MQ=255)                                           gcgcTGCCGAGACGCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTA‑CGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGAcgcg                                                                                                      >  6:338960/1‑148 (MQ=255)                                              cTGCCGAGACGCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTA‑CGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATg                                                                                                  >  5:189215/1‑149 (MQ=255)                                              cTGCCGAGACGCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTA‑CGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATg                                                                                                  <  3:398806/149‑1 (MQ=255)                                              cTGCCGAGACGCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGCGTA‑CACCAGCACCggt                                                                                                                                   >  6:273028/1‑114 (MQ=255)                                              cTGCCGAGACGCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGCGTA‑CACCAGCACCggt                                                                                                                                   <  5:273028/116‑3 (MQ=255)                                               tGCCGAGACGCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTA‑CGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATg                                                                                                  >  5:221507/1‑148 (MQ=255)                                                  cGAGACGCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTA‑CGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGAc                                                                                               >  2:291299/1‑148 (MQ=255)                                                      atgCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCCAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTA‑CGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATTGACGAAgg                                                                                          <  6:269721/147‑1 (MQ=255)                                                         cAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTA‑CGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGgc                                                                                         >  5:321212/1‑147 (MQ=255)                                                              cGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTA‑CGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGgccagcca                                                                                   <  6:221507/148‑1 (MQ=255)                                                               gAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTA‑CGCCAGCAGCGTCGCGGTGTACGGCGGCaa                                                                                                                  <  3:364143/116‑1 (MQ=255)                                                               gAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTA‑CGCCAGCAGCGTCGCGGTGTACGGCGGCaa                                                                                                                  >  4:364143/1‑116 (MQ=255)                                                                aaCGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGGGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTA‑CGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCg                                                                                                           >  5:329286/1‑122 (MQ=255)                                                                aaCGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGCGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTA‑CGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCg                                                                                                           <  6:329286/122‑1 (MQ=255)                                                                    gggTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTA‑CGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTcgc                                                                            <  6:189215/149‑1 (MQ=255)                                                                          ccAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTA‑CGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGc                                                                      >  6:354061/1‑149 (MQ=255)                                                                            aGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTA‑CGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAg                                                                                           >  3:304148/1‑126 (MQ=255)                                                                            aGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTA‑CGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAg                                                                                           <  4:304148/126‑1 (MQ=255)                                                                             ggTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTA‑CGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGt                                                                   >  4:89845/1‑149 (MQ=255)                                                                              gTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTA‑CGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCGGCTGTc                                                                  <  6:84637/149‑1 (MQ=255)                                                                                 aaCCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTA‑CGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCt                                                                <  3:89845/148‑1 (MQ=255)                                                                                    cTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTA‑CGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATgc                                                            >  3:376258/1‑149 (MQ=255)                                                                                        aagcaagcCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTA‑CGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATgcggcg                                                        >  2:109677/1‑149 (MQ=255)                                                                                         agcaagcCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTA‑c                                                                                                                                               >  8:385035/1‑61 (MQ=255)                                                                                         agcaagcCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTA‑c                                                                                                                                               <  7:385035/61‑1 (MQ=255)                                                                                          gcaagcCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTA‑CGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGca                                                      >  1:231146/1‑149 (MQ=255)                                                                                              gcCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTA‑CGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACaaa                                                  <  3:86976/149‑1 (MQ=255)                                                                                              gcCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTA‑CGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACaaa                                                  <  3:81377/149‑1 (MQ=255)                                                                                               ccTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTA‑CGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAAc                                                 >  1:361253/1‑149 (MQ=255)                                                                                                cTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTA‑CGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAACg                                                <  4:376258/149‑1 (MQ=255)                                                                                                     gctgctCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTA‑CGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAACGCAtgg                                           <  5:338960/149‑1 (MQ=255)                                                                                                        gctCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTA‑CGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAACGCAtggtg                                         >  3:262582/1‑148 (MQ=255)                                                                                                         ctCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTA‑CGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAACGCAtggtg                                         >  5:49272/1‑147 (MQ=255)                                                                                                            gACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTA‑CGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAACGCAt                                             <  5:322065/140‑1 (MQ=255)                                                                                                            gACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTA‑CGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAACGCAt                                             >  6:322065/1‑140 (MQ=255)                                                                                                            gACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTA‑CGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAACGCATGGTGgaga                                     <  1:109677/148‑1 (MQ=255)                                                                                                               cAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTA‑CGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAACGCATGGTGGAGATcgc                                 >  5:77336/1‑149 (MQ=255)                                                                                                               cAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTA‑CGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAACGCATGGTGGAGATcgc                                 <  1:347217/149‑1 (MQ=255)                                                                                                                aaCTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTA‑CGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGca                                                      <  5:384296/127‑1 (MQ=255)                                                                                                                aaCTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTA‑CGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGca                                                      >  6:384296/1‑127 (MQ=255)                                                                                                                     cGTAATCCCCTCTGTCCACCGGTGCTGGTGTA‑CGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAACGCATGGTGGAGATCGCGTTg                             >  4:273508/1‑147 (MQ=255)                                                                                                                     cGTAATCCCCTCTGTCCACCGGTGCTGGTGTA‑CGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAACGCATGGTGGAGATCGCGTTGc                            <  4:262582/148‑1 (MQ=255)                                                                                                                     cGTAATCCCCTCTGTCCACCGGTGCTGGTGTA‑CGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAACGCATGGTGGAGATCGCGTTGCa                           <  6:321212/149‑1 (MQ=255)                                                                                                                       tAATCCCCTCTGTCCACCGGTGCTGGTGTA‑CGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAACGCATGGTGGAGATCGCGTTGCAgg                         <  6:77336/149‑1 (MQ=255)                                                                                                                            cccTCTGTCCACCGGTGCTGGTGTA‑CGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAACGCATGGTGGAGATCGCGTTGCAGGATCtt                    >  7:358918/1‑149 (MQ=255)                                                                                                                              ctctGTCCACCGGTGCTGGTGTA‑CGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAACGCATGGTGGAGATCGCGTTGCAGGATCTTGc                  <  3:293773/149‑1 (MQ=255)                                                                                                                                ctGTCCACCGGTGCTGGTGTA‑CGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAACGCATGGTGGAGATCGCGTTGCAGGATCTTGccc                >  3:92787/1‑149 (MQ=255)                                                                                                                                      gacgctGCTGGCGTA‑CACCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAACGCATGGTGGAGATCGCGTTGCAGGATCTTGCCCGCCGTg          >  1:315404/6‑149 (MQ=255)                                                                                                                                           tGCTGGTGTA‑CGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAACGCATGGTGGAGATCGCGTTGCAGGATCTTGCCCGCCGTGGCGAg     <  2:231146/149‑1 (MQ=255)                                                                                                                                             cTGGCGTA‑CACCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAACGCATGGTGGAGATCGCGTTGCAGGATCTTGCCCGCCGTGGCGAg     <  2:315404/147‑1 (MQ=255)                                                                                                                                                gtgtA‑CGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAACGCATGGTGGAGATCGCGTTGCAGGATCTTGCCCGCCGTGGCGAGGtt  <  4:92787/147‑1 (MQ=255)                                                                                                                                                 tgtA‑CGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAACGCATGGTGgaga                                     >  2:188860/1‑111 (MQ=255)                                                                                                                                                 tgtA‑CGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAACGCATGGTGgaga                                     <  1:188860/111‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                   ACGGTGTCGATGTGGTCTTCCATCTCGTCAGCGTGCCTGGCGGCGCTGCCGAGACGCAGTACGAACGGGGTTACCAGGTAAACCTGCAAGCAAGCCTGGAGCTGCTCGACCAACTGCGTAATCCCCTCTGTCCACCGGTGCTGGTGTA‑CGCCAGCAGCGTCGCGGTGTACGGCGGCAAGTTGCCGACGCGCATGGACGAAGGCCAGCCAGCGTCGCCGCAGCTGTCCTATGCGGCGCACAAACGCATGGTGGAGATCGCGTTGCAGGATCTTGCCCGCCGTGGCGAGGTT  >  NC_002947/3639193‑3639482

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