Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,539,013	T→G	6.2%	K658T (AAG→ACG)	PP_4872 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,539,013	0	T	G	6.2%	203.6 / 3.0	64	K658T (AAG→ACG)	PP_4872	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (35/25);  new base G (2/2);  total (37/27)
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.92e-01

GCGCTCTGATAAGTATCGGTCTGCGGCTGGTGCAGCGGGATCAGTTGGTGCTTGCCCGGCCCGATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAACGCCTGTACCTGGTCGGCGC  >  NC_002947/5538869‑5539160                                                                                                                                                  |                                                                                                                                                    gcgcTCTGATAAGTATCGGTGTGCGGGTGGTGCAGCGGGATCAGTTGGTGCTTGCGCGGCCCGATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTtgtg                                                                                                                                                 <  5:267446/149‑1 (MQ=255)         aTAAGTATCGGTCTGCGGCTGGTGCAGCGGGATCAGTTGGTGCTTGCCCGGCCCGATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGc                                                                                                                                         >  7:17262/1‑149 (MQ=255)              tATCGGTCTGCGGCTGGTGCAGCGGGATTAGTTGGTGCTTGCCCGGCCCGATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGc                                                                                                                                    >  2:143069/1‑149 (MQ=255)              tATCGGTCTGCGGCTGGTGCAGCGGGATCAGTTGGTGCTTGCCCGGCCCGATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGc                                                                                                                                    >  8:48847/1‑149 (MQ=255)               aTCGGTCTGCGGCTGGTGCAGCGGGATCAGTTGGTGCTTGCCCGGCCCGATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGc                                                                                                                                    <  2:274598/148‑1 (MQ=255)                tCGGTCTGCGGCTGGTGCAGCGGGATCAGTTGGTGCTTGCCCGGCCCGATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCt                                                                                                                                   <  1:136412/148‑1 (MQ=255)                    tCTGCGGCTGGTGCAGCGGGATCAGTTGGTGCTTGCCCGGCCCGATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCgcgtcg                                                                                                                              >  4:298772/1‑149 (MQ=255)                          gCTGGTGCAGCGGGATCAGTTGGTGCTTGCCCGGCCCGATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCtt                                                                                                                         >  6:107167/1‑148 (MQ=255)                          gCTGGTGCAGCGGGATCAGTTGGTGCTTGCCCGGCCCGATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAAC‑GCCGTCTGCACAAGGCGttctt                                                                                                                                   <  4:318045/138‑4 (MQ=255)                          gCTGGTGCAGCGGGATCAGTTGGTGCTTGCCCGGCCCGATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAAC‑GCCGTCTGCACAAGGCGttctt                                                                                                                                   >  3:318045/1‑135 (MQ=255)                             ggTGCAGCGGGATCAGTTGGTGCTTGCCCGGCCCGATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGAc                                                                                                                      >  3:13427/1‑148 (MQ=255)                                 cAGCGGGATCAGTTGGTGCTTGCCCGGCCCGATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGAtc                                                                                                                  >  8:116860/1‑148 (MQ=255)                                  aGCGGGATCAGTTGGTGCTTGCCCGGCCCGATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGTtctc                                                                                                                >  1:194631/1‑149 (MQ=255)                                   gCGGGATCAGTTGGTGCTTGCCCGGCCCGATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGAtctc                                                                                                                >  4:327062/1‑148 (MQ=255)                                            gTTGGTGCTTGCCCGGCCCGATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCg                                                                                                      >  8:28384/1‑149 (MQ=255)                                                gTGCTTGCCCGGCCCGATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGctc                                                                                                   <  2:13461/148‑1 (MQ=255)                                                   cTTGCCCGGCCCGATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGt                                                                                               <  5:235322/149‑1 (MQ=255)                                                    ttGCCCGGCCCGATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTg                                                                                                >  2:66581/1‑147 (MQ=255)                                                      gcccggcccgATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGTCCAGCCCTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTa                                                                                              <  3:207959/147‑1 (MQ=255)                                                      gcccggcccgATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAg                                                                                             <  4:188021/148‑1 (MQ=255)                                                               gATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGctc                                                                                                   >  6:75019/1‑133 (MQ=255)                                                                aTCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCAcc                                                                                       <  3:180707/144‑1 (MQ=255)                                                                aTCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCAcc                                                                                       >  4:180707/1‑144 (MQ=255)                                                                      tcggctcggcCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGTTCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCgg                                                                            <  2:194631/149‑1 (MQ=255)                                                                          ctcggcCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCgg                                                                            <  6:54517/145‑1 (MQ=255)                                                                          ctcggcCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCgg                                                                            >  5:54517/1‑145 (MQ=255)                                                                          ctcggcCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTc                                                                        <  8:17262/149‑1 (MQ=255)                                                                             ggcCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCtt                                                                      <  7:48847/148‑1 (MQ=255)                                                                              gcCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGc                                                                    <  3:327062/149‑1 (MQ=255)                                                                               cccATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCACTTGTAGGTAACCTTGCGCAGCGGGTTCTTGcc                                                                   >  6:148766/1‑149 (MQ=255)                                                                                  aTACGCTGCAGCGCTTCACGCAGCATCGGCGAGCCCTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCt                                                                  >  1:30796/1‑147 (MQ=255)                                                                                  aTACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACTGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCt                                                                  >  1:267324/1‑148 (MQ=255)                                                                                   tACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAg                                                               >  6:73817/1‑149 (MQ=255)                                                                                   tACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAg                                                               >  7:62866/1‑149 (MQ=255)                                                                                                ttCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACAtggccg                                                  >  4:312729/1‑149 (MQ=255)                                                                                                  cACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACAtggccgt                                                 <  4:13427/148‑1 (MQ=255)                                                                                                     gcagcaTCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAg                                                                                             >  6:178657/1‑101 (MQ=255)                                                                                                     gcagcaTCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACAtggccgtggcc                                             >  6:162427/1‑149 (MQ=255)                                                                                                        gcaTCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAg                                          <  5:162427/149‑1 (MQ=255)                                                                                                        gcaTCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAg                                          <  3:312729/149‑1 (MQ=255)                                                                                                           tCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCTCGGCCACGGcc                                       <  2:30796/149‑1 (MQ=255)                                                                                                              gCCAGCCCTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCAtt                                    <  8:62866/149‑1 (MQ=255)                                                                                                                cAGCCCTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTgg                                  >  5:127588/1‑149 (MQ=255)                                                                                                                   cccTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGa                               >  6:341480/1‑149 (MQ=255)                                                                                                                    ccTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGAc                              >  1:86921/1‑149 (MQ=255)                                                                                                                    ccTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGAc                              <  5:341480/149‑1 (MQ=255)                                                                                                                     cTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACg                             >  4:70998/1‑149 (MQ=255)                                                                                                                      ttCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCAttt                                   <  1:224035/142‑2 (MQ=255)                                                                                                                      ttCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCATTCGCTCTTGTAGGTCACCTTGAGAAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTg                                   >  2:224035/1‑142 (MQ=255)                                                                                                                        cGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACgg                            <  4:304492/147‑1 (MQ=255)                                                                                                                         gggTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGt                          <  5:73817/148‑1 (MQ=255)                                                                                                                           gTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCa                                      >  2:312945/1‑134 (MQ=255)                                                                                                                           gTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCa                                      <  1:312945/134‑1 (MQ=255)                                                                                                                                 gTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAACGCCt                 <  2:86921/149‑1 (MQ=255)                                                                                                                                      gcAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAACGCCTGTAcc            <  7:28384/149‑1 (MQ=255)                                                                                                                                      cgagcaaC‑GCCGTCTGCACAAGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACAtggccgtg                                                <  3:268735/108‑1 (MQ=255)                                                                                                                                      cgagcaaC‑GCCGTCTGCACAAGGCCTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACAtggccgtg                                                >  4:268735/6‑113 (MQ=255)                                                                                                                                       cAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAACGCCTGTACCt           >  5:328950/1‑149 (MQ=255)                                                                                                                                          gAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAACGCCTGTACCTGGt        >  6:12437/1‑149 (MQ=255)                                                                                                                                           aaC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAACGCCTGTACCTGGTc       >  5:80988/1‑149 (MQ=255)                                                                                                                                           aaC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAACGCCTGTACCTGGTc       >  4:202527/1‑149 (MQ=255)                                                                                                                                                 cTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAACGCCTGTACCTGGTCGgcgc  >  2:151500/1‑149 (MQ=255)                                                                                                                                                 cTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAACGCCTGTACCTGGTCGgcgc  >  4:189671/1‑149 (MQ=255)                                                                                                                                                 cTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTAGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAACGCCTGTACCTGGTCGgcgc  >  7:166132/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                                    GCGCTCTGATAAGTATCGGTCTGCGGCTGGTGCAGCGGGATCAGTTGGTGCTTGCCCGGCCCGATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAAC‑GCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAACGCCTGTACCTGGTCGGCGC  >  NC_002947/5538869‑5539160

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