Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,566,212	C→G	9.2%	G124R (GGT→CGT)	glgP ←	glycogen phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,566,212	0	C	G	9.2%	227.0 / 11.4	75	G124R (GGT→CGT)	glgP	glycogen phosphorylase
Reads supporting (aligned to +/- strand):  ref base C (30/38);  new base G (4/3);  total (34/41)
Fisher's exact test for biased strand distribution p-value = 6.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.23e-01

GCCAGTAGTCTGGCGACTCTTTCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGT‑CCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAATCAAAAACTCCATCGAC  >  NC_000913/3566065‑3566357                                                                                                                                                     |                                                                                                                                                  gccagTAGTCTGGCGACTCTTTCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGT‑CCCAGGCCACCg                                                                                                                                                  >  2:1027857/1‑149 (MQ=255)   cagTAGTCTGGCGCCTCTTTCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTGCCCGGGCCCCCgg                                                                                                                                                 >  1:325978/1‑148 (MQ=255)        gtCTGGCGACTCTTTCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGT‑CCCAGGCCACCGGTACCgc                                                                                                                                           >  1:1002134/1‑148 (MQ=255)           tGGCGACTCTTTCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCATAGAATCAAGGAAGCAAGCCGCCAGACGT‑CCCAGGCCACCGTTACCGAGGc                                                                                                                                        <  1:933692/149‑1 (MQ=255)           tGGCGACTCTTTCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGT‑CACAGGCCACCGTTACCGAGGc                                                                                                                                        >  1:1164038/1‑149 (MQ=255)           gggCGACTCTTTCTGGCTACCGTTACCGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGT‑CCCAGGCCACCGTTACCGAGGc                                                                                                                                        <  1:929248/148‑1 (MQ=255)             gCGACTCTTTCTGGCTACCGTTAACGATGTGCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGTAAGCAAGCCGCCAGACGT‑CCCAGGCCTCGGTAACGGTGGCCt                                                                                                                                      <  2:791534/149‑1 (MQ=255)              cGACTCTTTCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGT‑CCCAGGCCTCGGTAACGGTGGCCTg                                                                                                                                     <  2:391807/149‑1 (MQ=255)              cGACTCTTTCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGT‑CCCAGGCCCCCGTTACCGCGGccgg                                                                                                                                     >  2:1348808/1‑147 (MQ=255)               gACTCTTTCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGT‑CCCAGGCCTCGGTAACGGTGGCCTgg                                                                                                                                    <  1:609032/149‑1 (MQ=255)                  tctTTCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCTAGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGT‑GCCCGCCCCCCCTTTCCCAGGCCGGGGtg                                                                                                                                 >  2:1117210/1‑148 (MQ=255)                  tctTTCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGT‑CCCCCGCCCCCCGGGCCGAGGCCTggggg                                                                                                                                 >  2:1111749/1‑147 (MQ=255)                  tctTTCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGT‑CCCAGGCCACCGTTACCGAGGCCTGGGTc                                                                                                                                 <  1:133148/149‑1 (MQ=255)                      tCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGT‑CCCAGGCCACCGTTACCGAGGCCTGGGTCAttt                                                                                                                             <  1:388342/149‑1 (MQ=255)                        tGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGT‑CCCAGGCCACCCTTAACGAGGCCTGGGTCAatttt                                                                                                                           >  2:1162766/1‑148 (MQ=255)                            aaCCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGT‑CCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTtca                                                                                                                       <  1:137449/148‑1 (MQ=255)                             aCCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGCCGT‑CCCAGGCCCCCCGTCACGAGGCCGGGGTGAttttagtcct                                                                                                                      >  2:153983/1‑143 (MQ=255)                               cGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGT‑CCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTtcatca                                                                                                                    <  1:1245448/149‑1 (MQ=255)                               cGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGT‑CCCAGGCCACCGGTACCGAGGCCTGGGTCATTTACTtcatca                                                                                                                    >  2:958678/1‑149 (MQ=255)                                 ttAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGT‑CCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCatcaaa                                                                                                                  >  1:814271/1‑148 (MQ=255)                                           tCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGT‑CCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCg                                                                                                        <  1:66727/149‑1 (MQ=255)                                             tGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGT‑CCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCgag                                                                                                      <  1:777058/149‑1 (MQ=255)                                             tGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGT‑CCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCgag                                                                                                      <  1:1024453/149‑1 (MQ=255)                                                 tGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGT‑CCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGAttt                                                                                                  <  1:792204/149‑1 (MQ=255)                                                      aTACCGTAGTCATAGCGTAGGCCGAAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGT‑CCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAAccc                                                                                             <  2:896479/149‑1 (MQ=255)                                                      aTACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGCCGT‑CCCACGCCCCCGGTACCGAGGCCGGGGTCAGTTTAGTAAGAAATAAGCACTTCGAGATTTAAccc                                                                                             >  2:120897/1‑149 (MQ=255)                                                         ccGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGCGAATCAAGGAAGCAAGCCGCCAGACGT‑CCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCAt                                                                                          <  2:71903/149‑1 (MQ=255)                                                          cGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGT‑CCCCGGCCCCCGGTACCGCGGCCGGGGGCATTTTATTAATCAATCAGCTCGTCGAGATTTAACCCCaac                                                                                         >  1:566418/1‑147 (MQ=255)                                                          cGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGT‑CCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATc                                                                                         >  2:1423978/1‑149 (MQ=255)                                                           gTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGT‑CCCCGGCCCCCGTTACCGAGGCCGGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCg                                                                                        >  2:447120/1‑149 (MQ=255)                                                             aGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGT‑CCCAGGCCCCCGGTACCGAGGCCGGGGTCATTTTGTTCATCAATCAGCGCGTCGAGATTTAACCCCCTCGcg                                                                                      >  1:1344658/1‑148 (MQ=255)                                                               tCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGT‑CCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTc                                                                                    <  1:688340/149‑1 (MQ=255)                                                               tCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGT‑CCCAGGCCACCGTTACCGAGGCCGGGGTCATTTTGTTCATCAATCAGCTCGTCGGGTTTTAACCCCCTCGCGTc                                                                                    >  2:250585/1‑149 (MQ=255)                                                               tCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCAAGACGT‑CCCCGGCCACCGTTTCCGCGGCCTGGGGGAGGTTCTTGATAAAACCGCTCTTCGAGATTTAAACCCAACGCTTc                                                                                    >  2:833704/1‑149 (MQ=255)                                                                 aTAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGT‑CCCAGGCCCCCGTTACCGAGGCCGGGGTGATTTTATTCATCAATCAGCTCTTCGAGGTTTAACCCCAACGCTTCCa                                                                                  >  2:973239/1‑149 (MQ=255)                                                                 aTAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGT‑CCCAGGCCACCGTTTCCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCa                                                                                  <  2:1362917/149‑1 (MQ=255)                                                                 aTAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGT‑CCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCa                                                                                  <  2:1266063/149‑1 (MQ=255)                                                                   aGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGT‑CCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGt                                                                                >  1:889885/1‑149 (MQ=255)                                                                   aGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGT‑CCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGt                                                                                <  2:475553/149‑1 (MQ=255)                                                                     cGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGT‑CCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTgc                                                                              <  2:889885/149‑1 (MQ=255)                                                                     cGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGT‑CCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTAGCGATTTAACCCCATCGCTTCCAGTgg                                                                              >  1:888053/1‑148 (MQ=255)                                                                     cGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGT‑CCCAGGCCACCGTTACCGAGGCCGGGGTCATTTTCTTAATCAATCAGCTCTTCGAGATTTAACCCCATAGCTTCCAgggg                                                                              >  2:1367016/1‑146 (MQ=255)                                                                            cGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGT‑CCCAGGCCACCGGGACCGCGGCGGGGGTCAATTTGTTCATCAATCACCGCTTCGAGGTTTAACCCCAGCGCTTCCAGAGCGccgggg                                                                       >  2:1160442/1‑145 (MQ=255)                                                                            cGTAACCGCGCCCCGGCAACCCTAACGTCGCAAGAGAATCAAGGAATCAAGCCGCCAGACGT‑CCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGt                                                                       <  2:987047/149‑1 (MQ=255)                                                                              tAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGT‑CCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTAc                                                                     <  2:1311683/149‑1 (MQ=255)                                                                               aaCCGCGCCCCGGCAACCCTAACGTCGCCAGAGAAGCAAGGAAGCAAGCCGCCAGACGT‑CCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCACTCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACa                                                                    <  2:1128756/149‑1 (MQ=255)                                                                                 ccGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGT‑CCCACGCCCCCCTTACCGACGCCAGGGTCATTTTATTAATAAAACAGCGCTTAGAGAATTAAACCCCACGCTTCCCGTGCGGCCGGTAccac                                                                  >  2:1000933/1‑146 (MQ=255)                                                                                   gcgcCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCATACGT‑CCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCtt                                                                >  1:391061/1‑149 (MQ=255)                                                                                       ccgtgCAACCCTAACGTCGCCAGAGAATCTAGGAAGCATGCAGCCAGACGT‑CCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCTATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACTTCTTCGTa                                                            <  2:1034805/145‑1 (MQ=255)                                                                                          ggCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGT‑CACAGGCCTCCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAAt                                                         <  2:1164038/149‑1 (MQ=255)                                                                                                ccTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGT‑CCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAg                                                   <  1:1083589/149‑1 (MQ=255)                                                                                                   aaCGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGT‑CCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGa                                                <  2:1185539/149‑1 (MQ=255)                                                                                                    aCGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGT‑CCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGac                                               <  1:1162766/149‑1 (MQ=255)                                                                                                        cGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGT‑CCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGacaaca                                           <  1:1423978/149‑1 (MQ=255)                                                                                                        cGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGT‑CCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGacaaca                                           <  2:566418/149‑1 (MQ=255)                                                                                                            agagAATCAAGGAAGCAAGCCGCCAGACGT‑CCCAGGCCACCCTTACCGCCGGCGGGGTGATTTTATTCATCAATAAGCTCTTAGAGATTTTACCCCCTCGCTTACCGTGCGCCCGGTACAACTTCGTAAATTACTAGCGccaaaaaggg                                       >  2:706657/1‑143 (MQ=255)                                                                                                                    aaGGAAGAAAGCCGCCAGACGT‑CCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGgaga                               <  2:1354775/149‑1 (MQ=255)                                                                                                                          gcaagcCGCCAGACGT‑CCCACGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTAc                         <  2:263770/149‑1 (MQ=255)                                                                                                                             agcCGCCAGACGT‑CCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGc                      <  1:1111749/149‑1 (MQ=255)                                                                                                                             agcCGCCAGACGT‑CCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGAGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGc                      <  1:1117210/149‑1 (MQ=255)                                                                                                                              gccgccAGACGT‑CCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGcc                     <  1:824030/149‑1 (MQ=255)                                                                                                                                cgccAGACGT‑CCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCaa                   <  1:924/149‑1 (MQ=255)                                                                                                                                  ccAGACGT‑CCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAATc                 <  1:48088/149‑1 (MQ=255)                                                                                                                                  ccAGACGT‑CCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAATc                 <  1:1348808/149‑1 (MQ=255)                                                                                                                                       cGT‑CCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTCGCGACAACATGGCGTTGGAGAGCGTACGGCCAATCaaaaa            <  1:425549/149‑1 (MQ=255)                                                                                                                                        tgacccAGGCCTCGGTAACGGTGGCCTGGGACATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTAGTAAATTGCTAGCGACAACATGGCGTTGGATAGCGTACGGCCAATCaaaaa            >  1:1399815/4‑149 (MQ=255)                                                                                                                                        tgacccAGGCCTCGGTAACGGTGGCCTGGGACATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCACATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTACTAGCGACAACATGGCGTTGGAGCGCGTACGGCCAATCaaaaa            >  1:543735/4‑149 (MQ=255)                                                                                                                                           cccAGGCCTCGGTAACGGTGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAATCAAAAACTc         >  1:859567/1‑149 (MQ=255)                                                                                                                                           cccAGGCCTCGGTAACGGTGGCCTGGGTCATTTTCTTCATCAATCAGCTATTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGGACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAATCAAAAACTc         >  2:566740/1‑149 (MQ=255)                                                                                                                                              aGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATGTTCGGAAATTCCTAGCGACAACATGGCGTGGGAGAGCGTACGGCCAATCAAAAACTCCAt      >  2:1229983/1‑149 (MQ=255)                                                                                                                                              aGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAATCAAAAACTCCAt      >  2:847770/1‑149 (MQ=255)                                                                                                                                               ggCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAATCAAAAACTCCATc     <  2:1344658/149‑1 (MQ=255)                                                                                                                                                gCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAATCAAAAACTCCATCg    >  2:1193184/1‑149 (MQ=255)                                                                                                                                                  cACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAATCAAAAACTCCATCGAc  <  1:1027857/149‑1 (MQ=255)                                                                                                                                                  cACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGATGGAGAGCGTACGGCCAATCAAAAACTCCATCGAc  >  2:633217/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                  GCCAGTAGTCTGGCGACTCTTTCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGT‑CCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAATCAAAAACTCCATCGAC  >  NC_000913/3566065‑3566357

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