Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,562,314	A→T	6.8%	V73E (GTG→GAG)	PP_3145 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,562,314	0	A	T	6.8%	169.1 / 2.3	59	V73E (GTG→GAG)	PP_3145	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (28/27);  new base T (2/2);  total (30/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.01e-01

GCTCGACGAACACCTTGCTGTACTGGTTGATATTGACCTTGGGGTCGATCCAGCGCATCACCGGGTCACCCGACGGGCTCTTGGCTTCTTGCAGCCGGCTGTAGTCTTTGAGGAAGCCGGAGTAATCCTTGGGGTCCA‑‑CGCGGTTGCTGGAGCA‑GGCGGCCACGGCAAGCAGCGCGGCGCACAACAGGGTAATGCGAGGCAATGATTTCATGATGAAAAGGCATCCGTGAGAGGTCGGCCTCAACAACGCTAGCAGGTGCGAGGCAGTTGGCCAGTGGACGGGG  >  NC_002947/3562177‑3562460                                                                                                                                          |                                                                                                                                                      gctCGACGATCACCTTGCTGTACTGGTTGATATTGACCTTGGGGTCGATCCAGCGCCTCACCGGGTCACCCGACGGGCTCTTGGCTTCTTGCAGCCGGCTGTAGTCTTTGAGGAAGCCGGAGTAATCCTTGGGGTCCA‑‑CGCGGTTGCTg                                                                                                                                          <  6:343660/149‑1 (MQ=255) gctCGACGAACACCTTGCTGTACTGGTTGATATTGACCTTGGGGTCGATCCAGCGCATCACCGGGTCACCCGACGGGCTCTTGGCTTCTTGCAGCCGGCTGTAGTCTTTGAGGAAGCCGGAGTAATCCTTGGGGTCCA‑‑CGCGGTTGCTg                                                                                                                                          <  4:7096/149‑1 (MQ=255) gctCGACGAACACCTTGCTGTACTGGTTGATATTGACCTTGGGGTCGATCCAGCGCATCACCGGGTCACCCGACGGGCTCTTGGCTTCTTGCAGCCGGCTGTAGTCTTTGAGGAAGCCGGAGTAATCCTTGGGGTCCA‑‑CGCGGTTGCTg                                                                                                                                          <  7:183030/149‑1 (MQ=255)   tCGACGAACACCTTGCTGTACTGGTTGATATTGACCTTGGGGTCGATCCAGCGCATCACCGGGTCACCCGACGGGCTCTTGGCTTCTTGCAGCCGGCTGTAGTCTTTGAGGAAGCCGGAGTAATCCTTGGGGTCCT‑‑CGCGGTTGCTg                                                                                                                                          >  1:333868/1‑147 (MQ=255)   tCGACGAACACCTTGCTGTACTGGTTGATATTGACCTTGGGGTCGATCCAGCGCATCACCGGGTCACCCGACGGGCTCTTGGCTTCTTGCAGCCGGCTGTAGTCTTTGAGGAAGCCGGAGTAATCCTTGGGGTCCT‑‑CGCGGTTGCTg                                                                                                                                          >  1:338314/1‑147 (MQ=255)        gaACACCTTGCTGTACTGGTTGATATTGACCTTGGGGTCGATCCAGCGCATCACCGGGTCACCCGACGGGCTCTTGGCTTCTTGCAGCCGGCTGTAGTCTTTGAGGAAGCCGGAGTAATCCTTGGGGTCCA‑‑CGCGGTTGCTGGAGCA‑gg                                                                                                                                  >  3:306932/1‑149 (MQ=255)              ctagCTGTACTGGTTGATATTGACCTTGGGGTCGATCCAGCGCATCACCGGGACACACGACGGGCTCTTCGCTTCTTGCAGCCGGCTGTAGTCTTTGAGGAAGCCGGAGTAATCCTTGGGGTCCA‑‑CGCTCTTGCTGCAGAA‑GGCGACCa                                                                                                                            >  4:25289/4‑149 (MQ=255)              cTTGCTGTACTGGTTGATATTGACCTTGGGGTCGATCCAGCGCATCACCTGGTCACCCTACGGGCTCTTGGCTTCTTGCAGCCGGCTGTAGTCTTTGAGGAAGCCGGAGTAATCCTTGGGGTCCA‑‑CGCGGTTGCTGGAGCA‑GGCGGCCa                                                                                                                            <  1:305787/149‑1 (MQ=255)              cTTGCTGTACTGGTTGATATTGACCTTGGGGTCGATCCAGCGCATCACCGGGTCACCCGACGGGCTCTTGGCTTCTTGCAGCCGGCTGTAGTCTTTGAGGAAGCCGGAGTAATCCTTGGGGTCCA‑‑CGCGGTTGCTGGAGCA‑GGCGGCCa                                                                                                                            <  4:306932/149‑1 (MQ=255)              cTTGCTGAACTGGTTGATATTGACCTTGGGGTCGATCAAGCCCATCACCGGGTCACACGACGGGCTCATGGCTTCATGCAGCCCGCAGTAGTCTTTGAGGAAGCCGGAGTAATACTTCCGCTCCA‑‑CGCGGTTCCTGGAGCA‑GGCGGCCa                                                                                                                            >  7:287725/1‑149 (MQ=255)                   tGTACTGGTTGATATTGACCTTGGGGTCGATCCAGCGCATCACCGGGTCACCCGACGGGCTCTTGGCTTCTTGCAGCCGGCTGTAGTCTTTGAGGAAGCCGGAGTAATCCTTGGGGTCCA‑‑CGCGGTTGCTGGAGCA‑GGCGGCCACGGc                                                                                                                        >  8:256543/1‑148 (MQ=255)                   tGTACTGGTTGATATTGACCTTGGGGTCGATCCAGCGCATCACCGGGTCACCCGACGGGCTCTTGGCTTCTTGCAGCCGGCTGTAGTCTTTGAGGAAGCCGGAGTAATCCTTGGGGTCCA‑‑CGCGGTTGCTGGAGCA‑GGCGGCCACGGc                                                                                                                        >  6:184531/1‑148 (MQ=255)                          gTTGATATTGACCTTGGGGTCGATCCAGCGCATCACCGGGTCACCCGACGGGCTCTTGGCTTCTTGCAGCCGGCTGTAGTCATTGAGGAAGCCGGAGTAATCCTTGGGGTCCA‑‑CGCGGTTGCTGGAGCA‑GGCGGCCACGGCAAGCAgcg                                                                                                                >  2:293649/1‑149 (MQ=255)                            tGATATTGACCTTGGGGTCGATCCAGCGCATCACCGGGTCACCCGACGGGCTCTTGGCTTCTTGCAGCCGGCTGTAGTCTTTGAGGAAGCCGGAGTAATCCTTGGGGTCCA‑‑CGCGGTTGCTGGAGCA‑ggcg                                                                                                                                >  6:325707/1‑131 (MQ=255)                            tGATATTGACCTTGGGGTCGATCCAGCGCATCACCGGGTCACCCGACGGGCTCTTGGCTTCTTGCAGCCGGCTGTAGTCTTTGAGGAAGCCGGAGTAATCCTTGGGGTCCA‑‑CGCGGTTGCTGGAGCA‑ggcg                                                                                                                                <  5:325707/131‑1 (MQ=255)                             gATATTGACCTTGGGGTCGATCCAGCGCATCACCGGGTCACCCGACGGGCTCTTGGCTTCTTGCAGCCGGCTGTAGTCTTTGAGGAAGCCGGAGTAATCCTTGGGGTCCA‑‑CGCGGTTGCTGGAGCA‑GGCGGCGACGGCAGGCAGCgcgg                                                                                                             >  4:198117/1‑149 (MQ=255)                              atatTGACCTTGGGGTCGATCCAGCGCATCACCGGGTCACCCGACGGGCTCTTGGCTTCTTGCAGCCGGCTGTAGTCTTTGAGGAAGCCGGAGTAATCCTTGGGGTCCA‑‑CGCGGTTGCTGGAGCA‑GGCGGCCACGGCAAGCAGCgcggc                                                                                                            <  7:256543/149‑1 (MQ=255)                                 ttGACCTTGGGGTCGATCCAGCGCATCACCGGGTCACCCGACGGGCTCTTGGCTTCTTGCAGCCGGCTGTAGTCTTTGAGGAAGCCGGAGTAATCCTTGGGGTCCA‑‑CGCGGTTGCTGGAGCA‑GGCGGCCACGGCAAGCAGCGCGGCGca                                                                                                         >  5:349235/1‑149 (MQ=255)                                 ttGACCTTGGGGTCGATCCAGCGCATCACCGGGTCACCCGACGGGCTCTTGGCTTCTTGCAGCCGGCTGTAGTCTTTGAGGAAGCCGGAGTAATCCTTGGGGTCCA‑‑CGCGGTTGCTGGAGCA‑GGCGGCCACGGCAAGCAGCGCGGCGca                                                                                                         >  2:192742/1‑149 (MQ=255)                                       ttGGGGTCGATCCAGCGCATCACCGGGTCACCCGACGGGCTCTTGGCTTCTTGCAGCCGGCTGTAGTCTTTGAGGAAGCCGGAGTAATCCTTGGGGTCCA‑‑CGCGGTTGCTGGAGCA‑GGCGGCCACGGCAAGCAGCGCGGCGCACAACAg                                                                                                   >  5:102534/1‑149 (MQ=255)                                            gTCGATCCAGCGCATCACCGGGTCACCCGACGGGCTCTTGGCTTCTTGCAGCCGGCTGTAGTCTTTGAGGAAGCCGGAGTAATCCTTGGGGTCCA‑‑CGCGGTTGCTGGAGCA‑GGCGGCCACGGCAAGCAGCGCGGCGCACAACAGGGTaa                                                                                              >  8:207139/1‑149 (MQ=255)                                             tCGATCCAGCGCATCACCGGGTCATCCGACGGGCTCTTGGCTTCTTGCAGCCGGCTGTAGTCTTTGAGGAAGCCGGAGTAATCCTTGGGGTCCA‑‑CGCGGTTGCTGGAGCA‑GGCGGCCACGGCAAGCAGCGCGGCGCACAACAGGGTAAt                                                                                             <  1:349595/149‑1 (MQ=255)                                               gATCCAGCGCATCACCGGGTCACCCGACGGGCTCTTGGCTTCTTGCAGCCGGCTGTAGTCTTTGAGGAAGCCGGAGTAATCCTTGGGGTCCA‑‑CGCGGTTGCTGGAGCA‑GGCGGCCACGGCAAGCAGCGCGGCGCACAAGAGGGTAAtca                                                                                           >  4:68159/1‑147 (MQ=255)                                               gATCCAGCGCATCACCGGGTCACCCGACGGGCTCTTGGCTTCTTGCAGCCGGCTGTAGTCTTTGAGGAAGCCGGAGTAATCCTTGGGGTCCA‑‑CGCGGTTGCTGGAGCA‑GGCGGCCACGGCAAGCAGCGCGGCGCACAACAGGGTAATGc                                                                                           >  8:282724/1‑149 (MQ=255)                                                 tCCAGCGCATCACCGGGTCACCCGACGGGCTCTTGGCTTCTTGCAGCCGGCTGTAGTCTTTGAGGAAGCCGGAGTAATCCTTGGGGTCCA‑‑CGCGGTTGCTGGAGCA‑GGCGGCCACGGCAAGCAGCGCGGCGCACAACAGGGTAATGCGa                                                                                         <  6:102534/149‑1 (MQ=255)                                                   cAGCGCATCACCGGGTCACCCGACGGGCTCTTGGCTTCTTGCAGCCGGCTGTAGTCTTTGAGGAAGCCGGAGTAATCCTTGGGGTCCA‑‑CGCGGTTGCTGGAGCA‑GGCGGCCACGGCAAGCAGCGCGGCGCACAACAGGGTAATGCGAgg                                                                                       <  2:369826/149‑1 (MQ=255)                                                    aGCGCATCACCGGGTCACCCGACGGGCTCTTGGCGTCTTGCAGCCGGCTGTAGTCTTTGAGGCAGCCGGAGTAATCCTTGGGGTCCA‑‑CGCGGTTGCTGGAGAA‑GGCGGCCACGGCAAGCAGCGCGGCGCACAACAGGGTAATGCGAgg                                                                                       >  6:81516/1‑148 (MQ=255)                                                         aTCACCGGGTCACCCGACGGGCTCTTGGCTTCTTGCAGCCGGCTGTAGTCTTTGAGGAAGCCGGAGTAATCCTTGGGGTCCA‑‑CGCGGTTGCTGGAGCA‑GGCGGCCACGGCAAGCAGCGCGGCGCACAACAGGGTAATGCGAGGCAATGa                                                                                 >  4:176018/1‑149 (MQ=255)                                                          tCACCGGGTCACCCGACGGGCTCTTGGCTTCTTGCAGCCGGCTGTAGTCTTTGAGGAAGCCGGAGTAATCCTTGGGGTCCA‑‑CGCGGTTGCTGGAGCA‑GGCGGCCACGGCAAGCAGCGCGGCGCACAACAGGGTAATGCGAGGCAATGAt                                                                                >  3:205961/1‑149 (MQ=255)                                                                ggTCACCCGACGGGCTCTTGGCTTCTTGCAGCCGGCTGTAGTCTTTGAGGAAGCCGGAGTAATCCTTGGGGTCCA‑‑CGCGGTTGCTGGAGCA‑GGCGGCCACGGCAAGCAGCGCGGCGCACAACAGGGTAATGCGAGGCAATGATTTCatg                                                                          <  7:282724/149‑1 (MQ=255)                                                                 gTCACCCGACGGGCTCTTGGCTTCTTGCAGCCGGCTGTAGTCTTTGAGGAAGCCGGAGTAATCCTTGGGGTCCA‑‑CGCGGTTGCTGGAGCATTGCGGCCACGGCAAGCAGCGCGGCGCACAACAGGGTAATGCGAGGCAATGATTTCatg                                                                          <  1:293649/149‑1 (MQ=255)                                                                  tCACCCGACGGGCTCTTGGCTTCTTGCAGCCGGCTGTAGTCTTTGAGGAAGCCGGAGTAATCCTTGGGGTCCA‑‑CGCGGTTGCTGGAGCA‑GGCGGCCACGGCAAGCAGCGCGGCGCACAACAGGGTAATGCGAGGCAATGATTTCatgat                                                                        <  6:349235/149‑1 (MQ=255)                                                                      ccGACGGGCTCTTGGCTTCTTGCAGCCGGCTGTAGTCTTTGAGGAAGCCGGAGTAATCCTTGGGGTCCT‑‑CGCGGTTGCTGGAGCA‑GGCGGCCACGGCAAGCAGCGCGGCGCACAACAGGGTAATGCGAGGCAATGATTTCATGATGaaa                                                                    <  2:333868/149‑1 (MQ=255)                                                                      ccGACGGGCTCTTGGCTTCTTGCAGCCGGCTGTAGTCTTTGAGGAAGCCGGAGTAATCCTTGGGGTCCT‑‑CGCGGTTGCTGGAGCA‑GGCGGCCACGGCAAGCAGCGCGGCGCACAACAGGGTAATGCGAGGCAATGATTTCATGATGaaa                                                                    <  2:338314/149‑1 (MQ=255)                                                                      ccGACGGGCTCTTGGCTTCTTGCAGCCGGCTGTAGTCTTTGAGGAAGCCGGAGTAATCCTTGGGGTCCAAGCGCGGTTGCTGGAGCA‑GGCGGCCACGGCAAGCAGCGCGGCGCACAACAGGGTAATGCGAGGCAATGATTTCatgatg                                                                       >  1:115915/1‑148 (MQ=255)                                                                       cGACGGGCTCTTGGCTTCTTGCAGCCGGCTGTAGTCTTTGAGGAAGCCGGAGTAATCCTTGGGGTCCA‑‑CGCTGTTGCTGGAGCA‑GGCGGCCACGGCAAGCAGCGCGGCGCACAACAGGGTAATGCGAGGCAATGATTTCATGATGaaaa                                                                   <  5:184531/149‑1 (MQ=255)                                                                          cGGGCTCTTGGCTTCTTGCAGCCGGCTGTAGTCTTTGAGGAAGCCGGAGTAATCCTTGCGGGCCA‑‑CGCGGTTGCTGGAGCA‑GGCGGCCACGGCAAGCAGCGCGGCGCACAACAGGGTAATGCGAGGCAATGATTTCATGATGAAAAGGc                                                                <  6:306502/149‑1 (MQ=255)                                                                                 ttGGCTTCTTGCAGCCGGCTGTAGTCTTTGAGGAAGCCGGAGTAATCCTTGGGGTCCA‑‑CGCGGTTGCTGGAGCA‑GGCGGCCACGGCAAGCAGCGCGGCGCACAACAGGGTAATGCGAGGCAATGATTTCATGATGAAAAGGCATCCg                                                           <  5:301888/147‑1 (MQ=255)                                                                                   ggCTTCTTGCAGCCGGCTGTAGTCTTTGAGGAAGCCGGAGTAATCCTTGGGGTCCA‑‑CGCGGTTGCTGGAGCA‑GGCGGCCACGGCAAGCAGCGCGGCGCACAACAGGGTAATGCGAGGCAATGATTTCATGATGAAAAGGCATCCGt                                                          >  6:164535/1‑146 (MQ=255)                                                                                               ccGGCTGTAGTCTTTGAGGAAGCCGGAGTAATCCTTGGGGTCCA‑‑CGCGGTTGCTGGAGCA‑GGCGGCCACGGCAAGCAGCGCGGCGCACAACAGGGTAATGCGAGGCAATGATTTCATGATGAAAAGGCTTCCGTGAGAGGTCGGCCTc                                            >  3:258117/1‑148 (MQ=255)                                                                                               ccGGCTGTAGTCTTTGAGGAAGCCGGAGTAATCCTTGGGGTCCA‑‑CGCGGTTGCTGGAGCA‑GGCGGCCACGGCAAGCAGCGCGGCGCACAACAGGGTAATGCGAGGCAATGATTTCATGATGAAAAGGCATCCGTGAGAGGTCGGCCTc                                            >  5:230739/1‑148 (MQ=255)                                                                                                cGGCTGTAGTCTTTGAGGAAGCCGGAGTAATCCTTGGGGTCCA‑‑CGCGGTTGCTGGAGCA‑GGCGGCCACGGCAAGCAGCGCGGCGCACAACAGGGTAATGCGAGGCAATGATTTCATGATGAAAAGGCATCCGTGAGAGGTCGGCCTcaa                                          <  1:192742/149‑1 (MQ=255)                                                                                                  gCTGTAGTCTTTGAGGAAGCCGGAGTAATCCTTGGGGTCCA‑‑CGCGGTTGCTGGAGCA‑GGCGGCCACGGCAAGCAGCGCGGCGCACAACAGGGTAATGCGAGGCAATGATTTCATGATGAAAAGGCATCCGTGAGAGGTCGGCCTcaac                                         >  4:94193/1‑148 (MQ=255)                                                                                                      tAGTCTTTGAGGAAGCCGGAGTAATCCTTGGGGTCCA‑‑CGCGGTTGCTGGAGCA‑GGCGGCCACGGCAAGCAGCGCGGCGCACAACAGGGTAATGCGAGGCAATGATTTCATGATGAAAAGGCATCCGTGAGAGGTCGGCCTCaacaac                                      <  6:230739/147‑1 (MQ=255)                                                                                                      tAGTCTTTGAGGAAGCCGGAGTAATCCTTGGGGTCCA‑‑CGCGGTTGCTGGAGCA‑GGCGGCCACGGCAAGCAGCGCGGCGCACAACAGGGTAATGCGAGGCAATGATTTCATGATGAAAAGGCATCCGTGAGAGGTCGGCCTCAACAACGc                                    >  4:226053/1‑149 (MQ=255)                                                                                                         tCTTTGAGGAAGCCGGAGTAATCCTTGGGGTCCA‑‑CGCGGTTGCTGGAGCA‑GGCGGCCACGGCAAGCAGCGCGGCGCACAACAGGGTAATGCGAGGCATTGATTTCATCATGAAAAGGCATCCGTGACAGGACGGCCTCAACAACGc                                    >  6:384817/1‑146 (MQ=255)                                                                                                          cTTTGAGGAAGCCGGAGTAATCCTTGGGGTCCA‑‑CGCGGTTGCTGGAGCA‑GGCGGCCACGGCAAGCAGCGCGGCGCACAACAGGGTAATGCGAGGCAATGATTTCATGATGAAAAGGCATCCGTGAGAGGTCGGCCTCAACAACGCTAg                                 <  7:207139/148‑1 (MQ=255)                                                                                                           tttGAGGAAGCCGGAGTAATCCTTGGGGTCCA‑‑CGCGGTTGCTGGAGCA‑GGCGGCCACGGCAAGCAGCGCGGCGCACAACAGGGTAATGCGAGGCAATGATTTCATGATGAAAAGGCATCCGTGAGAGGTCGGCCTCAACAACGCTAg                                 >  6:376198/1‑147 (MQ=255)                                                                                                             tGAGGAAGCCGGAGTAATCCTTGGGGTCCA‑‑CGCGGTTGCTGGAGCA‑GGCGGCCACGGCAAGCAGCGCGGCGCACAACAGGGTAATGCGAGGCAATGATTTCATGATGAAAAGGCATCCGTGAGAGGTCGGCCTCAACAACGCTAGCAg                              <  4:205961/148‑1 (MQ=255)                                                                                                               aGGAAGCCGGAGTAATCCTTGGGGTCCA‑‑CGCGGTTGCTGGAGCA‑GGCGGCCACGGCAAGCAGCGCGGCGCACAACAGGGTAATGCGAGGCAATGATTTCATGATGAAAAGGCATCCTTGagag                                                     <  4:341470/123‑1 (MQ=255)                                                                                                               aGGAAGCCGGAGTAATCCTTGGGGTCCA‑‑CGCGGTTGCTGGAGCA‑GGCGGCCACGGCAAGCAGCGCGGCGCACAACAGGGTAATGCGAGGCAATGATTTCATGATGAAAAGGCATCCGTGagag                                                     >  3:341470/1‑123 (MQ=255)                                                                                                                      cGGAGTAATCCTTGGGGTCCA‑‑CGCGGTTGCTGGAGCA‑GGCGGCCACGGCAAGCAGCGCGGCGCACAACAGGGTAATGCGAGGCAATGATTTCATGATGAAAAGGCATCCGTGAGAGGTCGGCCTCAACAACGCTAGCAGGTGCGAg                       <  4:47746/146‑1 (MQ=255)                                                                                                                      cGGAGTAATCCTTGGGGTCCA‑‑CGCGGTTGCTGGAGCA‑GGCGGCCACGGCAAGCAGCGCGGCGCACAACAGGGTAATGCGAGGCAATGATTTCATGATGAAAAGGCATCCGTGAGAGGTCGGCCTCAACAACGCTAGCAGGTGCGAg                       >  3:47746/1‑146 (MQ=255)                                                                                                                             aTCCTTGGGGTCCA‑‑GGCGGTGGCTGGAGGA‑GGGGGCGACGGCAAGCAGCGCGGCGTACAACAGGGTAATGCGAGGCAATGTTTTCATGATGAAAAGGCATCCGTGTGAGGTCGGCCTGAACTACGGTAGCTGGTGCGAGGCAGTTGGc              <  8:287725/148‑1 (MQ=255)                                                                                                                             aTCCTTGGGGTCCA‑‑CGCGGTTGCTGGAGCA‑GGCGGCCACGGCAAGCAGCGCGGCGCACAACAGGGTAATGCGAGGCAATGATTTCATGATGAAAAGGCATCCGTGAGAGGTCGGCCTCAACAACGCTAGCAGGTGCGAGGCAGTTGGc              <  4:258117/148‑1 (MQ=255)                                                                                                                                   ggggTCCA‑‑CGCGGTTGCTGGAGCA‑GGCTGACACGGCAAGCAGCGCGTCGCACAACAGTTTAATGCGATGCAATGATTTCATGATGAAAAGGCATCCTTGAGAGTTCTTCCTCAACAACTCTAGCAGGTGCGATGCAGTTGTCCAGTg         <  2:67/147‑1 (MQ=255)                                                                                                                                    gggTCCA‑‑CGCGGTTGCTGGAGCA‑GGCGGCCACGGCAAGCAGCGCGGCGCACAACAGGGTAATGCGAGGCAATGATTTCATGATGAAAAGGCATCCGTGAGAGGTCGGCCTCAACAACGCTAGCAGGTGCGAGGCAGTTGGCCAGTGGAc      >  2:233021/1‑149 (MQ=255)                                                                                                                                         cA‑‑CGCGGTTGCTGGAGCA‑GGCGGCCACGGCAAGCAGCGCGGCGCACAACAGGGTAATGCGAGGCAATGATTTCATGATGAAAAGGCATCCGTGAGAGGTCGGCCTCAACAACGCTAGCAGGTGCGAGGCAGTTGGCCAGTGGACgggg  <  5:376198/148‑1 (MQ=255)                                                                                                                                         cA‑‑CGCGGTTGCTGGAGCA‑GGCGGCCACGGCAAGCAGCGCGGCGCACAACAGGGTAATGCGAGGCAATGATTTCATGATGAAAAGGCATCCGTGAGAGGTCGGCCTCAACAACGCTAGCAGGTGCGAGGCAGTTGGCCAGTGGACgggg  <  2:115915/148‑1 (MQ=255)                                                                                                                                          a‑‑CGCGGTTGCTGGAGCA‑GGCGGCCACGGCAAGCAGCGCGGCGCACAACAGGGTAATGCGAGGCAATGATTTCATGATGAAAAGGCATCCGTGAGAGGTCGGCCTCAACAACGCTAGCAGGTGCGAGGCAGTTGGCCAGTGGACgg    >  2:343929/1‑145 (MQ=255)                                                                                                                                          a‑‑CGCGGTTGCTGGAGCA‑GGCGGCCACGGCAAGCAGCGCGGCGCACAACAGGGTAATGCGAGGCAATGATTTCATGATGAAAAGGCATCCGTGAGAGGTCGGCCTCAACAACGCTAGCAGGTGCGAGGCAGTTGGCCAGTGGACgg    <  1:343929/145‑1 (MQ=255)                                                                                                                                          |                                                                                                                                                      GCTCGACGAACACCTTGCTGTACTGGTTGATATTGACCTTGGGGTCGATCCAGCGCATCACCGGGTCACCCGACGGGCTCTTGGCTTCTTGCAGCCGGCTGTAGTCTTTGAGGAAGCCGGAGTAATCCTTGGGGTCCA‑‑CGCGGTTGCTGGAGCA‑GGCGGCCACGGCAAGCAGCGCGGCGCACAACAGGGTAATGCGAGGCAATGATTTCATGATGAAAAGGCATCCGTGAGAGGTCGGCCTCAACAACGCTAGCAGGTGCGAGGCAGTTGGCCAGTGGACGGGG  >  NC_002947/3562177‑3562460

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