Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,622,933	T→A	6.8%	M269L (ATG→TTG)	PP_4941 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,622,933	0	T	A	6.8%	198.5 / 3.2	69	M269L (ATG→TTG)	PP_4941	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (30/34);  new base A (2/3);  total (32/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.59e-01

AAGGCTCGCAGGTCGGCCAGCAGGTTGTGCCAGGTGGGTGCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGAT‑ATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGCGGCCTCGGCCTTTGCCAGGCCCATCTGCTGG  >  NC_002947/5622785‑5623078                                                                                                                                                      |                                                                                                                                                  aaGGCTCGCAGGTCGGCCAGCAGGTTGTGCCAGGTGGGTGCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGAT‑ATCGGCCAGGTCGACTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCAt                                                                                                                                                   <  6:384939/149‑1 (MQ=255)     cTCGCAGGTCGGCCAGCAGGTTGTGCCAGGTGGGTGCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGAT‑ATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATggca                                                                                                                                               <  5:37570/149‑1 (MQ=255)       cGCAGGTCGGCCAGCAGGTTGTGCCAGGTGGGTGCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGAT‑ATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATggcagg                                                                                                                                             >  3:341334/1‑149 (MQ=255)        gCAGGTCGGCCAGCAGGTTGTGCCAGGTGGGTGCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGAT‑ATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCCAtgcaggc                                                                                                                                            <  7:96491/149‑1 (MQ=255)        gCAGGTCGGCCAGCAGGTTGTGCCAGGTGGGTGCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGAT‑ATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATggcagg                                                                                                                                             >  6:373709/1‑148 (MQ=255)        gCAGGTCGGCCAGCAGGTTGTGCCAGGTGGGTGCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGAT‑ATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATggcagg                                                                                                                                             >  4:78799/1‑148 (MQ=255)        gCAGGGCGGCCAGCAGGTTGTGCCAGGTGGGTGCGCCATCGCGGTCGCCGGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGAT‑ATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATggcag                                                                                                                                              <  3:71327/147‑1 (MQ=255)          aGGTCGGCCAGCAGGTTGTGCCAGGTGGGTGCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGAT‑ATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCCATgcaggcag                                                                                                                                          <  6:84998/149‑1 (MQ=255)          aGGTCGGCCAGCAGGTTGTGCCAGGTGGGTGCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGAT‑ATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCCATgcaggcag                                                                                                                                          <  5:213010/149‑1 (MQ=255)               ggCCAGCAGGTTGTGCCAGGTGGGTGCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGAT‑ATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATggcagg                                                                                                                                             >  5:60703/1‑141 (MQ=255)               ggCCAGCAGGTTGTGCCAGGTGGGTGCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGAT‑ATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCa                                                                                                                                     <  5:350081/149‑1 (MQ=255)               ggCCAGCAGGTTGTGCCAGGTGGGTGCGCCATCGCGGTCGCCCGGTAACGGGAACGGGGTGAACTGGCGGAACACGCGGAT‑ATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATggcagg                                                                                                                                             <  6:60703/141‑1 (MQ=255)                  cagcagGTTGTGCCAGGTGGGTGCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGAT‑ATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAgg                                                                                                                                   >  7:324205/1‑148 (MQ=255)                        gTTGTGCCAGGTGGGTGCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGAGGAACACGCGGAT‑ATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAg                                                                                                                             >  7:58197/1‑148 (MQ=255)                                     ggtgCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGAT‑ATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTAt                                                                                                               <  2:251238/149‑1 (MQ=255)                                      gtgCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGAT‑ATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCCGCTGCAGGCAAAAGTAGCCCAGCTGTATg                                                                                                              >  1:369631/1‑149 (MQ=255)                                          gcCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGAT‑ATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAAc                                                                                                          <  5:373709/149‑1 (MQ=255)                                           ccATCGCGGTCGCCCGGTTACGGGAACGGGTTGAACTGGCGGAACACGCGGAT‑ATCGTCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCa                                                                                                            <  1:164106/146‑1 (MQ=255)                                           ccATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGAT‑ATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCa                                                                                                            >  2:164106/1‑146 (MQ=255)                                                gcgGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGAT‑ATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGAc                                                                                                     <  5:194467/148‑1 (MQ=255)                                                   gTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGAT‑ATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCg                                                                                                 >  8:221885/1‑149 (MQ=255)                                                      gCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGAT‑ATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCa                                                                                              <  8:157994/149‑1 (MQ=255)                                                      gCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGAT‑ATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCa                                                                                              >  6:123246/1‑149 (MQ=255)                                                      gCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGAT‑ATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCa                                                                                              >  6:9965/1‑149 (MQ=255)                                                      gCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGAT‑ATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCa                                                                                              <  4:13750/149‑1 (MQ=255)                                                      gCCCGGTAACGGGAACGGGGTGAACTGGCGGAACACGCGGAT‑ATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCa                                                                                              <  8:162987/149‑1 (MQ=255)                                                        ccGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGAT‑ATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACg                                                                                            >  7:272743/1‑149 (MQ=255)                                                           gTAACGGGAACGGGTTGAACTGGCGGAACACGCGGAT‑ATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCg                                                                                         >  4:305987/1‑149 (MQ=255)                                                            tAACGGGAACGGGTTGAACTGGCGGAACACGCGGATGATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCg                                                                                         >  1:63780/1‑149 (MQ=255)                                                                   aacgggTTGAACTGGCGGAACACGCGGAT‑ATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAAcg                                                                                 >  6:284879/1‑149 (MQ=255)                                                                       ggTTGAACTGGCGGAACACGCGGAT‑ATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCgg                                                                              >  3:17549/1‑148 (MQ=255)                                                                         ttGAACTGGCGGAACACGCGGAT‑ATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCac                                                                           >  8:69173/1‑149 (MQ=255)                                                                         ttGAACTGGCGGAACACGCGGAT‑ATCGGCCAGGTCGGCTTCACGCGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTc                                                                                                   >  5:300402/1‑125 (MQ=255)                                                                         ttGAACTGGCGGAACACGCGGAT‑ATCGGCCAGGTCGGCTTCACGCGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTc                                                                                                   <  6:300402/125‑1 (MQ=255)                                                                           gAACTGGCGGAACACGGGGAT‑ATCGGCGAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGAAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCGAAcgcg                                                                               <  7:224939/143‑1 (MQ=255)                                                                           gAACTGGCGGAACACGCGGAT‑ATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTAAGCAACGCGACACCCGCACGCCCGCCCAAcgcg                                                                               >  8:224939/1‑143 (MQ=255)                                                                             aCTGGCGGAACACGCGGAT‑ATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCacgac                                                                        >  5:239990/1‑148 (MQ=255)                                                                                   ggAACACGCGGAT‑ATCGGCCAGGTCGCCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCAGCAAGGCAGGCAGCTGCACGCAAAAGTAGCCCAGATGTATGCACCGCGACGCCGGCACGCCGGCCAAACGCGTCaccacgcttct                                                                  >  5:316600/1‑143 (MQ=255)                                                                                          gcgGAT‑ATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAg                                                                                                                                    >  5:57375/1‑75 (MQ=255)                                                                                          gcgGAT‑ATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAg                                                                                                                                    <  6:57375/75‑1 (MQ=255)                                                                                          gcgGAT‑ATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAgg                                                          >  4:368805/1‑149 (MQ=255)                                                                                            ggAT‑ATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCa                                                        >  7:90080/1‑149 (MQ=255)                                                                                                 atCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACg                                                      >  2:204855/1‑147 (MQ=255)                                                                                                   cGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAgg                                                  >  2:137846/1‑149 (MQ=255)                                                                                                    ggCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGgc                                                 >  2:254046/1‑149 (MQ=255)                                                                                                        aGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTTCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGTCGGCCCAACGCGGTACGACGATGCTGTCCCAGGCACGCAGGCGccc                                             <  8:58197/149‑1 (MQ=255)                                                                                                           tCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCtt                                          <  1:204855/149‑1 (MQ=255)                                                                                                             ggCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCa                                        <  4:351402/149‑1 (MQ=255)                                                                                                             ggCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCGGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCa                                        <  7:214215/149‑1 (MQ=255)                                                                                                               cTTCACGTGAGGCGGCCGGCTGGTCGGGGGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAgg                                      <  4:318586/149‑1 (MQ=255)                                                                                                                  cACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGgcgc                                   >  1:206436/1‑149 (MQ=255)                                                                                                                   aCGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGgcgcg                                  <  7:221885/149‑1 (MQ=255)                                                                                                                    cGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGgcgcgc                                 >  1:379584/1‑149 (MQ=255)                                                                                                                    cGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGgcgcgc                                 <  1:137846/149‑1 (MQ=255)                                                                                                                    cGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGgcgcgc                                 <  8:324205/149‑1 (MQ=255)                                                                                                                           cggccggcTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGcggcctcg                          <  5:123246/149‑1 (MQ=255)                                                                                                                           cggccggcTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGcggcctcg                          <  6:239990/149‑1 (MQ=255)                                                                                                                             gccggcTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGcggcctcggc                        <  4:341334/149‑1 (MQ=255)                                                                                                                             gccggcTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGcggcctcggc                        <  8:272743/149‑1 (MQ=255)                                                                                                                             gccggcTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGcggcctcggc                        <  3:368805/149‑1 (MQ=255)                                                                                                                                 gcTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGCGGCCTCGGCCttt                    <  8:90080/149‑1 (MQ=255)                                                                                                                                 gcTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGCGGCCTCGGCCttt                    <  2:350083/149‑1 (MQ=255)                                                                                                                                    ggTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGCGGCCTCGGCCTTTGcc                 <  5:9965/149‑1 (MQ=255)                                                                                                                                        gggTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGCGGCCTCGGCCTTTGCCAGGc             <  7:69173/149‑1 (MQ=255)                                                                                                                                         ggTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGCGGCCTCGGCCTTTGCCAGGcc            <  3:305987/149‑1 (MQ=255)                                                                                                                                          gTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGCGGCCTCGGCCTTTGCCAGGcc            >  8:345837/1‑148 (MQ=255)                                                                                                                                             ctgcctgcCATGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGCGGCCTCGGCCTTTGCCAGGCCCATc        >  3:172523/1‑149 (MQ=255)                                                                                                                                               gcctgcATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGCGGCCTCGGCCTTTGCCAGGCCCATctg      >  2:269767/1‑149 (MQ=255)                                                                                                                                                    cATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGCGGCCTCGGCCTTTGCCAGGCCCATctgctg   <  3:78799/147‑1 (MQ=255)                                                                                                                                                    cATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGCGGCCTCGGCCTTTGCCAGGCCCATCTGCTgg  <  7:345837/148‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                  AAGGCTCGCAGGTCGGCCAGCAGGTTGTGCCAGGTGGGTGCGCCATCGCGGTCGCCCGGTAACGGGAACGGGTTGAACTGGCGGAACACGCGGAT‑ATCGGCCAGGTCGGCTTCACGTGAGGCGGCCGGCTGGTCGGGTGCTGCCTGCATGGCAGGCAGCTGCAGGCAAAAGTAGCCCAGCTGTATGCAACGCGACTCCGGCACGCCGGCCCAACGCGGCACGACGATGCTGTCCCAGGCACGCAGGCGCCCCTTCAGGCGCGCGGCCTCGGCCTTTGCCAGGCCCATCTGCTGG  >  NC_002947/5622785‑5623078

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑