Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,991,164	C→G	7.2%	A805G (GCC→GGC) ‡	PP_1780 →	mannosyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,991,164	0	C	G	7.2%	259.7 / 6.7	81	A805G (GCC→GGC) ‡	PP_1780	mannosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (36/39);  new base G (2/4);  total (38/43)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

ATACCTTGGCCCTCAAGATGGGCGAACTGCGCGGCAAAATCGTCGAAGGCGACGATGCACTCGAGCAGC‑GCCGGCGCCAAGGCATT‑GAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGC‑GGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTGGGCAACCAAATGCGGCATTGCTGAATACACTCAGGCATTGCTCAGTGAGTTCAAGCAGGACGATTACGACGTT  >  NC_002947/1991020‑1991306                                                                                                                                                    |                                                                                                                                               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                                                                                                                         <  6:20400/148‑1 (MQ=255)                     ggCGAACTGCGCGGCAAAATCGTCGAAGGCGACGATGCACTCGAGCAGC‑GCCGGCGCCAAGGCATT‑GAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGC‑GGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGc                                                                                                                         <  8:314200/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                                                                                                                                      <  4:337782/105‑1 (MQ=255)                                                   gacgaTGCACTCGAGCAGC‑GCCGGCGCCAAGGCATT‑GAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGC‑GGTTTTAGCCAGTGCCAACTGGCTc                                                                                                                                      >  3:337782/1‑105 (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                                                              >  3:354883/1‑148 (MQ=255)                                                                                 aGGCATT‑GAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGC‑GGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTGGGCAACCAAATGc                                                              <  4:354883/148‑1 (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                                                       <  7:352347/148‑1 (MQ=255)                                                                                         gAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGC‑GGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTGGGCAACCAAATGCGGCATTg                                                       <  1:218961/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                                   >  5:53669/1‑149 (MQ=255)                                                                                                            acacTTGGGAGCAGAGTGCGCAGGC‑GGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTGGGCAACCAAATGCGGCATTGCTGAATACACTCAGGCATTg                                   >  8:335632/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                               <  2:302790/146‑1 (MQ=255)                                                                                                                   ggAGCAGAGTGCGCAGGC‑GGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTGGGCAACCAAATGCGGCATTGCTGAATACACTCAGGCATTGCTCa                               >  1:302790/1‑146 (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)                                                                                                                                              ccAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTGGGCAACCAAATGCGGCATTGCTGAATACACTCAGGCATTGCTCAGTGAGTTCAAGCAGGACGATTACGACGtt  >  5:62728/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                               ATACCTTGGCCCTCAAGATGGGCGAACTGCGCGGCAAAATCGTCGAAGGCGACGATGCACTCGAGCAGC‑GCCGGCGCCAAGGCATT‑GAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGC‑GGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTGGGCAACCAAATGCGGCATTGCTGAATACACTCAGGCATTGCTCAGTGAGTTCAAGCAGGACGATTACGACGTT  >  NC_002947/1991020‑1991306

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