Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,882,451	C→G	5.6%	S13R (AGC→AGG)	PP_4291 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,882,451	0	C	G	5.6%	223.4 / 3.7	71	S13R (AGC→AGG)	PP_4291	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (33/34);  new base G (2/2);  total (35/36)
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.68e-01

TTCACCTCTGCAGTAACCGCTACGCCGTCGACCCGCCGGAATGGGTCGCCCGGGGCGCCTGCGCGCCAAAAAAACGCCAACAAAAACAATAAGTCCGGGAAT‑CACAACATGAAGCGCATCACCACGTCCCTGTTGTTGGGCAGCAGCCTGCT‑GGCCACCCTACCCTCGCACGCAGGCGAATGGCTGCAATGGCACGGCGAAAGCCTGACCTACCTGTACGGCAAGGACTTCAAGGTCAACCCTGATATCCAGCAGACCATCACCTTCGAGCACGCCAACAAATGGAAGTACGGC  >  NC_002947/4882305‑4882598                                                                                                                                                    |                                                                                                                                                     ttCACCTCTGCAGTAACCGCTACGCCGTCGACCCGCCGGAATGGGTCGCCCGGGGCGCCTGCGCGCCAAAAAAACGCCAACAAAAACAATAAGTCCGGGAAT‑CACAACATGAAGCGCATCACCACGTCCCTGTTGTTGGGCAGCAGCct                                                                                                                                                    >  4:152840/1‑149 (MQ=255)    aCCTCTGCAGTAACCGCTACGCCGTCGACCCGCCGGAATGGGTCGCCCGGGGCGCCTGCGCGCCAAAAAAACGCCAACAAAAACAATAAGTCCGGGAAT‑CACAACATGAAGCGCATCACCACGTCCCTGTTGTTGGGCAGCAGCctgct                                                                                                                                                 >  4:192383/1‑149 (MQ=255)      ctctGCAGTAACCGCTACGCCGTCGACCCGCCGGAATGGGTCGCCCGGGGCGCCTGCGCGCCAAAAAAACGCCAACAAAAACAATAAGTCCGGGAAT‑CACAACATGAAGCGCATCACCACGTCCCTGTTGTTGGGCAGCAGCCTGCT‑gg                                                                                                                                              >  2:256914/1‑149 (MQ=255)           cAGTAACCGCTACGCCGTCGACCCGCCGGAATGGGTCGCCCGGGGCGCCTGCGCGCCAAAAAAACGCCAACAAAAACAATAAGTCCGGGAAT‑CACAACATGAAGCGCATCACCACGTCCCTGTTGTTGGGCAGCAGCCTGCT‑GGCCac                                                                                                                                          <  7:208958/148‑1 (MQ=255)           cAGTAACCGCTACGCCGTCGACCCGCCGGAATGGGTCGCCCGGGGCGCCTGCGCGCCAAAAAAACGCCAACAAAAACAATAAGTCCGGGAAT‑CACAACATGAAGCGCATCACCACGTCCCTGTTGTTGGGCAGCAGCCTGCT‑GGCCa                                                                                                                                           <  5:204624/147‑1 (MQ=255)           cAGTAACCGCTACGCCGTCGACCCGCCGGAATGGGTCGCCCGGGGCGCCTGCGCGCCAAAAAAACGCCAACAAAAACAATAAGTCCGGGAAT‑CACAACATGAAGCGCATCACCACGTCCCTGTTGTTGGCCAGCAGGCTGCT‑Gccc                                                                                                                                            <  8:33543/146‑1 (MQ=255)           cAGTAACCGCTACGCCGTCGACCCGCCGGAATGGGTCGCCCGGGGCGCCTGCGCGCCAAAAAAACGCCAACAAAAACAATAAGTCCGGGAAT‑CACAACATGAAGCGCATCACCACGTCCCTGTTGTTGGCCAGCAGGCTGCT‑GCCCaa                                                                                                                                          <  1:216559/148‑2 (MQ=255)                  cGCTACGCCGTCGACCCGCCGGAATGGGTCGCCCGGGGCGCCTGCGCGCCAAAAAAACGCCAACAAAAACAATAAGTCCGGGAAT‑CACAACATGAAGCGCATCACCACGTCCCTGTTGTTGGGCAGCAGCCTGCT‑GGCCaccctaccct                                                                                                                                  >  5:101121/1‑149 (MQ=255)                   gCTACGCCGTCGACCCGCCGGAATGGGTCGCCCGGGGCGCCTGCGCGCCAAAAAAACGCCAACAAAAACAATAAGTCCGGGAAT‑CACAACATGAAGCGCATCACCACGTCCCTGTTGTTGGGCAGCAGCCTGCT‑GGCCAccctcccctc                                                                                                                                 >  7:224413/1‑149 (MQ=255)                   gCTACGCCGTCGACCCGCCGGAATGGGTCGCCCGGGGCGCCTGCGCGCCAAAAAAACGCCAACAAAAACAATAAGTCCGGGAAT‑CACAACATGAAGCGCATCACCACGTCCCTGTTGTTGGGCAGCAGCCTGCT‑GGCCACCCTACCCTc                                                                                                                                 <  8:296350/149‑1 (MQ=255)                     tACGCCGTCGACCCGCCGGAATGGGTCGCCCGGGGCGCCTGCGCGCCAAAAAAACGCCAACAAAAACAATAAGTCCGGGAAT‑CACAACATGAAGCGCATCACCACGTCCCTGTTGTTGGGCAGCAGCCTGCT‑GGCCACCCTACCCTcgc                                                                                                                               >  5:290957/1‑149 (MQ=255)                     tACGCCGTCGACCCGCCGGAATGGGTCGCCCGGGGCGCCTGCGCGCCAAAAAAACGCCAACAAAAACAATAAGTCCGGGAAT‑CACAACATGAAGCGCATCACCACGTCCCTGTTGTTGGGCAGCAGCCTGCT‑GGCCACCCTACCCTcgc                                                                                                                               >  5:235216/1‑149 (MQ=255)                           gTCGACCCGCCGGAATGGGTCGCCCGGGGCGCCTGCGCGCCAAAAAAACGCCAACAAAAACAATAAGTCCGGGAAT‑CACAACATGAAGCGCATCACCACGTCCCTGTTGGTGGGCAGCAGCCTGCT‑GGCCACCCTACCCTCGCACGCAg                                                                                                                         >  2:244861/1‑149 (MQ=255)                                 ccgccgGAATGGGTCGCCCGGGGCGCCTGCGCGCCAAAAAAACGCCAACAAAAACAATAAGTCCGGGAAT‑CACAACATGAAGCGCATCACCACGTCCCTGTTGTTGGGCAGCAGCCTGCT‑GGCCACCCTACCCTCGCACGCAGGCGAAt                                                                                                                   <  3:152840/149‑1 (MQ=255)                                 ccgccgGAATGGGTCGCCCGGGGCGCCTGCGCGCCAAAAAAACGCCAACAAAAACAATAAGTCCGGGAAT‑CACAACAAGAAGCGCATCACAACGGCCGTGTTGTTGGGCAGCAGCCAGCA‑GGCCACCCTACCCTCCCACGCAGGCGAAt                                                                                                                   >  4:314341/1‑149 (MQ=255)                                      ggAATGGGTCGCCCGGGGCGCCTGCGCGCCAAAAAAACGCCAACAAAAACAATAAGTCCGGGAAT‑CACAACATGAAGCGCATCACCACGTCCCTGTTGTTGGGCAGCAGCCTGCT‑GGCCACCCTACCCTCGCACGCAGGCGAATGGCTg                                                                                                              >  2:249523/1‑149 (MQ=255)                                         aTGTGTCGCCCGGGGCGCCTGCGCGCCAAAAAAACGACTACAAAAACAATAAGTCCGGGATT‑CACAACATGAAGCGCATCACCACGTCCCTGTTGTTGGGCAGCAGCCTGCT‑GGCCACCCTACCCTCGCACGCAGGCGAATGGCTGCaa                                                                                                           <  7:178257/149‑1 (MQ=255)                                         aTGGGTCGCCCGGGGCGCCTGCGCGCCAAAAAAACGCCAACAAAAACAATAAGTCCGGGAAT‑CACAACATGAAGCGCATCACCACGTCCCTGTTGTTGGGCAGCAGCCTGCT‑GGCCACCCTACCCTCGCACGCAGGCGAATGGCTGCaa                                                                                                           <  1:100485/149‑1 (MQ=255)                                              tCGCCCGGGGCGCCTGCGCGCCAAAAAAACGCCAACAAAAACAATAAGTCCGGGAAT‑CACAACATGAAGCGCATCACCACGTCCCTGTTGTTGGGCAGCAGCCTGCT‑GGCCACCCTACCCTCGCACGCAGGCGAATGGCTGCAATGGCa                                                                                                      >  3:73142/1‑149 (MQ=255)                                                  ccGGGGCGCCTGCGCGCCAAAAAAACGCCAACAAAAACAATAAGTCCGGGAAT‑CACAACATGAAGCGCATCACCACGTCCCTGTTGTTGGGCAGCAGCCTGCT‑GGCCACCCTACCCTCGCACGCAGGCGAAt                                                                                                                   >  6:44153/1‑132 (MQ=255)                                                  ccGGGGCGCCTGCGCGCCAAAAAAACGCCAACAAAAACAATAAGTCCGGGAAT‑CACAACATGAAGCGCATCACCACGTCCCTGTTGTTGGGCAGCAGCCTGCT‑GGCCACCCTACCCTCGCACGCAGGCGAAt                                                                                                                   <  5:44153/132‑1 (MQ=255)                                                  ccGGGGCGCCTGCGCGCCAAAAAAACGCCAACAAAAACAATAAGTCCGGGAAT‑CACAACATGAAGCGCATCACCACGTCCCTGTTGTTGGGCAGCAGCCTGCT‑GGCCACCCTACCCTCGCACGCAGGCGAATGGCTGCAATGGCACGGc                                                                                                  <  2:71754/149‑1 (MQ=255)                                                     gggCGCCTGCGCGCCAAAAAAACGCCAACAAAAACAATAAGTCCGGGAATGCACAACATGAAGCGCATCACCACGTCCCTGTTGTTGGGCAGCAGCCTGCT‑GGCCACCCTACCCTCGCACGCAGGCGAATGGCTGCAATGGCACGGCg                                                                                                 >  1:243990/1‑148 (MQ=255)                                                     gggCGCCTGCGCGCCAAAAAAACGCCAACAAAAACAATAAGTCCGGGAAT‑CACAACATGAAGCGCATCACCACGTCCCTGTTGTTGGGCAGCAGCCTGCT‑GGCCACCCTACCCTCGCACGCAGGCGAATGGCTGCAATGGCACGGCGa                                                                                                >  6:262740/1‑148 (MQ=255)                                                      ggCGCCTGCGCGCCAAAAAAACGCCAACAAAAACAATAAGTCCGGGAAT‑CACAACATGAAGCGCATCACCACGTCCCTGTTGTTGGGCAGCAGCCTGCT‑GGCCACCCTACCCTCGCACGCAGGCGAATGGCTGCAATGGCACGGCGaaa                                                                                              <  1:317046/149‑1 (MQ=255)                                                          ccTGCGCGCCAAAAAAACTCCAACAAAAACAATCTGTCCGTGAAT‑CACAACATGAAGCGCATCACCACGTCCCTGTTGTTGGGCAGCAGCCTGCT‑GGCCACCCTACCCTCGCACGCAGGCGAATGGCTGCAATGGCACGGCGAAAGCCt                                                                                          <  2:222276/149‑1 (MQ=255)                                                              cgcgcCAAAAAAACGCCAACAAAAACAATAAGTCCGGGAAT‑CACAACATGAAGCGCATCACCACGTCCCTGTTGTTGGGCAGCAGCCTGCT‑GGCCACCCTACCCTCGCACGCAGGCGAATGGCTGCTATGGCACGGCGAAAGCCTGa                                                                                        <  1:249523/147‑1 (MQ=255)                                                              cgcgcCAAAAAAACGCCAACAAAAACAATAAGTCCGGGAAT‑CACAACATGAAGCGCATCACCACGTCCCTGTTGTTGGGCAGCAGCCTGCT‑GGCCACCCTACCCTCGCACGCAGGCGAATGGCTGCAATGGCACGGCGAAAGCCTGac                                                                                       <  6:235216/148‑1 (MQ=255)                                                                 gcCAAAAAAACGCCAACAAAAACAATAAGTCCGGGAAT‑CACAACATGAAGCGCATCACCACGTCCCTGTTGTTGGGCAGCAGCCTGCT‑GGCCACCCTACCCTCGCACGCAGGCGAATGGCTGCAATGGCACGGCGAAAGCCTGacctac                                                                                   <  2:279701/149‑1 (MQ=255)                                                                  ccAAAAAAACGCCAACAAAAACAATAAGTCCGGGAAT‑CACAACATGAAGCGCATCACCACGTCCCTGTTGTTGGGCAGCAGCCTGCTCGGCCACCCTACCCTCGCACGCAGGCGAATGGCTGCAATGGCACGGCGAAAGCCTGacctac                                                                                   <  1:256914/149‑1 (MQ=255)                                                                     aaaaaaCGCCAACAAAAACAATAAGTCCGGGAAT‑CACAACATGAAGCGCATCACCACGTCCCTGTTGTTGGGCAGCAGCCTGCT‑GGCCACCCTACCCTCGCACGCAGGCGAATGGCTGCAATGGCACGGCGAAAGCCTGACCTACCTg                                                                                >  4:312846/1‑148 (MQ=255)                                                                     aaaaaaCGCCAACAAAAACAATAAGTCCGCGAAT‑CACAACATGAAGCGCATCACCACGTCCCTGATCATGGGCAGCAGCCTGCT‑GGCCACCCTACCCACGCACCCACGCGAATGGCTGCAATGGCACGGCGAAAGCCTGACCTACCTg                                                                                >  3:261645/1‑148 (MQ=255)                                                                       aaaaCGCCAACAAAAACAATAAGTCCGGGAAT‑CACAACATGAAGCGCATCACCACGTCCCTGTTGTTGGGCAGCAGCCTGCT‑GGCCACCCTACCCTCGCACGCAGGCGAATGGCTGCAATGGCACGGCGAAAGCCTGACCTACCTGTAc                                                                             <  3:192383/149‑1 (MQ=255)                                                                         aaCGCCAACAAAAACAATAAGTCCGGGAAT‑CACAACATGAAGCGCATCACCACGTCCCTGTTGTTGGGCAGCAGCCTGCT‑GGCCACCCTACCCTCGCACGCAGGCGAATGGCTGCAATGGCACGGCGAAAGCCTGACCTACCTg                                                                                >  6:30105/1‑144 (MQ=255)                                                                         aaCGCCAACAAAAACAATAAGTCCGGGAAT‑CACAACATGAAGCGCATCACCACGTCCCTGTTGTTGGGCAGCAGCCTGCT‑GGCCACCCTACCCTCGCACGCAGGCGAATGGCTGCAATGGCACGGCGAAAGCCTGACCTACCTg                                                                                <  5:30105/144‑1 (MQ=255)                                                                            gCCAACAAAAACAATAAGTCCGGGAAT‑CACAACATGAAGCGCATCACCACGTCCCTGTTGTTGGGCAGCAGCCTGCT‑GGCCACCCTACCCTCGCACGCAGGCGAATGGCTGCAATGGCACGGCGAAAGCCTGACCTACCTGTACGGCaa                                                                        <  2:241922/149‑1 (MQ=255)                                                                             cgatcaaAAAGAATAAGTCCGGGAAT‑CACTACATGAAGCGGATGACCACGTCCCTGTTGTTGGGCAGCAGCCTGCT‑GGCCACCCTACCCTCGCACTCATGCGAATGGCTGCAATGTCACTGCGTAAGCCTGACCTACCTGTAGGGCAAg                                                                       <  4:273209/145‑1 (MQ=255)                                                                             ccAACAAAAACAATAAGTCCGGGAAT‑CACAACATGAAGCGCATCACCACGTCCCTGTTGTTGGGCAGCAGCCTGCT‑GGCCACCCTACCCTCGCACGCAGGCGAATGGCTGCAATGGCACGGCGAAATCCTGACCTACCTGTACGGCAAg                                                                       <  6:290957/149‑1 (MQ=255)                                                                             ccAACAAAAACAATAAGTCCGGGAAT‑CACAACATGAAGCGCATCACCACGTCCCAGTAGTTGGGCAGCACCCTGCT‑GGCCAACCTACCCTCGCACGCAGCCCAATGGCAGCAAAGGAACGGCGAAAGCCTGACCTACCAGTACGGCAAg                                                                       >  7:189662/1‑149 (MQ=255)                                                                              caacaaAAACAATAAGTCCGGGAAT‑CACAACATGAAGCGCATCACCACGTCCCTGTTGTTGGGCAGCAGCCTGCT‑GGCCACCCTACCCTCGCACGCAGGCGAATGGCTGCAATGGCACGGCGAAAGCCTGACCTGCCTGTACGGCAAg                                                                       >  8:248521/1‑148 (MQ=255)                                                                                acaaAAACAATAAGTCCGGGAAT‑CACAACATGAAGCGCATCACCCCGTCCCAGTTGTTGGGCAGCAGCCAGCT‑GGCCCCCCTCCCCTCGCACGCAGGCGAATGGCTGCAATGGCACGGCGAAAGCCAGACCTACCTGTACGGCAAggcc                                                                    >  1:135710/1‑147 (MQ=255)                                                                                      aCAATAAGTCCGGGAAT‑CACAACATGAAGCGCATCACCACGTCCCTGTTGGTGGGCAGCAGCCTGCT‑GGCCACCCTACCCTCGCACGCAGGCGAATGGCTGCAATGGCACGGCGAAAGCCTGACCTACCTGTACGGCAAGGACTTCAAg                                                              <  1:244861/149‑1 (MQ=255)                                                                                             gTCCGGGAAT‑CACGACATGAAGCGCATCACCACGTCCCTGTTGTTGGGCAGCAGCCTGCT‑GGCCACCCTACCCTCGCACGCAGGCGAATGGCTGCAATGGCATGGCGAAAGCCTGACCTACCTGTACGGCAAGGACTTAAAGGTAAAcc                                                       >  5:232512/1‑149 (MQ=255)                                                                                                cGGGAAT‑CACAACATGAAGCGCATCACCACGTCCCTGTTGTTGGGCAGCAGCCTGCT‑GGCCACCCTACCCTCGCACGCAGGCGAATGGCTGCAATGGCACGGCGAAAGCCTGACCTACCTGTACGGCAAGGACTTCAAGGTCAACCCt                                                     >  2:161529/1‑148 (MQ=255)                                                                                                      t‑CACAACATGAAGCGCATCACCACGTCCCTGTTGTTGGGCAGCAGCCTGCT‑GGCCACCCTACCCTCGCACGCAGGCGAATGGCTGCAATGGCACGGCGAAAGCCTGACCTACCTGTACg                                                                            <  7:56183/119‑1 (MQ=255)                                                                                                      t‑CACAACATGAAGCGCATCACCACGTCCCTGTTGTTGGGCAGCAGCCTGCT‑GGCCACCCTACCCTCGCACGCAGGCGAATGGCTGCAATGGCACGGCGAAAGCCTGACCTACCTGTACg                                                                            >  8:56183/1‑119 (MQ=255)                                                                                                         acaacaTGAAGCGCATCACCACGTCCCTGTTGTTGGGCAGCAGCCTGCT‑GGCCACCCTACCCTCGCACGCAGGCGAATGGCTGCAATGGCACGGCGAAAGCCTGACCTACCTGTACGGCAAGGACTTCAAGGTCAACCCTGATATCcag                                            <  6:99972/149‑1 (MQ=255)                                                                                                         acaacaTGAAGCGCATCACCACGTCCCTGTTGTTGGGCAGCAGCCTGCT‑GGCCACCCTACCCTCGCACGCAGGCGAATGGCTGCAATGGCACGGCGAAAGCCTGACCTACCTGTACGGCAAGGACTTCAAGGTCAACCCTGATATCca                                             <  5:199981/148‑1 (MQ=255)                                                                                                         acaacaTGAAGCGCATCACCACGTCCATGTTGTTGGGCAGCAGCCTGCT‑TGCCACCCTACCCTCGCACGCAGGCGAATGGCTGCAATGGCACGGCGAAAGCCTGACCTACCTGTACGGCAAGGACTTCAAGGTCAACCCTGATATCca                                             <  4:76709/148‑1 (MQ=255)                                                                                                          caacaTGAAGCGCATCACCACGTCCCTGTTGTTGGGCAGCAGCCTGCT‑GGCCACCCTACCCTCGCACGCAGGCGAATGGCTGCAATGGCACGGCGAAAGCCTGacctac                                                                                   <  5:49789/109‑1 (MQ=255)                                                                                                          caacaTGAAGCGCATCACCACGTCCCTGTTGTTGGGCAGCAGCCTGCT‑GGCCACCCTACCCTCGCACGCAGGCGAATGGCTGCAATGGCACGGCGAAAGCCTGacctac                                                                                   >  6:49789/1‑109 (MQ=255)                                                                                                          caacaTGAAGCGCATCACCACGTCCCTGTTGTTGGGCAGCAGCCTGCT‑GGCCACCCTACCCTCGCACGCAGGCGAATGGCTGCAATGGCACGGCGAAAGCCTGACCTACCTGTACGGCAAGGACTTCAAGGTCAACCCTGATATCcagc                                           <  1:161529/149‑1 (MQ=255)                                                                                                           aacaTGAAGCGCATCACCACGTCCCTGTTGTTGGGCAGCAGCCTGCT‑GGCCACCCTACCCTCGCACGCAGGCGAATGGCTGCAATGGCACGGCGAAAGCCTGACCTACCTGTACGGCAAGGACTTCAAGGTCAACCCTGATATCcag                                            >  4:179879/1‑147 (MQ=255)                                                                                                              aTGAAGCGCATCACCACGTCCCTGTTGTTGGGCAGCAGCCTGCT‑GGCCACCCTACCCTCGCACGCAGGCGAATGGCTGCAATGGCACGGCGAAAGCCTGACCTACCTGTACGGCAAGGACTTCAAGGTCAACCCTGATATCCAGCAGAc                                       <  3:158164/149‑1 (MQ=255)                                                                                                                gAAGCGCATCACCACGTCCCTGTTGTTGGGCAGCAGCCTGCT‑GGCCACCCTACCCTCGCACGCAGGCGAATGGCTGCAATGGCACGGCGAAAGCCTGACCTACCTGTACGGCAAGGACTTCAAGGTCAACCCTGATATCCAGCAGACCa                                     >  3:247631/1‑149 (MQ=255)                                                                                                                        tCACCACGTCCCTGTTGTTGGGCAGCAGCCTGCT‑GGCCACCCTACCCTCGCACGCAGGCGAATGGCTGCAATGGCACGGCGAAAGCCTGACCTACCTGTACGGCAAGGACTTCAAGGTCAACCCTGATATCCAGCAGACCATCACCTTc                             >  2:66337/1‑149 (MQ=255)                                                                                                                        tCACCAAGTCCCTGTTGTTGGGCAGCAGCCTGCT‑GGCCACCCTACCCTCGCACGCAGTCGAATGGGTGCAATGTCACGGCGAAAGCCTTACCTACCTGTACGGCAAGGACTTCAAGGTCAACCCTGATATCCAGCAGAGCATCACCTTc                             <  4:261645/149‑1 (MQ=255)                                                                                                                         caccacGTCCCTGTTGTTGGGCAGCAGCCTGCT‑GGCCACCCTACCCTCGCACGCAGGCGAATGGCTGCAATGGCACGGCGAAAGCCTGACCTACCTGTACGGCAAGGACTTCAAGGTCAACCCTGATATCCAGCAGACCATCACCTTCg                            >  5:236643/1‑149 (MQ=255)                                                                                                                           ccacGTCCCTGTTGTTGGGCAGCATCCTGCT‑GGCCACCCTACCCTCGCACGCAGGCGAATGGCTGCAATGGCACGGCGAAAGCCTGACCTACCTGTACGGCAAGGACTTCAAGGTCAACCCTGATATCCAGCAGACCATCACCTTCGAg                          >  7:200768/1‑149 (MQ=255)                                                                                                                                 cccTGTTGTTGGGCAGCAGCCTGCT‑GGCCACCCTACCCTCGCACGCAGGCGAATGGCTGCAATGGCACGGCGAAAGCCTGACCTACCTGTACGGCAAGGACTTCAAGGTCAACCCTGATATCCAGCAGACCATCACCTTCGAGCACGc                     >  6:134246/1‑148 (MQ=255)                                                                                                                                  ccTGTTGTTGGGCAGCAGCCTGCT‑GGCCACCCTACCCTCGCACGCAGGCGAATGGCTGCAATGGCACGGCGAAAGCCTGACCTACCTGTACGGCAAGGACTTCAAGGTCAACCCTGATATCCAGCAGACCATCACCTTCGAGCACGcc                    <  2:243990/148‑1 (MQ=255)                                                                                                                                  ccTGTTGTTGGGCAGCAGCCTGCT‑GGCCACCCTACCCTCGCACGCAGGCGAATGGCTGCAATGGCACGGCGAAAGCCTGACCTACCTGTACGGCAAGGACTTCAAGGTCAACCCTGATATCCAGCAGACCATCACCTTCGAGCACGcc                    >  6:20380/1‑148 (MQ=255)                                                                                                                                   cTGTTGTTGGGCAGCAGCCTGCT‑GGCCACCCTACCCTCGCACGCAGGCGAATGGCTGCAATGGCACGGCGAAAGCCTGACCTACCTGTACGGCAAGGACTTCAAGGTCAACCCTGATATCCAGCAGACCATCACCTTCGAGCACGCca                   >  5:248701/1‑148 (MQ=255)                                                                                                                                     tagggtgGCCAGCAGGCTGCT‑GCCCACCCTACCCTCGCACGCAGCCGAATGGCTGCAATGGCACGGCGAAAGCCTGACCTACCTGTACGGCAAGGACTTCAAGGTCAACCCAGATATCCAGCAGACCATCACCTTCGAGCACGCcaac                 >  2:315217/6‑148 (MQ=255)                                                                                                                                         ttgGGCAGCAGCCTGCT‑GGCCACCCTACCCTCGCACGCAGGCGAATGGCTGCAATGGCACGGCGAAAGCCTGACCTACCTGTACGGCAAGGACTTCAAGGTCAACCCTGATATCCAGCAGACCATCACCTTCGAGCACGCCAACAAAt             <  1:66337/148‑1 (MQ=255)                                                                                                                                         ttgGGCAGCAGCCTGCT‑GGCCACCCTACCCTCGCACGCAGGCGAATGGCTGCAATGGCACGGCGAAAGCCTGACCTACCTGTACGGCAAGGACTTCAAGGTCAACCCTGATATCCAGCAGACCATCACCTTCGAGCACGCCAACAAATg            >  2:103700/1‑149 (MQ=255)                                                                                                                                          tgGCCAGCAGGCTGCT‑GCCCACCCTACCCTCGCACGCAGGCGAATGGCTGCAATGGCACGGCGAAAGCCTGACCTACCTGTACGGCAAGGACTTCAAGGTCAACCCTGATATCCAGCAGACCATCACCTTCGAGCACGCCAACAAATg            >  6:13383/1‑148 (MQ=255)                                                                                                                                           gggCAGCAGCCTGCT‑GGCCACCCTACCCTCGCACGCAGGCGAATGGCTGCAATGGCACGGCGAAAGCCTGACCTACCTGTACGGCAAGGACTTCAAGGTCAACCCTGATATCCAGCAGACCATCACCTTCGAGCACGCCAACAAATgg           <  6:236643/148‑1 (MQ=255)                                                                                                                                            ggCAGCAGCCTGCT‑TGACACCCTACCCTCGCAAGCAGGCGAATGGCTGCAATGGCACGGCGAAAGCCTGACCTACCTGTACTGCAACGACTTCAAGGTCAACCCTGATATCCAGCAGACGATCACCTTCGAGCACGCCAACAAATgtta         <  8:189662/149‑4 (MQ=255)                                                                                                                                              cagcagCCTGCT‑GGCCACCCTACCCTCGCACGCAGGCGAATGGCTGCAATGGCACGGCGAAAGCCTGACCTACCTGTACGGCAAGGACTTCAAGGTCAACCCTGATATCCAGCAGACCATCACCTTCGAGCACGCCAACAAATGGAAGt       <  5:134246/149‑1 (MQ=255)                                                                                                                                              cagcagCCTGCT‑GGCCACCCTACCCTCGCACGCAGGCGAATGGCTGCAATGGCACGGCGAAAGCCTGACCTACCTGTACGGCAAGGACTTCAAGGTCAACCCTGATATCCAGCAGACCATCACCTTCGAGCACGCCAACAAATGGAAGt       <  4:73142/149‑1 (MQ=255)                                                                                                                                                   gcCTGCT‑GGCCACCCTACCCTCGCACGCAGGCGAATGGCTGCAATGGCACGGCGAAAGCCTGACCTACCTGTACGGCAAGGACTTCAAGGTCAACCCTGATATCCAGCAGACCATCACCTTCGAGCACGCCAACAAATGGAAGTACGGc  >  3:248347/1‑149 (MQ=255)                                                                                                                                                   gcCTGCT‑GGCCACCCTACCCTCGCACGCAGGCGAATGGCTGCAATGGCACGGCGAAAGCCTGACCTACCTGTACGGCAAGGACTTCAAGGTCAACCCTGATATCCAGCAGACCATCACCTTCGAGCACGCCAACAAATGGAAGTACGGc  <  4:247631/149‑1 (MQ=255)                                                                                                                                                    ccTGCT‑GGCCACCCTACCCTCGCACGCAGGCGAATGGCTGCAATGGCACGGCGAAAGCCTGACCTACCTGTACGGCAAGGACTTCAAGGTCAACCCTGATATCCAGCAGACCATCACCTTCGAGCACGCCAACAAATGGAAGTACGGc  >  6:34770/1‑148 (MQ=255)                                                                                                                                                    |                                                                                                                                                     TTCACCTCTGCAGTAACCGCTACGCCGTCGACCCGCCGGAATGGGTCGCCCGGGGCGCCTGCGCGCCAAAAAAACGCCAACAAAAACAATAAGTCCGGGAAT‑CACAACATGAAGCGCATCACCACGTCCCTGTTGTTGGGCAGCAGCCTGCT‑GGCCACCCTACCCTCGCACGCAGGCGAATGGCTGCAATGGCACGGCGAAAGCCTGACCTACCTGTACGGCAAGGACTTCAAGGTCAACCCTGATATCCAGCAGACCATCACCTTCGAGCACGCCAACAAATGGAAGTACGGC  >  NC_002947/4882305‑4882598

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