Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,786,024	T→A	13.1%	intergenic (+40/‑146)	rpsB → / → tsf	30S ribosomal protein S2/elongation factor Ts

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,786,024	0	T	A	13.1%	155.8 / 13.4	61	intergenic (+40/‑146)	rpsB/tsf	30S ribosomal protein S2/elongation factor Ts
Reads supporting (aligned to +/- strand):  ref base T (25/28);  new base A (4/4);  total (29/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 = 9.80e-01

AGCTGTACATGACTTCGATGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTGAAAAACTTTCCAAGAGGATTTTGAAA  >  NC_002947/1785878‑1786170                                                                                                                                                   |                                                                                                                                                   aGCTGTACATGACTTCGATGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGc                                                                                                                                                  >  7:258391/1‑149 (MQ=255)     gTACATGACTTCGATGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGccccc                                                                                                                                              >  3:244430/1‑149 (MQ=255)         aTGACTTCGATGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCtttt                                                                                                                                          >  6:295165/1‑149 (MQ=255)          tGACTTCGATGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCAAAGCCCCCttttt                                                                                                                                         <  3:82203/149‑1 (MQ=255)          tGACTTCGATGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCAAAGCCCCCttttt                                                                                                                                         <  4:233832/149‑1 (MQ=255)          tGACTTCGATGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCAAAGCCCCCttttt                                                                                                                                         <  6:36588/149‑1 (MQ=255)                 gATGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTATTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCAcc                                                                                                                                   <  7:313688/148‑1 (MQ=255)                 gATGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTATTTGGCACGCAAAAAGCGGGCTTTGCCCCCTTTTTGCCAcc                                                                                                                                   <  7:308486/148‑1 (MQ=255)                       tGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCtt                                                                                                                                 <  5:90561/144‑1 (MQ=255)                       tGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCtt                                                                                                                                 >  6:90561/1‑144 (MQ=255)                         aCGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATcc                                                                                                                          >  3:17738/1‑149 (MQ=255)                          cGCACTCCTCCGCGGCCGCAACCACGTCGCCGGCGGCACCGAACTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCct                                                                                                                         <  8:258391/149‑1 (MQ=255)                                 aTCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGt                                                                                                                    >  5:373285/1‑147 (MQ=255)                                             aacaacGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCg                                                                                                      >  2:323557/1‑149 (MQ=255)                                             aacaacGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCg                                                                                                      >  1:792/1‑149 (MQ=255)                                                aacGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGca                                                                                                    >  1:244418/1‑148 (MQ=255)                                                 acGTTGCCGGCTGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGcatc                                                                                                  <  6:185045/149‑1 (MQ=255)                                                   gTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGcatcat                                                                                                <  2:8555/149‑1 (MQ=255)                                                       ccGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGc                                                                                            >  1:334279/1‑149 (MQ=255)                                                          gcggcACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCctg                                                                                         <  8:203425/149‑1 (MQ=255)                                                               aCCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCtga                                                                                     >  3:111110/1‑148 (MQ=255)                                                               aCCGAAGTATACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCtga                                                                                     >  4:272605/1‑148 (MQ=255)                                                                 cGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGc                                                                                            <  1:135058/139‑1 (MQ=255)                                                                 cGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGc                                                                                            >  2:135058/1‑139 (MQ=255)                                                                     gTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCtgaaatgaaa                                                                              <  3:326542/149‑1 (MQ=255)                                                                      tttACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGaaaa                                                                             <  6:203442/149‑1 (MQ=255)                                                                          cGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGc                                                                         <  1:323557/149‑1 (MQ=255)                                                                           gCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCctgct                                                                                                                      <  6:1551/103‑1 (MQ=255)                                                                           gCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCctgct                                                                                                                      >  5:1551/1‑103 (MQ=255)                                                                           gCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCg                                                                        >  7:184621/1‑149 (MQ=255)                                                                               aagaagCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGAtt                                                                    <  2:792/149‑1 (MQ=255)                                                                               aagaagCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGAt                                                                     <  6:177864/148‑1 (MQ=255)                                                                                     cGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCa                                                                <  1:134683/147‑1 (MQ=255)                                                                                     cGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCa                                                                >  2:134683/1‑147 (MQ=255)                                                                                      ggCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGaa                                                             >  4:223606/1‑149 (MQ=255)                                                                                            aCCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAAc                                                       <  5:295165/149‑1 (MQ=255)                                                                                                 ctgctgAGTAATAGGCGCTATCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGccc                                                   <  6:347508/148‑1 (MQ=255)                                                                                                    ctgAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGa                                               <  3:272605/149‑1 (MQ=255)                                                                                                     tgAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGAt                                              >  5:143625/1‑149 (MQ=255)                                                                                                       aGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGAtt                                             <  4:244430/148‑1 (MQ=255)                                                                                                       aGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGAt                                              <  2:244418/147‑1 (MQ=255)                                                                                                        gTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGa                                           <  4:17738/149‑1 (MQ=255)                                                                                                              ggCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGt                                      <  6:373285/148‑1 (MQ=255)                                                                                                                     gCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATg                                <  6:143625/147‑1 (MQ=255)                                                                                                                      cGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGtt                             >  3:52356/1‑149 (MQ=255)                                                                                                                       gTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGtt                             >  7:8068/1‑148 (MQ=255)                                                                                                                           aCTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGt                              <  3:3134/143‑1 (MQ=255)                                                                                                                           aCTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGt                              >  4:3134/1‑143 (MQ=255)                                                                                                                            cTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCTGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGTATGGTTGaaaaa                       <  6:40469/149‑1 (MQ=255)                                                                                                                            cTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGAt                                              <  4:347058/126‑1 (MQ=255)                                                                                                                            cTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGAt                                              >  3:347058/1‑126 (MQ=255)                                                                                                                               ggCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTGAAAAACtt                    >  2:229497/1‑149 (MQ=255)                                                                                                                                gCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTGAAAAACttt                   <  1:229497/149‑1 (MQ=255)                                                                                                                                  aCGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTGAAAAACTTTc                  >  8:129036/1‑148 (MQ=255)                                                                                                                                      aaaaaGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTAAGCATCGGCCCAGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAAAGCCCGTGATTGACGGGGGGAATGGTTGAAAAACTTTCCaa               >  3:248075/1‑147 (MQ=255)                                                                                                                                      aaaaaGGGGGCAAAGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGtt                             >  5:26201/1‑133 (MQ=255)                                                                                                                                      aaaaaGGGGGCAAAGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGtt                             <  6:26201/133‑1 (MQ=255)                                                                                                                                      aaaaaGGGGGCAAAGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTGAAAAACTTTCCaa               >  5:71699/1‑147 (MQ=255)                                                                                                                                      aaaaaGGGGGCAAAGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTGAAAAACTTTCCaa               >  4:329683/1‑147 (MQ=255)                                                                                                                                      aaaaaGGGGGCAAAGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGAAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTGAAAAACTTTCCaa               >  5:77624/1‑147 (MQ=255)                                                                                                                                             gggCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTGAAAAACTTTCCAAGAGGAtttt      <  2:334279/149‑1 (MQ=255)                                                                                                                                                   tGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTGAAAAACTTTCCAAGAGGATTTTGaaa  <  8:184621/147‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                   AGCTGTACATGACTTCGATGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTGAAAAACTTTCCAAGAGGATTTTGAAA  >  NC_002947/1785878‑1786170

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