Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,243,741	G→A	8.6%	intergenic (+30/+69)	tsaA → / ← rnt	peroxiredoxin/ribonuclease T

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,243,741	0	G	A	8.6%	266.4 / 10.2	94	intergenic (+30/+69)	tsaA/rnt	peroxiredoxin/ribonuclease T
Reads supporting (aligned to +/- strand):  ref base G (40/45);  new base A (4/4);  total (45/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.00e-01

GTTCCACGAAGAGCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGCTCAGCGGTTTTCTCG  >  NC_002947/1243594‑1243888                                                                                                                                                     |                                                                                                                                                       gTTCCACGAAGAGCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATAAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGt                                                                                                                                                       >  4:36473/1‑149 (MQ=255)  ttCCACGAAGAGCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTg                                                                                                                                                      <  7:258232/149‑1 (MQ=255)  ttCCACGAAGAGCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTg                                                                                                                                                      <  8:199152/149‑1 (MQ=255)    ccACGAAGAGCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGc                                                                                                                                                    >  4:29192/1‑149 (MQ=255)     cACGAAGAGCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGcc                                                                                                                                                   >  8:188570/1‑149 (MQ=255)     cACGAAGAGCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGcc                                                                                                                                                   <  7:393177/149‑1 (MQ=255)      aCGAAGAGCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAAACCGGCCTTCGTGGcc                                                                                                                                                   <  4:386770/149‑1 (MQ=255)      aCGAAGAGCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCg                                                                                                                                                  >  3:148525/1‑149 (MQ=255)      aCGAAGAGCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCg                                                                                                                                                  <  2:214189/149‑1 (MQ=255)          agagCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑ttt                                                                                                                                             >  3:298837/1‑149 (MQ=255)           gagCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCGGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCACGAAGGCCGG‑tttt                                                                                                                                            <  6:303874/149‑1 (MQ=255)            agCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑ttttt                                                                                                                                           >  5:231142/1‑149 (MQ=255)             gCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑tttttt                                                                                                                                          >  2:247845/1‑149 (MQ=255)             gCACGGCGAAGGGGGCCCTGCCAACTGGAAAACAGGCGACAAGGGCATGAACGCGTCGCCTGAAGGCGGTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑tttttt                                                                                                                                          <  3:34576/149‑1 (MQ=255)              cACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑tttttt                                                                                                                                          >  6:315553/1‑148 (MQ=255)              cACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑tttttt                                                                                                                                          >  1:195683/1‑148 (MQ=255)              cACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATA‑AAAAAACCGGCCACGAAGGCCGG‑ttttttc                                                                                                                                         <  1:419461/149‑2 (MQ=39)               aCGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCg                                                                                                                                                  >  4:49794/1‑140 (MQ=255)               aCGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCg                                                                                                                                                  <  3:49794/140‑1 (MQ=255)                cGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTAt                                                                                                                                       <  6:437085/149‑1 (MQ=255)                cGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTAt                                                                                                                                       <  1:2504/149‑1 (MQ=255)                  gCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAAAGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATg                                                                                                                                      >  7:206200/1‑148 (MQ=255)                      aGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTc                                                                                                                                 <  4:391456/149‑1 (MQ=255)                      aGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTc                                                                                                                                 <  1:99372/149‑1 (MQ=255)                        tgtgCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCt                                                                                                                               <  2:327585/149‑1 (MQ=255)                             ccTGCCAACTGGAAAAAAGGCGACAAGGCCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGAAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGcc                                                                                                                          <  5:339678/149‑1 (MQ=255)                              cTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑TAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCg                                                                                                                         <  7:257873/149‑1 (MQ=255)                                 ccAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGcc                                                                                                                      >  7:335645/1‑149 (MQ=255)                                   aaCTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTc                                                                                                                                 <  2:368022/136‑1 (MQ=255)                                   aaCTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTc                                                                                                                                 >  1:368022/1‑136 (MQ=255)                                   aaCTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTACATG‑AAAAAACCGGCCTTCCTGGACGG‑TTTTATTATGCCTTCCTGTGCCGTccac                                                                                                                    >  3:272033/1‑147 (MQ=255)                                         aaaaaaGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTcgcg                                                                                                              <  4:148525/149‑1 (MQ=255)                                          aaaaaGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCgg                                                                                                             <  8:58202/149‑1 (MQ=255)                                          aaaaaGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCgg                                                                                                             <  2:396788/149‑1 (MQ=255)                                                  gACAAGGGCATAAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑Gccc                                                                                                   <  3:36473/149‑1 (MQ=255)                                                   aCAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCcg                                                                                                  >  6:49445/1‑149 (MQ=255)                                                    cAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCcgc                                                                                                 <  1:310368/149‑1 (MQ=255)                                                        ggCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAAt                                                                                             <  1:422666/149‑1 (MQ=255)                                                           aTGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑ttttttt                                                                                                                                         >  4:413528/1‑104 (MQ=255)                                                           aTGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑ttttttt                                                                                                                                         <  3:413528/104‑1 (MQ=255)                                                           aTGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGcc                                                                                          >  7:305269/1‑149 (MQ=255)                                                             gAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCtt                                                                                        >  3:129118/1‑149 (MQ=255)                                                              aaCGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCtt                                                                                        <  6:231142/148‑1 (MQ=255)                                                               aCGCTTCGCCTGAAGGCGTTTCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCTTTa                                                                                      <  4:73810/149‑1 (MQ=255)                                                               aCGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCTTTa                                                                                      <  6:129443/149‑1 (MQ=255)                                                                cGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCTTTAg                                                                                     <  3:182885/149‑1 (MQ=255)                                                                      gCCTGAAGGCGTTGCTGCCTACCTGCGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCTTTAGCAGATc                                                                               <  1:203086/149‑1 (MQ=255)                                                                       ccTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCTTTAGCAGATc                                                                               >  3:181567/1‑148 (MQ=255)                                                                         tGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTATTTATGCCAACCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCTTTAGCAGATCAAt                                                                            >  1:356262/1‑149 (MQ=255)                                                                             ggCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCTTTAGCAGATCAATCGtt                                                                        >  4:167930/1‑149 (MQ=255)                                                                              gCGTTGCTGCCTACCTGAGCTAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCTTTAGCAGATCAATCGTTg                                                                       <  1:247845/149‑1 (MQ=255)                                                                               cGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCTTTAGCAGATCAATCGTTg                                                                       >  1:398707/1‑148 (MQ=255)                                                                                gTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAACCCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCTTTAGCAGATCAATCGTTGaa                                                                     <  8:289195/149‑1 (MQ=255)                                                                                gTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGGTTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCTTTAGCAGATCAATCGTTg                                                                       >  4:297052/1‑148 (MQ=255)                                                                                  tgctgcCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCTTTAGCAGATCAATCGTTGa                                                                      <  6:80531/146‑1 (MQ=255)                                                                                  tgctgcCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCTTTAGCAGATCAATCGTTGa                                                                      >  5:80531/1‑146 (MQ=255)                                                                                      gcCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATcgcg                                                               >  1:106326/1‑149 (MQ=255)                                                                                      gcCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATcgcg                                                               >  5:136908/1‑149 (MQ=255)                                                                                       ccTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATcgc                                                                <  2:195683/147‑1 (MQ=255)                                                                                          aCCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAg                                                           >  2:136947/1‑149 (MQ=255)                                                                                             tGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAgccg                                                        >  6:422786/1‑149 (MQ=255)                                                                                             tGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAgccg                                                        >  2:39711/1‑149 (MQ=255)                                                                                                      aCGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCAtttct                                               >  6:332688/1‑149 (MQ=255)                                                                                                          tGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCAtttctttc                                            >  5:343273/1‑148 (MQ=255)                                                                                                            gCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCa                                          <  5:422786/148‑1 (MQ=255)                                                                                                              aaGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCg                                       <  4:181567/149‑1 (MQ=255)                                                                                                                gCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATcgc                                                                >  8:404987/1‑122 (MQ=255)                                                                                                                gCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATcgc                                                                <  7:404987/122‑1 (MQ=255)                                                                                                                 cTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGtt                                    <  5:332688/149‑1 (MQ=255)                                                                                                                    tAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATg                                                                                            <  3:7682/90‑1 (MQ=255)                                                                                                                    tAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATg                                                                                            >  4:7682/1‑90 (MQ=255)                                                                                                                    tAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGa                                  <  5:49445/148‑1 (MQ=255)                                                                                                                          ctgcGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGcc                           <  4:298837/149‑1 (MQ=255)                                                                                                                          ctgcGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGcc                           <  1:47618/149‑1 (MQ=255)                                                                                                                             cGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGcag                          >  2:114695/1‑145 (MQ=255)                                                                                                                                gCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGaa                     >  8:385392/1‑149 (MQ=255)                                                                                                                                 cATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGaa                     <  4:129118/148‑1 (MQ=255)                                                                                                                                 cATA‑AAAAAACCGGCCACGAAGGCCGG‑TTTTTTCATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCa                                          >  1:270062/1‑127 (MQ=38)                                                                                                                                 cATA‑AAAAAACCGGCCACGAAGGCCGG‑TTTTTTCATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCa                                          <  2:270062/127‑1 (MQ=38)                                                                                                                                 cATA‑AAAAAACCGGCCACGAAGGCCGG‑TTTTTTCATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGaa                     >  8:324327/1‑148 (MQ=38)                                                                                                                                  aTG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACa                   <  7:385392/149‑1 (MQ=255)                                                                                                                                   ataaaaaaaCCGGCCACGAAGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGc                            >  1:341021/4‑140 (MQ=39)                                                                                                                                   ataaaaaaaCCGGCCACGAAGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGc                            <  2:341021/137‑1 (MQ=39)                                                                                                                                    g‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGCAGGCCCGGGAATGCCGTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGaa                     >  1:131320/1‑147 (MQ=255)                                                                                                                                    g‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCTGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGc                 >  3:310159/1‑149 (MQ=255)                                                                                                                                    g‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAgg                 >  5:133940/1‑148 (MQ=255)                                                                                                                                       aaaaaCCGGCCACGAAGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTATTTCCAGCGGTTGACGATg                             >  2:365120/1‑135 (MQ=255)                                                                                                                                         aaaCCGGCCTTCGTGGCCGGTTTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGCTCa              <  5:419776/149‑1 (MQ=255)                                                                                                                                         aaaCCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGCTCAg             >  6:437416/1‑149 (MQ=255)                                                                                                                                           cccGGCCTTCGTGGCCGG‑TTTTATTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGCTCAGCg           <  3:226170/148‑1 (MQ=255)                                                                                                                                              ggCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGCTCAGCGGtt        <  7:176691/149‑1 (MQ=255)                                                                                                                                                 cTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGCTCAGCGGTTTtc     >  2:266365/1‑149 (MQ=255)                                                                                                                                                    cGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCTCGGGCGC‑‑GCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAg                       >  7:302142/1‑128 (MQ=255)                                                                                                                                                    cGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCTCGGGCGC‑‑GCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAg                       <  8:302142/128‑1 (MQ=255)                                                                                                                                                     gTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGCTCAGCGGTTTTCTCg  <  8:335645/148‑1 (MQ=255)                                                                                                                                                     gTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGCTCAGCGGTTTTCTCg  <  6:343273/148‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                       GTTCCACGAAGAGCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATG‑AAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTCTTCGCGGGCGC‑‑GCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGCTCAGCGGTTTTCTCG  >  NC_002947/1243594‑1243888

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