Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,111,977	T→A	8.3%	H118Q (CAT→CAA)	PP_5362 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	6,111,977	0	T	A	8.3%	212.2 / 7.9	72	H118Q (CAT→CAA)	PP_5362	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (32/34);  new base A (4/2);  total (36/36)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01

GGGCACAAACCAACAGCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCT‑GCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAACCCTTTGTCAGACCGAAAAGCTATGCCGTAAC  >  NC_002947/6111829‑6112124                                                                                                                                                     |                                                                                                                                                     gggCACAAACCAACAGCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCAt                                                                                                                                                      >  8:165238/1‑149 (MQ=255)    cacaAACCAACAGCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGTGGGGAGCGCTTCAAGCTGCATTg                                                                                                                                                    >  7:253146/1‑148 (MQ=255)     acaAACCAACAGCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTg                                                                                                                                                    >  3:54297/1‑147 (MQ=255)     acaAACCAACAGCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCAATg                                                                                                                                                    >  7:110218/1‑147 (MQ=255)       aaaCCAACAGCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGc                                                                                                                                                >  6:271173/1‑149 (MQ=255)       aaaCCAACAGCCGCTATGGCACATTGACTTTTTTTACCGGCCTGACCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGc                                                                                                                                                >  3:45967/1‑149 (MQ=255)        aaCCAACAGCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCt                                                                                                                                               <  4:54297/149‑1 (MQ=255)         aCCAACAGCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCAGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCtt                                                                                                                                              >  2:66070/1‑149 (MQ=255)            atcAGCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCAATGCAGCTTGaa                                                                                                                                           <  8:110218/147‑2 (MQ=255)            aaCAGCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCAATGCAGCTTGaa                                                                                                                                           <  5:146754/149‑2 (MQ=255)               aGCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATc                                                                                                                                          >  8:281650/1‑147 (MQ=255)                gCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGTGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAg                                                                                                                                       <  8:253146/149‑1 (MQ=255)                gCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAg                                                                                                                                       <  6:308483/149‑1 (MQ=255)                 ccGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAgg                                                                                                                                      <  6:124534/149‑1 (MQ=255)                 ccGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAgg                                                                                                                                      >  2:181115/1‑149 (MQ=255)                 ccGCTATGGCAAATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAgg                                                                                                                                      <  8:225101/149‑1 (MQ=255)                  cGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAggg                                                                                                                                     >  6:181947/1‑149 (MQ=255)                    cTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGt                                                                                                                                   <  4:232432/149‑1 (MQ=255)                       tGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTg                                                                                                                                  <  5:154873/147‑1 (MQ=255)                          cacaTTGACTTTTTTTACCGGCCTGACCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGtt                                                                                                                             <  4:45967/149‑1 (MQ=255)                           acaTTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTg                                                                                                                            <  2:298514/149‑1 (MQ=255)                                 aCTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGcct                                                                                                                      >  5:143537/1‑149 (MQ=255)                                  cTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACTACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGcctc                                                                                                                     >  5:142517/1‑149 (MQ=255)                                  cTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGcctc                                                                                                                     >  8:261660/1‑149 (MQ=255)                                  cTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGcctc                                                                                                                     <  6:310356/149‑1 (MQ=255)                                   tttttttACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGcctcc                                                                                                                    >  2:147894/1‑149 (MQ=255)                                         tACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCa                                                                                                               <  7:281650/148‑1 (MQ=255)                                           ccGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCAttt                                                                                                            <  5:179584/149‑1 (MQ=255)                                            cGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCAt                                                                                                              <  8:37745/146‑1 (MQ=255)                                            cGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCAt                                                                                                              >  7:37745/1‑146 (MQ=255)                                             ggCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTg                                                                                                                                             <  5:240345/114‑1 (MQ=255)                                             ggCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTg                                                                                                                                             >  6:240345/1‑114 (MQ=255)                                              gCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTg                                                                                                          >  1:318281/1‑148 (MQ=255)                                                     cAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGtttt                                                                                                  >  4:192299/1‑149 (MQ=255)                                                      aGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGttttt                                                                                                 <  3:192299/149‑1 (MQ=255)                                                       gCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGtttttt                                                                                                <  7:261660/149‑1 (MQ=255)                                                         gATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTtct                                                                                              <  5:271173/149‑1 (MQ=255)                                                          aTGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTtctc                                                                                             >  1:131616/1‑149 (MQ=255)                                                             gccaGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTctgc                                                                                          <  5:339280/146‑1 (MQ=255)                                                                aGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTctgct‑g                                                                                       <  7:165238/148‑1 (MQ=255)                                                                           ggCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCT‑GCCTATTTTGTaa                                                                           <  2:131616/149‑1 (MQ=255)                                                                                cAGCCTACGCGGCGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCT‑GCCTATTTTGTaa                                                                           >  2:256954/1‑144 (MQ=255)                                                                                cAGCCTACGCGGCGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCT‑GCCTATTTTGTaa                                                                           <  1:256954/144‑1 (MQ=255)                                                                                     gccgcGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCT‑GCCTATTTTGt                                                                             >  4:226941/1‑137 (MQ=255)                                                                                     gccgcGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCT‑GCCTATTTTGt                                                                             <  3:226941/137‑1 (MQ=255)                                                                                             ggCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCT‑GCCTATTTTGTAAGTCAGCTTCGAGCAcgc                                                          <  1:181115/148‑1 (MQ=255)                                                                                               cTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCT‑GCCTATTTTGTAAGTCAGCTTCGAGCACGCGCt                                                       <  6:143537/149‑1 (MQ=255)                                                                                                tGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGtttttt                                                                                                <  6:138743/108‑1 (MQ=255)                                                                                                tGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGtttttt                                                                                                >  5:138743/1‑108 (MQ=255)                                                                                                tGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCT‑GCCTATTTTGTAAGTCAGCTTCGAGCACGCGCt                                                       >  5:333295/1‑148 (MQ=255)                                                                                                 ggACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCT‑GCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGa                                                     <  4:14086/149‑1 (MQ=255)                                                                                                    ccaccTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCT‑GCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGaaaa                                                  <  2:318281/149‑1 (MQ=255)                                                                                                     caccTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCT‑GCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGaaaa                                                  >  8:59226/1‑148 (MQ=255)                                                                                                     caccTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCT‑GCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAac                                                 <  5:181947/149‑1 (MQ=255)                                                                                                       ccTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCT‑GCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAacac                                               >  3:194792/1‑149 (MQ=255)                                                                                                       ccTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCT‑GCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAacac                                               >  4:171062/1‑149 (MQ=255)                                                                                                           gAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCGTTCAATTTGAACGTTTTTTCTCTGCTCGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCa                                             <  1:66070/148‑1 (MQ=255)                                                                                                                 cgTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGc                                                                                                                        >  1:319691/1‑67 (MQ=255)                                                                                                                 cgTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGc                                                                                                                        <  2:319691/67‑1 (MQ=255)                                                                                                                       cAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCT‑GCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAAcc                                <  6:142517/148‑1 (MQ=255)                                                                                                                             cAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCT‑GCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAACCCTTTGt                          <  6:333295/148‑1 (MQ=255)                                                                                                                              aGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCT‑GCCt                                                                                    >  6:304090/1‑89 (MQ=255)                                                                                                                              aGCGGGGAGCGCTTCAAGCTGAATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCT‑GCCt                                                                                    <  5:304090/89‑1 (MQ=255)                                                                                                                                         gatCAAGCTGCAATGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCT‑GCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAACCCTTTGTCAGACCGAAAAGc             >  2:297694/3‑149 (MQ=255)                                                                                                                                         cTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCT‑GCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAACCCTTTGTCAGACCGAAAAGc             >  5:120651/1‑149 (MQ=255)                                                                                                                                           tCAAGCTGCAATGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCT‑GCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAACACTTTGTCAGACCGAAAAGCTa           >  7:143698/1‑149 (MQ=255)                                                                                                                                            cAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCT‑GCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAACCCTTTGTCAGACCGAAAAGCTAt          >  4:62207/1‑149 (MQ=255)                                                                                                                                            cAAGCTGCAATGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCT‑GCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAACCCTTTGTCAGACCGAAAAGCTAt          >  4:231189/1‑149 (MQ=255)                                                                                                                                                cTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCT‑GCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAACCCTTTGTCAGACCGAAAAGCTATGCCg      >  4:77655/1‑149 (MQ=255)                                                                                                                                                 tGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCT‑GCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAACCCTTTGTCAGACCGAAAAGCTATGCCg      >  5:51240/1‑148 (MQ=255)                                                                                                                                                   cATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCT‑GCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAACCCTTTGTCAGACCGAAAAGCTATg         <  4:194051/143‑1 (MQ=255)                                                                                                                                                   cATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCT‑GCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAACCCTTTGTCAGACCGAAAAGCTATg         >  3:194051/1‑143 (MQ=255)                                                                                                                                                    aTTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCT‑GCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAACCCTTTGTCAGACCGAAAAGCTATGCCGTAAc  >  2:289164/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                     GGGCACAAACCAACAGCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCT‑GCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAACCCTTTGTCAGACCGAAAAGCTATGCCGTAAC  >  NC_002947/6111829‑6112124

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