Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,882,444	G→C	6.0%	G11A (GGC→GCC)	PP_4291 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,882,444	0	G	C	6.0%	213.8 / 3.3	67	G11A (GGC→GCC)	PP_4291	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (33/30);  new base C (2/2);  total (35/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.53e-01

TGCCCCTCGTTCACCTCTGCAGTAACCGCTACGCCGTCGACCCGCCGGAATGGGTCGCCCGGGGCGCCTGCGCGCCAAAAAAACGCCAACAAAAACAATAAGTCCGGGAAT‑CACAACATGAAGCGCATCACCACGTCCCTGTTGTTGGGCAGCAGCCTGCT‑GGCCACCCTACCCTCGCACGCAGGCGAATGGCTGCAATGGCACGGCGAAAGCCTGACCTACCTGTACGGCAAGGACTTCAAGGTCAACCCTGATATCCAGCAGACCATCACCTTCGAGCACGCCAACAAATGGAA  >  NC_002947/4882296‑4882591                                                                                                                                                      |                                                                                                                                                     tGCCCCTCGTTCACCTCTGCAGTAACCGCTACGCCGTCGACCCGCCGGAATGGGTCGCCCGGGGCGCCTGCGCGCCAAAAAAACGCCAACAAAAACAATAAGTCCGGGAAT‑CACAACATGAAGCGCATCACCACGTCCCTGTTGTgggg                                                                                                                                                      >  1:279701/1‑149 (MQ=255)      cTCGTTCACCTCTGCAGTAACCGCTACGCCGTCGACCCGCCGGAATGGGTCGCCCGGGGCGCCTGCGCGCCAAAAAAACGCCAACAAAAACAATAAGTCCGGGAAT‑CACAACATGAAGCGCATCACCACGTCCCTGTTGTTGGgcagca                                                                                                                                                 <  3:149364/149‑1 (MQ=255)          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:235216/1‑149 (MQ=255)                              tACGCCGTCGACCCGCCGGAATGGGTCGCCCGGGGCGCCTGCGCGCCAAAAAAACGCCAACAAAAACAATAAGTCCGGGAAT‑CACAACATGAAGCGCATCACCACGTCCCTGTTGTTGGGCAGCAGCCTGCT‑GGCCACCCTACCCTcgc                                                                                                                        >  5:290957/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                                                                     >  8:56183/1‑119 (MQ=255)                                                                                                               t‑CACAACATGAAGCGCATCACCACGTCCCTGTTGTTGGGCAGCAGCCTGCT‑GGCCACCCTACCCTCGCACGCAGGCGAATGGCTGCAATGGCACGGCGAAAGCCTGACCTACCTGTACg                                                                     <  7:56183/119‑1 (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)                                                                                                                                                      |                                                                                                                                                     TGCCCCTCGTTCACCTCTGCAGTAACCGCTACGCCGTCGACCCGCCGGAATGGGTCGCCCGGGGCGCCTGCGCGCCAAAAAAACGCCAACAAAAACAATAAGTCCGGGAAT‑CACAACATGAAGCGCATCACCACGTCCCTGTTGTTGGGCAGCAGCCTGCT‑GGCCACCCTACCCTCGCACGCAGGCGAATGGCTGCAATGGCACGGCGAAAGCCTGACCTACCTGTACGGCAAGGACTTCAAGGTCAACCCTGATATCCAGCAGACCATCACCTTCGAGCACGCCAACAAATGGAA  >  NC_002947/4882296‑4882591

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