Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,938,839	C→T	6.2%	E73K (GAG→AAG)	PP_1738 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,938,839	0	C	T	6.2%	189.0 / 2.1	64	E73K (GAG→AAG)	PP_1738	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (33/27);  new base T (2/2);  total (35/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.87e-01

AACAGCTTTTTCAGGTATTCATTGATGCTGCTGCTGTAGGGCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCT‑GGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGTAAGAGGGTAGCCATAGCTGGCAGCGT  >  NC_002947/1938691‑1938979                                                                                                                                                     |                                                                                                                                              aaCAGCTTTTTCAGGTATTCATTGATGCTGCTGCTGTAGGGCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTc                                                                                                                                               >  6:213188/1‑149 (MQ=255)  aCAGCTTTTTCAGGTATTCATTGATGCTGCTGCTGTAGGGCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTcg                                                                                                                                              <  5:213188/149‑1 (MQ=255)             aGGTATTCATTGATGCTGCTGCTGTAGGGCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGgcg                                                                                                                                   <  7:57819/149‑1 (MQ=255)             aGGTATTCATTGATGCTGCTGCTGTAGGGCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGgc                                                                                                                                    <  4:128626/148‑1 (MQ=255)              ggTATTCATTGATGCTGCTGCTGTAGGGCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGgcg                                                                                                                                   >  3:207004/1‑148 (MQ=255)               gTATTCATTGATGCTGCTGCTGTAGGGCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGgcgc                                                                                                                                  >  2:89555/1‑148 (MQ=255)                 aTTCATTGATGCTGCTGCTGTAGGGCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTTGCGCCCGGAGCGCGAg                                                                                                                               >  1:130928/1‑149 (MQ=255)                 aTTCATTGATGCTGCTGCTGTAGGGCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCgcgc                                                                                                                                           <  3:309260/137‑1 (MQ=255)                 aTTCATTGATGCTGCTGCTGTAGGGCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCgcgc                                                                                                                                           >  4:309260/1‑137 (MQ=255)                     aTTGATGCTGCTGCTGTAGGGCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTTGCGCCCGGAGCGCGAGTTCa                                                                                                                           <  2:130928/149‑1 (MQ=255)                     aTTGATGCTGCTGCTGTAGGGCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGtt                                                                                                                             >  6:344334/1‑147 (MQ=255)                              gctgctgTAGGGCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGt                                                                                                                  >  5:92682/1‑149 (MQ=255)                               ctgctgTAGGGCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTc                                                                                                                 <  2:332907/149‑1 (MQ=255)                                  ctgTAGGGCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCgg                                                                                                               >  3:68026/1‑148 (MQ=255)                                     tAGGGCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACt                                                                                                            <  6:92682/148‑1 (MQ=255)                                       gggCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGa                                                                                                          <  1:89555/148‑1 (MQ=255)                                       gggCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAg                                                                                                         >  4:49318/1‑149 (MQ=255)                                       gggCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAg                                                                                                         >  7:119423/1‑149 (MQ=255)                                            cgcCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGAtt                                                                                                    <  3:49318/149‑1 (MQ=255)                                               cGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCa                                                                                                 >  7:224867/1‑149 (MQ=255)                                                 gTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCAtc                                                                                               <  8:119423/149‑1 (MQ=255)                                                    ccGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCAtcttc                                                                                            >  8:311480/1‑149 (MQ=255)                                                     cGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCAtcgtcg                                                                                           >  5:103042/1‑149 (MQ=255)                                                          atgatgAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGAGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCT‑ggg                                                                                     <  4:68026/149‑1 (MQ=255)                                                           tgatgaAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCT‑ggg                                                                                     >  3:115633/1‑148 (MQ=255)                                                             atgaAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCT‑GGGCAg                                                                                  <  5:344334/149‑1 (MQ=255)                                                               gaAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCT‑GGGCAGCt                                                                                >  6:138666/1‑149 (MQ=255)                                                                  gATCAGCGGCGCTTCGTGATACTGGAGGGCCACGCGGTACTTGAGCTTCTTCACTGCACAGAAGTTGTCGGGAAACGTGAACTCGCGCTCCGGGCTCAAGTTCAGGATCTAGTCGGACTGAGAGATTTCATCGTAGCT‑GCGCAGCTCAg                                                                             >  2:1701/1‑149 (MQ=255)                                                                    tCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCT‑GGGCAGCGCAggg                                                                           <  7:311480/149‑1 (MQ=255)                                                                     cAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCT‑GGGCAGCTCAggg                                                                           <  4:207004/148‑1 (MQ=255)                                                                     cAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGATCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCT‑GGGCAGCTCAGGgc                                                                          >  8:25115/1‑149 (MQ=255)                                                                       gcggcgCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCT‑GGGCAGCTCAGGgcgc                                                                        <  5:138666/149‑1 (MQ=255)                                                                            gcTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCT‑GGGCAGCTCAGGGCGCTTGTc                                                                   >  6:169315/1‑149 (MQ=255)                                                                              ttCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCT‑GGGCAGCTCAGGGCGCTTGTCcg                                                                 >  3:193600/1‑149 (MQ=255)                                                                                  tGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCT‑GGGCCGAGCAGGGCGCGTGTCgggcgg                                                             >  8:146788/1‑149 (MQ=255)                                                                                       cTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCT‑GGGCAGCTCAGGGCGCTTGTCCGGCGGCgtgg                                                        >  8:64955/1‑149 (MQ=255)                                                                                        tGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCT‑GGGCAGCTCAGGGCGCTTGTCCGGCGGCgtggt                                                       >  4:267608/1‑149 (MQ=255)                                                                                        tGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCT‑GGGCAGCTCAGGGCGCTTGTCCGGCGGCgtggt                                                       >  3:164190/1‑149 (MQ=255)                                                                                              gCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCT‑GGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGAt                                                 <  4:164190/149‑1 (MQ=255)                                                                                                    cGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCT‑GGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGc                                           >  4:250241/1‑149 (MQ=255)                                                                                                        aCTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCATGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCT‑GGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTc                                        >  4:262879/1‑148 (MQ=255)                                                                                                            gAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCT‑GGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGaa                                     <  7:64955/147‑1 (MQ=255)                                                                                                              gCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCT‑GGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACggg                                 <  2:336479/149‑1 (MQ=255)                                                                                                                ttcttcACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCT‑GGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGtt                               >  2:275096/1‑149 (MQ=255)                                                                                                                  cttcACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCT‑GGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTg                              <  4:193600/148‑1 (MQ=255)                                                                                                                  cttcACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCT‑GGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTg                              <  3:250241/148‑1 (MQ=255)                                                                                                                        tGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCT‑GGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGTAagag                       >  2:317035/1‑149 (MQ=255)                                                                                                                          cccAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCT‑GGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGTAAGAgg                      <  8:224867/148‑1 (MQ=255)                                                                                                                              gAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCT‑GGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGTAAGAGGGTAGc                 <  1:317035/149‑1 (MQ=255)                                                                                                                                  ttGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCT‑GGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGTAAGAgg                      <  8:168476/140‑1 (MQ=255)                                                                                                                                  ttGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCT‑GGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGTAAGAgg                      >  7:168476/1‑140 (MQ=255)                                                                                                                                    gTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCT‑GGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGTAAGAGGGTAGCCATAGc           >  5:170132/1‑149 (MQ=255)                                                                                                                                     tCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCT‑GGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGTAAGAGGGTAGCCATAGCt          >  3:251696/1‑149 (MQ=255)                                                                                                                                       gggCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCT‑GGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGTAAGAGGGTAGCCATAGCTgg        <  3:267608/149‑1 (MQ=255)                                                                                                                                        ggCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGTGTGATTTCATCGTCGCT‑GGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGTAAGAGGGTAGCCATAGCTGGc       <  6:103042/149‑1 (MQ=255)                                                                                                                                         gCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTAAAGCAGGACAGGGCGGTTGTCCGGCGGGGTGGTGGGGATGGTGGCCTCGAACGGGTTGGTAAGAGGGTAGCCATTGCTgg        >  1:279557/1‑148 (MQ=255)                                                                                                                                         gCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCT‑GGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGTAaga                        >  3:139657/1‑131 (MQ=255)                                                                                                                                         gCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCT‑GGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGTAaga                        <  4:139657/131‑1 (MQ=255)                                                                                                                                             cGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCT‑GGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGTAAGAGGGTAGCCATAGCTGGCAGCGt  >  3:13429/1‑149 (MQ=255)                                                                                                                                             cGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCT‑GGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGTAAGAGGGTAGCCATAGCTGGCAGCGt  <  4:251696/149‑1 (MQ=255)                                                                                                                                               tGAACTTGCGCCCGGAGCGCGAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCT‑GGGCAGCTCAGGgcgc                                                                        <  6:258299/77‑1 (MQ=37)                                                                                                                                               tGAACTTGCGCCCGGAGCGCGAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCT‑GGGCAGCTCAGGgcgc                                                                        >  5:258299/1‑77 (MQ=37)                                                                                                                                                gAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCT‑GGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGTtggtg                                             <  5:248707/103‑1 (MQ=255)                                                                                                                                                gAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGACGGACTGAGTCATTTCATCGTCGCT‑GGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGAtggtg                                             >  6:248707/1‑103 (MQ=255)                                                                                                                                                 aaCTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCT‑GGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGTAAGAGGGTAGCCATAGCTgg        <  2:32042/139‑1 (MQ=255)                                                                                                                                                 aaCTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCT‑GGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGTAAGAGGGTAGCCATAGCTgg        >  1:32042/1‑139 (MQ=255)                                                                                                                                                     |                                                                                                                                              AACAGCTTTTTCAGGTATTCATTGATGCTGCTGCTGTAGGGCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGGGCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCT‑GGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGTAAGAGGGTAGCCATAGCTGGCAGCGT  >  NC_002947/1938691‑1938979

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