Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,991,155	T→G	5.5%	L802* (TTA→TGA)	PP_1780 →	mannosyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,991,155	0	T	G	5.5%	271.4 / 5.4	89	L802* (TTA→TGA)	PP_1780	mannosyltransferase
Reads supporting (aligned to +/- strand):  ref base T (39/45);  new base G (2/3);  total (41/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.73e-01

AACCGGACAAGCATACCTTGGCCCTCAAGATGGGCGAACTGCGCGGCAAAATCGTCGAAGGCGACGATGCACTCGAGCAGC‑GCCGGCGCCAAGGCATT‑GAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGC‑GGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTGGGCAACCAAATGCGGCATTGCTGAATACACTCAGGCATTGCTCAGTGAGTTCAAGCAGGACGATT  >  NC_002947/1991008‑1991298                                                                                                                                                       |                                                                                                                                                aaCCGGACAAGCATACCTTGGCCCTCAAGATGGGCGAACTGCGCGGCAAAATCGTCGAAGGCGACGATGCACTCGAGCAGC‑GCCGGCGCCAAGGCATT‑GAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGC‑GGTTTTa                                                                                                                                                <  3:143605/149‑1 (MQ=255)  aCCGGACAAGCATACCTTGGCCCTCAAGATGGGCGAACTGCGCGGCAAAATCGTCGAAGGCGACGATGCACTCGAGCAGC‑GCCGGCGCCAAGGCATT‑GAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGC‑GGTTTTAg                                                                                                                                               >  7:100636/1‑149 (MQ=255)   ccGGACAAGCATACCTTGGCCCTCAAGATGGGCGAACTGCGCGGCAAAATCGTCGAAGGCGACGATGCACTCGAGCAGC‑GCCGGCGCCAAGGCATT‑GAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGC‑GGTTTTAGc                                                                                                                                              <  6:194006/149‑1 (MQ=255)       aCAAGCATACCTTGGCCCTCAAGATGCGCGAACTGCGCGGCAAAATCGTCGAAGGCGACGATGCACTCGAGCAGC‑GCCGGCGCCAAGGCATT‑GAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGC‑GGTTTTAGCCAGt                                                                                                                                          <  2:70926/149‑1 (MQ=255)           gCATACCTTGGCCCTCAAGATGGGCGAACTGCGCGGCAAAATCGTCGAAGGCGACGATGCACTCGAGCAGC‑GCCGGCGCCAAGGCATT‑GAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGC‑GGTTTTAGCCAGttgg                                                                                                                                       >  1:129496/1‑145 (MQ=255)             aTACCTTGGCCCTCAAGATGGGCGAACTGCGCGGCAAAATCGTCGAAGGCGACGATGCACTCGAGCAGC‑GCCGGCGCCAAGGCATT‑GAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGC‑GGTTTTAg                                                                                                                                               <  4:281563/138‑1 (MQ=255)             aTACCTTGGCCCTCAAGATGGGCGAACTGCGCGGCAAAATCGTCGAAGGCGACGATGCACTCGAGCAGC‑GCCGGCGCCAAGGCATT‑GAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGC‑GGTTTTAg                                                                                                                                               >  3:281563/1‑138 (MQ=255)             aTACCTTGGCCCTCAAGATGGGCGAACTGCGCGGCAAAATCGTCGAAGGCGACGATGCACTCGAGCAGC‑GCCGGCGCCAAGGCATT‑GAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGC‑GGTTTTAGCCAGTGCCAAc                                                                                                                                    <  3:18252/149‑1 (MQ=255)                     gCCCTCAAGATGGGCGAACTGCGCGGCAAAATCGTCGAAGGCGACGATGCACTCGAGCAGC‑GCCGGCGCCAAGGCATT‑GAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGC‑GGTTTTAGCCAGTGcc                                                                                                                                       <  2:33620/138‑1 (MQ=255)                     gCCCTCAAGATGGGCGAACTGCGCGGCAAAATCGTCGAAGGCGACGATGCACTCGAGCAGC‑GCCGGCGCCAAGGCATT‑GAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGC‑GGTTTTAGCCAGTGcc                                                                                                                                       >  1:33620/1‑138 (MQ=255)                      cccTCAAGATGGGCGAACTGCGCGGCAAAATCGTCGAAGGCGACGATGCACTCGAGCAGC‑GCCGGCGCCAAGGCATT‑GAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGC‑GGTTTTAGCCAGTTGGCACTGGCTaaaa                                                                                                                           <  2:129496/149‑1 (MQ=255)                       ccTCAAGATGGGCGAACTGCGCGGCAAAATCGTCGAAGGCGACGATGCACTCGAGCAGC‑GCCGGCGCCAAGGCATT‑GAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGC‑GGTTTTAGCCAGTGCCAACTGGCTCaaa                                                                                                                           >  1:180473/1‑148 (MQ=255)                       ccTCAAGATGGGCGAACTGCGCGGC‑AAATCGTCGAAGGCGACGATGCACTCGAGCAGC‑GCCGGCGCCAAGGCATT‑GAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGCTGGTTTTAGCCAGTGCCAACTGGCTCAAac                                                                                                                          <  1:31306/149‑1 (MQ=255)                         tCAAGATGGGCGAACTGCGCGGCAAAATCGTCGAAGGCGACGATGCACTCGAGCAGC‑GCCGGCGCCAAGGCATT‑GAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGC‑GGTTTTAGCCAGTGCCAACTGGCTCAAacac                                                                                                                        >  5:426724/1‑149 (MQ=255)                           aaGCTGGGCGAACTGCGCGGCAAAATCGTCGAAGGCGACGATGCACTCGAGCAGC‑GCCGGCGCCAAGGCATT‑GAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGC‑GGTTTGAGCCAGTTGGCACTGGCTAAAaccgc                                                                                                                       <  4:76318/148‑4 (MQ=255)                           aaGATGGGCGAACTGCGCGGCAAAATCGTCGAAGGCGACGATGCACTCGAGCAGC‑GCCGGCGCCAAGGCATT‑GAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGC‑GGTTTTAGCCAGTGCCAACTGGCTCAAacac                                                                                                                        <  1:413267/147‑1 (MQ=255)                             gATGGGCGAACTGCGCGGCAAAATCGTCGAAGGCGACGATGCACTCGAGCAGC‑GCCGGCGCCAAGGCATT‑GAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGC‑GGTTTTAGCCAGTGCCAACTGGCTCAAACACGaaa                                                                                                                    >  2:326583/1‑149 (MQ=255)                              aTGGGCGAACTGCGCGGCAAAATCGTCGAAGGCGACGATGCACTCGAGCAGC‑GCCGGCGCCAAGGCATT‑GAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGC‑GGTTTGAGCCAGTTGGCACTGGCTAAAaccgcctgc                                                                                                                   <  3:93821/149‑8 (MQ=38)                                gggCGAACTGCGCGGCAAAATCGTCGAAGGCGACGATGCACTCGAGCAGCTGCCGGCGCCAAGGCATT‑GAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGC‑GGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATg                                                                                                                  >  1:308328/1‑149 (MQ=255)                                 ggCGAACTGCGCGGCAAAATCGTCGAAGGCGACGATGCACTCGAGCAGC‑GCCGGCGCCAAGGCATT‑GAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGC‑GGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGcc                                                                                                                >  2:87302/1‑149 (MQ=255)                                 ggCGAACTGCGCGGCAAAATCGTCGAAGGCGACGATGCACTCGAGCAGC‑GCCGGCGCCAAGGCATT‑GAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGC‑GGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGc                                                                                                                 <  8:314200/148‑1 (MQ=255)                                 ggCGAACTGCGCGGCAAAATCGTCGAAGGCGACGATGCACTCGAGCAGC‑GCCGGCGCCAAGGCATT‑GAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGC‑GGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGc                                                                                                                 <  6:20400/148‑1 (MQ=255)                                          cgcgGCAAAATCGTCGAAGGCGACGATGCACTCGAGCAGC‑GCCGGCGCCAAGGCATT‑GAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGC‑GGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGgcc                                                                                                         <  6:102479/147‑1 (MQ=255)                                               cAAAATCGTCGAAGGCGACGATGCACTCGAGCAGC‑GCCGGCGCCAAGGCATT‑GAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGC‑GGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGgccgccgc                                                                                                    >  6:210371/1‑147 (MQ=255)                                                aaaaTCGTCGAAGGCGACGATGCACTCGAGCAGC‑GCCGGCGCCAAGGCATT‑GAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGC‑GGTTTTAGCCAGTGCCAACTGGGTCAAACACGAAATGCCGGCGGCCGCCGCCAg                                                                                                 >  1:355597/1‑149 (MQ=255)                                                aaaaTCGTCGAAGGCGACGATGCACTCGAGCAGC‑GCCGGCGCCAAGGCATT‑GAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGC‑GGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAg                                                                                                 <  8:100636/149‑1 (MQ=255)                                                 aaaTCGTCGAAGGCGACGATGCACTCGAGCAGC‑GCCGGCGCCAAGGCATT‑GAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGC‑GGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCcggcgg                                                                                                           <  3:55410/138‑1 (MQ=255)                                                 aaaTCGTCGAAGGCGACGATGCACTCGAGCAGC‑GCCGGCGCCAAGGCATT‑GAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGC‑GGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCcggcgg                                                                                                           >  4:55410/1‑138 (MQ=255)                                                        cgAAGGCGACGATGCACTCGAGCAGC‑GCCGGCGCCAAGGCATT‑GAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGC‑GGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGAcc                                                                                          >  1:341413/1‑148 (MQ=255)                                                             gCGACGATGCACTCGAGCAGC‑GCCGGCGCCAAGGCATT‑GAACTTGTTCGCAGCATCTACACTAGGGAGCAGAGTGCGCAGGC‑GGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACtg                                                                                    >  1:122639/1‑149 (MQ=255)                                                               gacgaTGCACTCGAGCAGC‑GCCGGCGCCAAGGCATT‑GAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGC‑GGTTTTAGCCAGTGCCAACTGGCTc                                                                                                                              >  3:337782/1‑105 (MQ=255)                                                               gacgaTGCACTCGAGCAGC‑GCCGGCGCCAAGGCATT‑GAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGC‑GGTTTTAGCCAGTGCCAACTGGCTc                                                                                                                              <  4:337782/105‑1 (MQ=255)                                                               gacgaTGCACTCGAGCAGC‑GCCGGCGCCAAGGCATT‑GAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGC‑GGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACtggt                                                                                  >  6:277550/1‑149 (MQ=255)                                                                  gaTGCACTCGAGCAGC‑GCCGGCGCCAAGGCATT‑GAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGC‑GGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACtggtggt                                                                               <  5:277550/149‑1 (MQ=255)                                                                  gaTGCACTCGAGCAGC‑GCCGGCGCCAAGGCATT‑GAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGC‑GGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACtggtggt                                                                               <  1:326583/149‑1 (MQ=255)                                                                      cACTCGAGCAGC‑GCCGGCGCCAAGGCATT‑GAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGC‑GGTTTTAGCCAGTGCCAACTGGCTCaaa                                                                                                                           >  8:279289/1‑101 (MQ=255)                                                                      cACTCGAGCAGC‑GCCGGCGCCAAGGCATT‑GAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGC‑GGTTTTAGCCAGTGCCAACTGGCTCaaa                                                                                                                           <  7:279289/101‑1 (MQ=255)                                                                       aCTCGAGCAGC‑GCCGGCGCCAAGGCATT‑GAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGC‑GGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGttt                                                                          >  5:125360/1‑149 (MQ=255)                                                                       aCTCGAGCAGC‑GCCGGCGCCAAGACATTAGAACTTGTTCGCAGAATCTACACTTGGAAGCAGAGTGCGAAGGC‑GTTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGttt                                                                           >  1:24570/1‑149 (MQ=255)                                                                           gAGCAGC‑GCCGGCGCCAAGGCATT‑GAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGC‑GGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCaa                                                                                              <  1:306368/125‑1 (MQ=255)                                                                           gAGCAGC‑GCCGGCGCCAAGGCATT‑GAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGC‑GGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCaa                                                                                              >  2:306368/1‑125 (MQ=255)                                                                            agcagc‑GCCGGCGCCAAGGCATT‑GAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGC‑GGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTtc                                                                         >  5:214029/1‑145 (MQ=255)                                                                            agcagc‑GCCGGCGCCAAGGCATT‑GAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGC‑GGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTtc                                                                         <  6:214029/145‑1 (MQ=255)                                                                              cagc‑GCCGGCGCCAAGGCATT‑GAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGC‑GGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTg                                                                   >  2:358204/1‑149 (MQ=255)                                                                                    ccGGCGCCAAGGCATT‑GAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGC‑GGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCGAGACCAAACTGGTGGTGGTTTCTCCTTGGGCaa                                                              >  2:188981/1‑149 (MQ=255)                                                                                          ccAAGGCATT‑GAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGC‑GGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTGGGCAACCAAAt                                                        <  6:426724/149‑1 (MQ=255)                                                                                           cAAGGCATT‑GAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGC‑GGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTGGGCAACCAAATg                                                       >  6:305560/1‑149 (MQ=255)                                                                                           cAAGGCATT‑GAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGC‑GGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGACCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTGGGCAACCAAATg                                                       >  2:82049/1‑149 (MQ=255)                                                                                            aaGGCATT‑GAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGC‑GGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTGGGCAACCAAATGc                                                      <  2:308328/149‑1 (MQ=255)                                                                                            aaGGCATT‑GAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGC‑GGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGACCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTGGGCAACCAAATGc                                                      <  1:82049/149‑1 (MQ=255)                                                                                             aGGCATT‑GAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGC‑GGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTGGGCAACCAAATGc                                                      <  4:354883/148‑1 (MQ=255)                                                                                             aGGCATT‑GAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGC‑GGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTGGGCAACCAAATGc                                                      >  3:354883/1‑148 (MQ=255)                                                                                             aGGCATT‑GAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGC‑GGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTGGGCAACCAAATGCg                                                     <  2:189843/149‑1 (MQ=255)                                                                                              ggCATT‑GAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGC‑GGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTGGGCAACCAAATGCgg                                                    <  1:87302/149‑1 (MQ=255)                                                                                              ggCATT‑GAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGC‑GGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTGGGCAACCAAATGCgg                                                    <  6:125360/149‑1 (MQ=255)                                                                                                  tt‑GAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGC‑GGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTGGGCAACCAAATGCGGCAtt                                                >  2:396562/1‑149 (MQ=255)                                                                                                     gAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGC‑GGTTTTAGCCAGTGCCAACTGGGTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTGGGCAACCAAATGCGGCATTg                                               <  2:355597/148‑1 (MQ=255)                                                                                                     gAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGC‑GGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTGGGCAACCAAATGCGGCATTg                                               <  1:218961/148‑1 (MQ=255)                                                                                                     gAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGC‑GGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTGGGCAACCAAATGCGGCATTg                                               <  7:352347/148‑1 (MQ=255)                                                                                                     gAAATTGTTAGAAGCATCTACAATTGGGAGCAGAGTGCTCAGGC‑GGTTTTAGACACTGCCAACTGGCTCAAACACTAAATGCCGGCCGCCACCGCCACCAACACCAAACAGTTGTTGGTATCTACTTTGAAAAACAAATGCGGCATTg                                               >  4:238078/1‑148 (MQ=255)                                                                                                       aCTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGC‑GGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTGGGCAACCAAATGCGGCATTGCt                                             <  2:122639/148‑1 (MQ=255)                                                                                                         ttGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGC‑GGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTGGGCAACCAAATGCGGCATTGCTGa                                           >  1:254849/1‑148 (MQ=255)                                                                                                              cGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGC‑GGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTGGGCAACCAAATGCGGCATTGCTGAATaca                                      <  2:24570/148‑1 (MQ=255)                                                                                                               gcagcaTCTACACTTGGGAGCAGAGTGCGCAGGC‑GGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTGGGCAACCAAATGCGGCATTGCTGAATacac                                     <  5:305560/148‑1 (MQ=255)                                                                                                                cagcaTCTACACTTGGGAGCAGAGTGCGCAGGC‑GGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTGGGCAACCAAATGCGGCATTGCTGAATACACTc                                   <  2:341413/149‑1 (MQ=255)                                                                                                                    aTCTACACTTGGGAGCAGAGTGCGGAGGC‑GGTTTTAGTCAGTGCGTAGTGGCTCTATCTGTAAATGCCGGCGTCCGCCGCCAGCAAGAGCAAAGTGTTGGTGGTTtc                                                                         <  1:302233/107‑1 (MQ=255)                                                                                                                    aTCTACACTTGGGAGCAGAGTGCGCAGGC‑GGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTGGGCAACCAAATGCGGCATTGCTGAATACACTCAgg                                >  4:325823/1‑148 (MQ=255)                                                                                                                      cTACACTTGGGAGCAGAGTGCGCAGGC‑GGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTGGGCAACCAAATGCGGCATTGCTGAATACACTCAGGCAt                             >  8:44643/1‑149 (MQ=255)                                                                                                                       tACACTTGGGAGCAGAGTGCGCAGGC‑GGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTGGGCAACCAAATGCGGCATTGCTGAATACACTCAGGCAtt                            >  4:251566/1‑149 (MQ=255)                                                                                                                        acacTTGGGAGCAGAGTGCGCAGGC‑GGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTGGGCAACCAAATGCGGCATTGCTGAATACACTCAGGCATTg                           >  8:335632/1‑149 (MQ=255)                                                                                                                        acacTTGGGAGCAGAGTGCGCAGGC‑GGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTGGGCAACCAAATGCGGCATTGCTGAATACACTCAGGCATTg                           >  5:53669/1‑149 (MQ=255)                                                                                                                         cacTTGGGAGCAGAGTGCGCAGGC‑GGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTGGGCAACCAAATGCGGCATTGCTGAATACACTCAGGCa                              <  3:325823/145‑1 (MQ=255)                                                                                                                             tGGGAGCAGAGTGCGCAGGC‑GGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTGGGCAACCAAATGCGGCATTGCTGAATACACTCAg                                 >  7:279953/1‑138 (MQ=255)                                                                                                                             tGGGAGCAGAGTGCGCAGGC‑GGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTGGGCAACCAAATGCGGCATTGCTGAATACACTCAg                                 <  8:279953/138‑1 (MQ=255)                                                                                                                               ggAGCAGAGTGCGCAGGC‑GGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTGGGCAACCAAATGCGGCATTGCTGAATACACTCAGGCATTGCTCa                       >  1:302790/1‑146 (MQ=255)                                                                                                                               ggAGCAGAGTGCGCAGGC‑GGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTGGGCAACCAAATGCGGCATTGCTGAATACACTCAGGCATTGCTCa                       <  2:302790/146‑1 (MQ=255)                                                                                                                               ggAGCAGAGTGCGCAGGC‑GGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTGGGCAACCAAATGCGGCATTGCTGAATACACTCAGGCATTGCTCAGt                     <  3:331521/148‑1 (MQ=255)                                                                                                                               ggAGCAGAGTGCGCAGGC‑GGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTGGGCAACCAAATGCGGCATTGCTGAATACACTCAGGCATTGCTCAGt                     <  5:210371/148‑1 (MQ=255)                                                                                                                                       gTGCGCAGGC‑GGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTGGGCAACCAAATGCg                                                     <  3:400556/108‑1 (MQ=255)                                                                                                                                       gTGCGCAGGC‑GGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTGGGCAACCAAATGCg                                                     >  4:400556/1‑108 (MQ=255)                                                                                                                                          cgcAGGC‑GGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTGGGCAACCAAATGCGGCATTGCTGAATACACTCAGGCATTGCTCAGTGAGTTCAAGCAg         <  3:238078/149‑1 (MQ=255)                                                                                                                                               gC‑GGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTGGGCAACCAAATGCGGCATTGCTGAATACACTCAg                                 <  4:197917/120‑1 (MQ=255)                                                                                                                                               gC‑GGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTGGGCAACCAAATGCGGCATTGCTGAATACACTCAGGCATTGCTCAGTGa                   <  7:131043/134‑1 (MQ=255)                                                                                                                                               gC‑GGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTGGGCAACCAAATGCGGCATTGCTGAATACACTCAGGCATTGCTCAGTGa                   >  8:131043/1‑134 (MQ=255)                                                                                                                                               gC‑GGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTGGGCAACCAAATGCGGCATTGCTGAATACACTCAGGCATTGCTCAGTGAGTTCAAGCAGGACGa    >  4:191715/1‑149 (MQ=255)                                                                                                                                                  tgTTTGAGCCAGTTGGCACTGGCTAAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTGGGCAACCAAATGCGGCATTGCTGAATACACTCAGGCATTGCTCAGTGa                   >  7:332781/2‑132 (MQ=255)                                                                                                                                                  tgTTTGAGCCAGTTGGCACTGGCTAAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTGGGCAACCAAATGCGGCATTGCTGAATACACTCAGGCATTGCTCAGTGa                   <  8:332781/131‑1 (MQ=255)                                                                                                                                                  tgTTTGAGCCAGTTGGCACTGGCTAAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTGGGCAACCAAATGCGGCATTGCTGAATACACTCAGGCATTGCTCAGTGAGTTCAAGCAGGACGAt   >  2:200997/2‑148 (MQ=255)                                                                                                                                                  ggTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTGGGCAACCAAATGCGGCATTGCTGAATACACTCAGGCATTGCTCAGTGAGTTCAAGCAGGACGAtt  <  1:358204/149‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                AACCGGACAAGCATACCTTGGCCCTCAAGATGGGCGAACTGCGCGGCAAAATCGTCGAAGGCGACGATGCACTCGAGCAGC‑GCCGGCGCCAAGGCATT‑GAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGC‑GGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTGGGCAACCAAATGCGGCATTGCTGAATACACTCAGGCATTGCTCAGTGAGTTCAAGCAGGACGATT  >  NC_002947/1991008‑1991298

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