Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,066,882	G→A	7.2%	intergenic (+92/‑80)	PP_2678 → / → qedH‑II	hydrolase/quinoprotein ethanol dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,066,882	0	G	A	7.2%	158.0 / 2.2	56	intergenic (+92/‑80)	PP_2678/qedH‑II	hydrolase/quinoprotein ethanol dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (27/24);  new base A (2/2);  total (29/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.51e-01

AGTGTCAGCCACCTGTACCCGCGTTACGAGCGCCAGCAGTTGCAGCGCTTGCCATGACATTTTTTGTGCCTGTCAGCCCCCGTGCCGCCCACCCATCACTCAGGCCGCTACCAAGGAACTAGAAAACTCCGCACTGCGAGCAATTGTTGGCA‑ATCGGG‑‑‑TCGGGACCAAGAATTCCCGACACACCGGGAAATTTTCCCGGGCACAACAAAAAAACCAAAGGTAGCCGTCATGACCCGATCCCCACGTCGCCCCTTGTTCGCCGTGAGCCTGGTGCTCAGCGCCATGCTGCTTGCCG  >  NC_002947/3066734‑3067028                                                                                                                                                     |                                                                                                                                                       aGTGTCAGCCACCTGTACCCGCGTTACGAGCGCCAGCAGTTGCAGCGCTTGCCATGACATTTTTTGTGCCTGTCAGCCCCCGTGCCGCCCACCCATCACTCAGGCCGCTACCAAGGAACTAGAAAACTCCGCACTGCGAGCAAttgttg                                                                                                                                                        <  7:22125/149‑1 (MQ=255) aGTGTCAGCCACCTGTACCCGCGTTACGAGCGCCAGCAGTTGCAGCGCTTGCCATGACATTTTTTGTGCCTGTCAGCCCCCGTGCCGCCCACCCATCACTCAGGCCGCTACCAAGGAACTAGAAAACTCCGCACTGCGAGCAAttgttg                                                                                                                                                        <  5:188834/149‑1 (MQ=255)  gtgtCAGCCACCTGTACCCGCGTTACGAGCGCCAGCAGTTGCAGCGCTTGCCATGACATTTTTTGTGCCTGTCAGCCCCCGTGCCGCCCACCCATCACTCAGGCCGCTACCAAGGAACTAGAAAACTCCGCACTGCGAGCAATTGTTgg                                                                                                                                                       <  3:300060/149‑1 (MQ=255)       aGCCACCTGTACCCGCGTTACGAGCGCCAGCAGTTGCAGCGCTTGCCATGACA‑TTTTTGTGCCTGTCAGCCCCCGTGCCGCCCACCCATCACTCAGGCCGCTACCAAGGAACTAGAAAACTCCGCACTGCGAGCAATTGTTGGCA‑Atcg                                                                                                                                                >  4:210249/1‑149 (MQ=255)         ccACCTGTACCCGCGTTACGAGCGCCAGCAGTTGCAGCGCTTGCCATGACATTTTTTGTGCCTGTCAGCCCCCGTGCCGCCCACCCATCACTCAGGCCGCTACCAAGGAACTAGAAAACTCCGCACTGCGAGCAATTGTTGGCA‑Atcgg                                                                                                                                               <  8:92525/149‑1 (MQ=255)           aCCTGTACCCGCGTTACGAGCGCCAGCAGTTGCAGCGCTTGCCATGACATTTTTTGTGCCTGTCAGCCCCCGTGCCGCCCACCCATCACTCAGGCCGCTACCAAGGAACTAGAAAACTCCGCACTGCGAGCAATTGTTGGCA‑Atcggg‑‑‑t                                                                                                                                          >  4:31644/1‑149 (MQ=255)             cTGTACCCGCGTTACGAGCGCCAGCAGTTGCAGCGCTTGCCATGACATTTTTTGTGCCTGTCAGCCCCCGTGCCGCCCACCCATCACTCAGGCCGCTACCAAGGAACTAGAAAACTCCGCACTGCGAGCAATTGTTGGCA‑Atcggg‑‑‑tcg                                                                                                                                        <  2:278758/149‑1 (MQ=255)              tGTACCCGCGTTACGAGCGCCAGCAGTTGCAGCGCTTGCCATGACATTTTTTGTGCCTGTCAGCCCCCGTGCCGCCCACCCATCCCTCAGGCCGCTACCAAGGAACTAGAAAACTCCGCACTGCGAGCAATTGTTGGCA‑Atcggg‑‑‑tcgg                                                                                                                                       <  4:27773/149‑1 (MQ=255)                tACCCGCGTTACGAGCGCCAGCAGTTGCAGCGCTTGCCATGACATTTTTTGTGCCTGTCAGCCCCCGTGCCGCCCACCCATCACTCAGGCCGCTACCAAGGAACTAGAAAACTCCGCACTGCGAGCAATTGTTGGCA‑ATCGGG‑‑‑TCGGGa                                                                                                                                     <  4:27435/149‑1 (MQ=255)                tACCCGCGTTACGAGCGCCAGCAGTTGCAGCGCTTGCCATGACATTTTTTGTGCCTGTCAGCCCCCGTGCCGCCCACCCATCACTCAGGCCGCTACCAAGGAACTAGAAAACTCCGCACTGCGAGCAATTGTTGGCA‑ATCGGG‑‑‑TCGGGa                                                                                                                                     >  2:161287/1‑149 (MQ=255)                tACCCGCGTTACGAGCGCCAGCAGTTGCAGCGCTTGCCATGACATTTTTTGTGCCTGTCAGCCCCCGTGCCGCCCACCCATCACTCAGGCCGCTACCAAGGAACTAGAAAACTCCGCACTGCGAGCAATTGTTGGCA‑ATCGGG‑‑‑TCGGGa                                                                                                                                     <  3:31644/149‑1 (MQ=255)                 aCCCGCGTTACGAGCGCCAGCAGTTGCAGCGCTTGCCATGACATTTTTTGTGCCTGTCAGCCCCCGTGCCGCCCACCCATCACTCAGGCCGCTACCAAGGAACTAGAAAACTCCGCACTGCGAGCAATTGTTGGCA‑ATCGGG‑‑‑TCGGGAc                                                                                                                                    >  5:71705/1‑149 (MQ=255)                     gcgTTACGAGCGCCAGCAGTTGCAGCGCTTGCCATGACATTTTTTGTGCCTGTCAGCCCCCGTGCCGCCCACCCATCACTCAGGCCGCTACCAAGGAACTAGAAAACTCCGCACTGCGAGCAATTGTTGGCA‑ATCGGG‑‑‑TCGGGAcc                                                                                                                                   >  3:324185/1‑146 (MQ=255)                      cgTTACGAGCGCCAGCAGTTGCAGCGCTTGCCATGACATTTTTTGTGCCTGTCAGCCCCCGTGCCGCCCACCCATCACTCAGGCCGCTACCAAGGAACTAGAAAACTCCGCACTGCGAGCAATTGTTGGCA‑ATCGGG‑‑‑TCGGGACCAAGa                                                                                                                               <  1:161287/149‑1 (MQ=255)                       gTTACGAGCGCCAGCAGTTGCAGCGCTTGCCATGACATTTTTTGTGCCTGTCAGCCCCCGTGCCGCCCACCCATCACTCAGGCCGCTACCAAGGAACTAGAAAACTCCGCACTGCGAGCAATTGTTGGCA‑ATCGGG‑‑‑TCGGGACCAAGaa                                                                                                                              >  3:363708/1‑149 (MQ=255)                             aGCGCCAGCAGTTGCAGCGCTTGCCATGACATTTTTTGTGCCTGTCAGCCCCCGTGCCGCCCACCCATCACTCAGGCCGCTACCAAGGAACTAGAAAACTCCGCACTGCGAGCAATTGTTGGCA‑ATCGGG‑‑‑TCGGGACCAAGAATTccc                                                                                                                         >  6:190835/1‑148 (MQ=255)                              gcgcCAGCAGTTGCAGCGCTTGCCATGACATTTTTTGTGCCTGTCAGCCCCCGTGCCGCCCACCCATCACTCAGGCCGCTACCAAGGAACTAGAAAACTCCGCACTGCGAGCAATTGTTGGCA‑ATCGGG‑‑‑TCGGGACCAAGAATTCCCGa                                                                                                                       <  5:383425/149‑1 (MQ=255)                                   agcagTTGCAGCGCTTGCCATGACATTTTTTGTGCCTGTCAGCCCCCGTGCCGCCCACCCATCACTCAGGCCGCTACCAAGGAACTAGAAAACTCCGCACTGCGAGCAATTGTTGGCA‑ATCGGG‑‑‑TCGGGACCAAGAATTCCCGacacac                                                                                                                  >  7:7450/1‑149 (MQ=255)                                             gcgaTTGCCATGACATTTTTTGTGCCTGTCAGCCCCCGTGCCGCCCACCCATCACTCAGGCCGCTACCAAGGAACTAGAAAACTCCGCACTGCGAGCAATTGTTGGCA‑ATCGGG‑‑‑TCGGGACCAAGAATTCCCGACACACCGGGAAAttt                                                                                                        >  7:329323/1‑149 (MQ=255)                                                cTTGCCATGACATTTTTTGTGCCTGTCAGCCCCCGTGCCGCCCACCCATCACTCAGGCCGCTACCAAGGAACTAGAAAACTCCGCACTGCGAGCAATTGTTGGCA‑ATCGGG‑‑‑TCGGGACCAAGAATTCCCGACACCCCGGGAAATTTTcc                                                                                                     >  6:110318/1‑149 (MQ=255)                                                cTTGCCATGACATTTTTTGTGCCTGTCAGCCCCCGTGCCGCCCACCCATCACTCAGGCCGCTACCAAGGAACTAGAAAACTCCGCACTGCGAGCAATTGTTGGCA‑ATCGGG‑‑‑TCGGGACCAAGAATTCCCGACACACCGGGAAATTTTcc                                                                                                     <  5:190835/149‑1 (MQ=255)                                                 ttGCCATGACATTTTTTGTGCCTGTCAGCCCCCGTGCCGCCCACCCATCACTCAGGCCGCTACCAAGGAACTAGAAAACTCCGCACTGCGAGCAATTGTTAGCA‑ATCGGG‑‑‑TCGGGACCAAGAATTCCCGACACACCGGGAAATTTTccc                                                                                                    >  2:354537/1‑149 (MQ=255)                                                   gCCATGACATTTTTTGTGCCTGTCAGCCCCCGTGCCGCCCACCCATCACTCAGGCCGCTACCAAGGAACTAGAAAACTCCGCACTGCGAGCAATTGTTGGCA‑ATCGGG‑‑‑TCGGGACCAAGAATTCCCGACACACCGGGAAATTTTccc                                                                                                    <  4:363708/147‑1 (MQ=255)                                                    ccATGACATTTTTTGTGCCTGTCAGCCCCCGTGCCGCCCACCCATCACTCAGGCCGCTACCAAGGAACTAGAAAACTCCGCACTGCGAGCAATTGTTGGCA‑ATCGGT‑‑‑TCGGGACCAAGAATTCCCGa                                                                                                                       >  2:64628/1‑127 (MQ=255)                                                    ccATGACATTTTTTGTGCCTGTCAGCCCCCGTGCCGCCCACCCATCACTCAGGCCGCTACCAAGGAACTAGAAAACTCCGCACTGCGAGCAATTGTTGGCA‑ATCGGA‑‑‑TCGGGACCAAGAATTCCCGa                                                                                                                       <  1:64628/127‑1 (MQ=255)                                                      aTGACATTTTTTGTGCCTGTCAGCCCCCGTGCCGCCCACCCATCACTCAGGCCGCTACCAAGGAACTAGAAAACTCCGCACTGCGAGCAATTGTTGGCA‑ATCGGG‑‑‑TCGGGACCAAGAATTCCCGACACACCGGGAAATTTTCCCGGGca                                                                                               >  3:204355/1‑149 (MQ=255)                                                       tGACATTTTTTGTGCCTGTCAGCCCCCGTGCCGCCCACCCATCACTCAGGCCGCTACCAAGGAACTAGAAAACTCCGCACTGCGAGCAATTGTTGGCA‑ATCGGG‑‑‑TCGGGACCAAGAATTCCCGACACACCGGGAAATTTTCCCGGGcac                                                                                              <  5:344371/149‑1 (MQ=255)                                                           aTTTTTTGTGCCTGTCAGCCCCCGTGCCGCCCACCCATCACTCAGGCCGCTACCAAGGAACTAGAAAACTCCGCACTGCGAGCAATTGTTAGCA‑ATCGGG‑‑‑TCGGGACCAAGAATTCCCGACACACCGGGAAATTTTCCCGGGCacaac                                                                                           >  5:255739/1‑148 (MQ=255)                                                                ttGTGCCTGTCAGCCCCCGTGCCGCCCACCCATCACTCAGGCCGCTACCAAGGAACTAGAAAACTCCGCACTGCGAGCAATTGTTGGCA‑ATCGGG‑‑‑TCGGGACCAAGAATTCCCGACACACCGGGAAATTTTCCCGGGCACAACaaaaa                                                                                      <  8:7450/148‑1 (MQ=255)                                                                         tCAGCCCCCGTGCCGCCCACCCATCACTCAGGCCGCTACCAAGGAACTAGAAAACTCCGCACTGCGAGCAATTGTTAGCA‑ATCGGG‑‑‑TCGGGACCAAGAATTCCCGACACACCGGGAAATTTTCCCGGGCACAACAAAAAAACCAAAGGt                                                                            <  6:255739/149‑1 (MQ=255)                                                                             cccccGTGCCGCCCACCCATCACTCAGGCCGCTACCAAGGAACTAGAAAACTCCGCACTGCGAGCAATTGTTGGCA‑ATCGGG‑‑‑TCGGGACCAAGAATTCCCGACACACCGGGAAATTTTCCCGGGCACAACaaaaaaa                                                                                    <  3:174984/137‑1 (MQ=255)                                                                             cccccGTGCCGCCCACCCATCACTCAGGCCGCTACCAAGGAACTAGAAAACTCCGCACTGCGAGCAATTGTTGGCA‑ATCGGG‑‑‑TCGGGACCAAGAATTCCCGACACACCGGGAAATTTTCCCGGGCACAACaaaaaaa                                                                                    >  4:174984/1‑137 (MQ=255)                                                                             cccccGTGCCGCCCACCCATCACTCAGGCCGCTACCAAGGAACTAGAAAACTCCGCACTGCGAGCAATTGTTGGCA‑ATCGGG‑‑‑TCGGGACCAAGAATTCCCGACACACCGGGAAATTTTCCCGGGCACAACAAAAAAACCAAAGGTAGcc                                                                        >  1:149257/1‑149 (MQ=255)                                                                              ccccGTGCCGCCCACCCATCACTCAGGCCGCTACCAAGGAACTAGAAAACTCCGCACTGCGAGCAATTGTTGGCA‑ATCGGG‑‑‑TCGGGACCAAGAATTCCCGACACACCGGGAAATTTTCCCGGGCacaa                                                                                            >  8:26061/1‑128 (MQ=255)                                                                              ccccGTGCCGCCCACCCATCACTCAGGCCGCTACCAAGGAACTAGAAAACTCCGCACTGCGAGCAATTGTTGGCA‑ATCGGG‑‑‑TCGGGACCAAGAATTCCCGACACACCGGGAAATTTTCCCGGGCacaa                                                                                            <  7:26061/128‑1 (MQ=255)                                                                                 cGTGCCGCCCACCCATCACTCAGGCCGCTACCAAGGAACTAGAAAACTCCGCACTGCGAGCAATTGTTGGCA‑ATCGGG‑‑‑TCGGGACCAAGAATTCCCGACACACCGGGAAATTTTCCCGGGCACAACAAAAAAACCAAAGGTAGCCGTCa                                                                    <  3:210249/149‑1 (MQ=255)                                                                                      cgccCACCCATCACTCAGGCCGCTACCAAGGAACTAGAAAACTCCGCACTGCGAGCAATTGTTGGCA‑ATCGGG‑‑‑TCGGGACCAAGAATTCCCGACACACCGGGAAATTTTCCCGGGCACAACAAAAAAACCAAAGGTAGCCGTCATGAcc                                                               >  7:157460/1‑149 (MQ=255)                                                                                       gccCACCCATCACTCAGGCCGCTACCAAGGAACTAGAAAACTCCGCACTGCGAGCAATTGTTGGCA‑ATCGGG‑‑‑TCGGGACCAAGAATTCCCGACACACCGGGAAATTTTCCCGGGca                                                                                               <  1:44183/116‑1 (MQ=255)                                                                                       gccCACCCATCACTCAGGCCGCTACCAAGGAACTAGAAAACTCCGCACTGCGAGCAATTGTTGGCA‑ATCGGG‑‑‑TCGGGACCAAGAATTCCCGACACACCGGGAAATTTTCCCGGGca                                                                                               >  2:44183/1‑116 (MQ=255)                                                                                           acccaTCACTCAGGCCGCTACCAAGGAACTAGAAAACTCCGCACTGCGAGCAATTGTTGGCA‑ATCGGGAGAAGGGGACCAAGAATTCCCGACACACCGGGAAATTTTCCCGGGCACAACAAAAAAACCAAAGGTAGCCGTCATGAccc                                                              >  1:307542/1‑148 (MQ=255)                                                                                           acccaTCACTCAGGCCGCTACCAAGGAACTAGAAAACTCCGCACTGCGAGCAATTGTTGGCA‑ATCGGG‑‑‑TCGGGACCAAGAATTCCCGACACACCGGGAAATTTTCCCGGGCACAACaaaaa                                                                                      <  5:367453/121‑1 (MQ=255)                                                                                           acccaTCACTCAGGCCGCTACCAAGGAACTAGAAAACTCCGCACTGCGAGCAATTGTTGGCA‑ATCGGG‑‑‑TCGGGACCAAGAATTCCCGACACACCGGGAAATTTTCCCGGGCACAACaaaaa                                                                                      >  6:367453/1‑121 (MQ=255)                                                                                           acccaTCACTCAGGCCGCTACCAAGGAACTAGAAAACTCCGCACTGCGAGCAATTGTTGGCA‑ATCGGG‑‑‑TCGGGACCAAGAATTCCCGACACACCGGGAAATTTTCCCGGGCACAACAAAAAAACCAAAGGTAGCCGTCATGACCCGATc                                                          >  8:113677/1‑149 (MQ=255)                                                                                            cccaTCACTCAGGCCGCTACCAAGGAACTAGAAAACTCCGCACTGCGAGCAATTGTTAGCA‑ATCGGG‑‑‑TCGGGACCAAGAATTCCCGACACACCGGGAAATTTTCCCGGGCACAACAAAAAAACCAAAGGTAGCCGTCATGACCCGATcc                                                         <  1:354537/149‑1 (MQ=255)                                                                                                taccTCAGGACGCTAGCAAGGAACTCGAAAACTCCGCCCTGCGAGAAATTGTTGGCA‑ATCGGG‑‑‑TCGGGACCAAGAATTCCCGACACACCGGGAAATTTTCCCGGGCACAACAAAAAAACCAAAGGTAGCCGTCATGACCCGATCCCCAc                                                     <  8:329323/146‑1 (MQ=255)                                                                                                        gCCGCTACCAAGGAACTAGAAAACTCCGCACTGCGAGCAATTGTTGGCA‑ATCGGG‑‑‑TCGGGACCAAGAATTCCCGACACACCGGGAAATTTTCCCGGGCACAACAAAAAAACCAAAGGTAGCCGTCATGACCCGATCCCCACGTCGcccc                                             >  2:80957/1‑149 (MQ=255)                                                                                                        gCCGCTACCAAGGAACTAGAAAACTCCGCACTGCGAGCAATTGTTGGCA‑ATCGGG‑‑‑TCGGGACCAAGAATTCCCGACACACCGGGAAATTTTCCCGGGCACAACAAAAAAACCAAAGGTAGCCGTCATGACCCGATCCCCACGTCGcccc                                             >  6:261622/1‑149 (MQ=255)                                                                                                          cccaaccaaccaACAACAAAACACCCCACACCGAGCAATTGTTTGCA‑ATCGGG‑‑‑TCTGGACCAAGAATTCCCGACACACCGGGAAATTTTCCCGGGCACAACAAAAAAACCAAAGGTAGCCGTCATGACCCGATCCCCACGTCGCCCCtt                                           <  2:149257/145‑1 (MQ=255)                                                                                                                  aGGAACTAGAAAACTCCGCACTGCGAGCAATTGTTGGCA‑ATCGGG‑‑‑TCGGGACCAAGAATTCCCGACACACCGGGAAATTTTCCCGGGCACAACAAAAAAACCAAAGGTAGCCGTCATGACCCGATCCCCACGTCGCCCCTTGTTCGCCg                                   <  5:261622/149‑1 (MQ=255)                                                                                                                  aGGAACTAGAAAACTCCGCACTGCGAGCAATTGTTGGCA‑ATCGGG‑‑‑TCGGGACCAAGAATTCCCGACACACCGGGAAATTTTCCCGGGCACAACAAAAAAACCAAAGGTAGCCGTCATGACCCGATCCCCACGTCGCCCCTTGTTCGCCg                                   <  1:80957/149‑1 (MQ=255)                                                                                                                   ggAACTAGAAAACTCCGCACTGCGAGCAATTGTTGGCA‑ATCGGG‑‑‑TCGGGACCAAGAATTCCCGACACACCGGGAAATTTTCCCGGGCACAACAAAAAAACCAAAGGTAGCCGTCATGACCCGATCCCCACGTCGCCCCTTGTTCGCCGt                                  <  5:110318/149‑1 (MQ=255)                                                                                                                   ggAACTAGAAAACTCCGCACTGCGAGCAATTGTTGGCA‑ATCGGG‑‑‑TCGGAACCAAGAATTCCCGACACACCGGGAAATTTTCCCGGGCACAACAAAAAAACCAAAGGTAGCCGTCATGACCCGATCCCCACGTCGCCCCTTGTTCGCCGt                                  >  1:91071/1‑149 (MQ=255)                                                                                                                   gagcttAGAAAACTCCGCACTGCGAGCAATTGTTGGCA‑ATCGGG‑‑‑TCGGGACCAAGAATTCCCGACACACCGGGAAATTTTCCCGGGCACAACAAAAAAACCAAAGGTAGCCGTCATGACCCGATCCCCACGTCGCCCCTTGTTCGCCGt                                  <  7:237146/144‑1 (MQ=255)                                                                                                                         aGAAAACTCCGCACTGCGAGCAATTGTTGGCA‑ATCGGG‑‑‑TCGGGACCAAGAATTCCCGACACACCGGGAAATTTTCCCGGGCACAACAAAAAAACCAAAGGTAGCCGTCATGACCCGATCCCCACGTCGCCCCTTGTTCGCCGt                                  >  8:237146/1‑143 (MQ=255)                                                                                                                          gAAAACTCCGCACTGCGAGCAATTGTTGGCA‑ATCGGG‑‑‑TCGGGACCAAGAATTCCCGACACACCGGGAAATTTTCCCGGGCACAACAAAAAAACCAAAGGTAGCCGTCATGACCCGATCCCCACGTCGCCCCTTGTTCGCCGTGAGCCTg                           >  6:15828/1‑149 (MQ=255)                                                                                                                             aaCTCCGCACTGCGAGCAATTGTTGGCA‑ATCGGG‑‑‑TCGGGACCAAGAATTCCCGACACACCGGGAAATTTTCCCGGGCACAACAAAAAAACCAAAGGTAGCCGTCATGACCCGATCCCCACGTCGCCCCTTGTTCGCCGTGAGCCTGGTg                        <  4:204355/149‑1 (MQ=255)                                                                                                                               cTCCGCACTGCGAGCAATTGTTGGCAGATCGGG‑‑‑TCGGGACCAAGAATTCCCGACACACCGGGAAATTTTCCCGGGCACAACAAAAAAACCAAAGGTAGCCGTCATGACCCGATCCCCACGTCGCCCCTTGTTCGCCGTGAGCCTGGTGc                       >  1:5191/1‑149 (MQ=255)                                                                                                                               cTCCGCACTGCGAGCAATTGTTGGCA‑ATCGGG‑‑‑TCGGGACCAAGAATTCCCGACACACCGGGAAATTTTCCCGGGCACAACAAAAAAACCAAAGGTAGCCGTCATGACCCGATCCCCACGTCGCCCCTTGTTCGCCGTGAGCCTGGTGCt                      <  5:15828/149‑1 (MQ=255)                                                                                                                                  cGCACTGCGAGCAATTGTTGGCA‑ATCGGG‑‑‑TCGGGACCAAGAATTCCCGACACACCGGGAAATTTTCCCGGGCACAACAAAAAAACCAAAGGTAGCCGTCATGACCCGATCCCCACGTCGCCCCTTGTTCGCCGTGAGCCTGGTGCTCAg                   >  3:311212/1‑149 (MQ=255)                                                                                                                                                  gttgGCA‑ATCGGG‑‑‑TCGGGACCAAGAATTCCCGACACACCGGGAAATTTTCCCGGGCACAACAAAAAAACCAAAGGTAGCCGTCATGACCCGATCCCCACGTCGCCCCTTGTTCGCCGTGAGCCTGGTGCTCAGCGCCATGCTGCTTGc    <  4:324185/148‑1 (MQ=255)                                                                                                                                                   ttgGCA‑ATCGGG‑‑‑TCGGGACCAAGAATTCCCGACACACCGGGAAATTTTCCCGGGCACAACAAAAAAACCAAAGGTAGCCGTCATGACCCGATCCCCACGTCGCCCCTTGTTCGCCGTGAGCCTGGTGCTCAGCGACATGCTGAtttcag  >  8:259216/1‑145 (MQ=255)                                                                                                                                                     |                                                                                                                                                       AGTGTCAGCCACCTGTACCCGCGTTACGAGCGCCAGCAGTTGCAGCGCTTGCCATGACATTTTTTGTGCCTGTCAGCCCCCGTGCCGCCCACCCATCACTCAGGCCGCTACCAAGGAACTAGAAAACTCCGCACTGCGAGCAATTGTTGGCA‑ATCGGG‑‑‑TCGGGACCAAGAATTCCCGACACACCGGGAAATTTTCCCGGGCACAACAAAAAAACCAAAGGTAGCCGTCATGACCCGATCCCCACGTCGCCCCTTGTTCGCCGTGAGCCTGGTGCTCAGCGCCATGCTGCTTGCCG  >  NC_002947/3066734‑3067028

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