Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,991,174	C→A	6.0%	L808L (CTC→CTA)	PP_1780 →	mannosyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,991,174	0	C	A	6.0%	270.2 / 6.3	83	L808L (CTC→CTA)	PP_1780	mannosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (39/39);  new base A (2/3);  total (41/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.52e-01

CCCTCAAGATGGGCGAACTGCGCGGCAAAATCGTCGAAGGCGACGATGCACTCGAGCAGC‑GCCGGCGCCAAGGCATT‑GAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGC‑GGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTGGGCAACCAAATGCGGCATTGCTGAATACACTCAGGCATTGCTCAGTGAGTTCAAGCAGGACGATTACGACGTTCGCATCTACTGCGAT  >  NC_002947/1991029‑1991321                                                                                                                                                     |                                                                                                                                                    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                                                  >  8:335632/1‑149 (MQ=255)                                                                                                   acacTTGGGAGCAGAGTGCGCAGGC‑GGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTGGGCAACCAAATGCGGCATTGCTGAATACACTCAGGCATTg                                                  >  5:53669/1‑149 (MQ=255)                                                                                                    cacTTGGGAGCAGAGTGCGCAGGC‑GGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTGGGCAACCAAATGCGGCATTGCTGAATACACTCAGGCa                                                     <  3:325823/145‑1 (MQ=255)                                                                                                        tGGGAGCAGAGTGCGCAGGC‑GGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTGGGCAACCAAATGCGGCATTGCTGAATACACTCAg                                                        <  8:279953/138‑1 (MQ=255)                                                                                                        tGGGAGCAGAGTGCGCAGGC‑GGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTGGGCAACCAAATGCGGCATTGCTGAATACACTCAg                                                        >  7:279953/1‑138 (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                                          >  8:131043/1‑134 (MQ=255)                                                                                                                          gC‑GGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTGGGCAACCAAATGCGGCATTGCTGAATACACTCAGGCATTGCTCAGTGa                                          <  7:131043/134‑1 (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)                                                                                                                                                  gCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTGGGCAACCAAATGCGGCATTGCTGAATACACTCAGGCATTGCTCAGTGAGTTCAAGCAGGACGATTACGACGTTCGCATCTACTGc     >  3:161945/1‑148 (MQ=255)                                                                                                                                                  gCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTGGGCAACCAAATGCGGCATTGCTGAATACACTCAGGCATTGCTCAGTGAGTTCAAGCAGGACGATTACGACGTTCGCATCTACTGc     >  4:389990/1‑148 (MQ=255)                                                                                                                                                   cTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTGGGCAACCAAATGCGGCATTGCTGAATACACTCAGGCATTGCTCAGTGAGTTCAAGCAGGACGATTACGACGTTCGCATCTACTGCGa   <  3:191715/149‑1 (MQ=255)                                                                                                                                                    tCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTGGGCAACCAAATGCGGCATTGCTGAATACACTCAGGCATTGCTCAGTGAGTTCAAGCAGGAc                             >  7:287344/1‑122 (MQ=255)                                                                                                                                                    tCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTGGGCAACCAAATGCGGCATTGCTGAATACACTCAGGCATTGCTCAGTGAGTTCAAGCAGGAc                             <  8:287344/122‑1 (MQ=255)                                                                                                                                                    tCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTGGGCAACCAAATGCGGCATTGCTGAATACACTCAGGCATTGCTCAGTGAGTTCAAGCAGGACGATTACGACGTTCGCATCTACTGCGAt  >  8:152463/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                    CCCTCAAGATGGGCGAACTGCGCGGCAAAATCGTCGAAGGCGACGATGCACTCGAGCAGC‑GCCGGCGCCAAGGCATT‑GAACTTGTTCGCAGCATCTACACTTGGGAGCAGAGTGCGCAGGC‑GGTTTTAGCCAGTGCCAACTGGCTCAAACACGAAATGCCGGCGGCCGCCGCCAGCAAGACCAAACTGGTGGTGGTTTCTCCTTGGGCAACCAAATGCGGCATTGCTGAATACACTCAGGCATTGCTCAGTGAGTTCAAGCAGGACGATTACGACGTTCGCATCTACTGCGAT  >  NC_002947/1991029‑1991321

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