Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,243,755	T→C	6.6%	intergenic (+44/+55)	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,755	0	T	C	6.6%	217.7 / 3.7	75	intergenic (+44/+55)	tsaA/rnt	peroxiredoxin/ribonuclease T
Reads supporting (aligned to +/- strand):  ref base T (35/35);  new base C (2/3);  total (37/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-01

CGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTG‑AAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGCTCAGCGGTTTTCTCGGTGTCATAA  >  NC_002947/1243609‑1243897                                                                                                                                                    |                                                                                                                                                  cGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTAt                                                                                                                                                 <  8:23049/149‑1 (MQ=255)  ggCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATg                                                                                                                                                <  6:130047/149‑1 (MQ=255)  ggCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATAAAAAAACCGGCCACGAAGGCCGG‑TTTTTTCATg                                                                                                                                                <  5:262059/149‑1 (MQ=38)  ggCGAAGTGTCCCCTGCCAACTCGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGAGCCCGG‑TTTTTTTATg                                                                                                                                                >  8:287032/1‑149 (MQ=255)          gtgCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCtg                                                                                                                                        >  4:209825/1‑149 (MQ=255)              ccTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGcc                                                                                                                                    >  8:41427/1‑149 (MQ=255)                  ccAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGc                                                                                                                                 >  3:243066/1‑148 (MQ=255)                       tGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTc                                                                                                                         >  4:24197/1‑149 (MQ=255)                               aGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGgcgc                                                                                                                 >  8:95374/1‑149 (MQ=255)                                   gACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGccc                                                                                                             >  7:162949/1‑149 (MQ=255)                                   gACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGccc                                                                                                             >  1:297134/1‑149 (MQ=255)                                     cAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCcgc                                                                                                           <  8:162949/149‑1 (MQ=255)                                     cAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCcgc                                                                                                           <  8:29526/149‑1 (MQ=255)                                        gggCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATAAAAAAACCGGCCACGAAGGCCGG‑TTTTTTCATGc                                                                                                                                               <  4:65354/112‑1 (MQ=37)                                        gggCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATAAAAAAACCGGCCACGAAGGCCGG‑TTTTTTCATGc                                                                                                                                               <  4:69525/112‑1 (MQ=37)                                        gggCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATAAAAAAACCGGCCACGAAGGCCGG‑TTTTTTCATGc                                                                                                                                               >  3:69525/1‑112 (MQ=37)                                        gggCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATAAAAAAACCGGCCACGAAGGCCGG‑TTTTTTCATGc                                                                                                                                               >  3:65354/1‑112 (MQ=37)                                        gggCAGGAAGGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAAGGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGaa                                                                                                        <  5:118849/149‑1 (MQ=255)                                         ggCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAAt                                                                                                       >  6:180530/1‑149 (MQ=255)                                         ggCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAAt                                                                                                       <  5:176677/149‑1 (MQ=255)                                           cATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCcgcg                                                                                                          <  6:192466/144‑1 (MQ=255)                                           cATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGaagat                                                                                                     >  5:192466/1‑146 (MQ=255)                                           cATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGc                                                                                                     >  1:266611/1‑149 (MQ=255)                                           cATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGc                                                                                                     <  3:147280/149‑1 (MQ=255)                                            aTGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGcc                                                                                                    >  5:281556/1‑149 (MQ=255)                                            aTGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGcc                                                                                                    >  2:193712/1‑149 (MQ=255)                                            aTGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGG‑GTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGcg                                                                                                    >  6:241336/1‑148 (MQ=255)                                              gAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCtt                                                                                                  >  6:290077/1‑149 (MQ=255)                                               aaCGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCtt                                                                                                  <  3:278678/148‑1 (MQ=255)                                                 cGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTAg                                                                                               <  7:41427/149‑1 (MQ=255)                                                    ttCGCCTGAAGGCGTTGCTGCCTACCTGAGCCAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTAg                                                                                               >  7:198636/1‑146 (MQ=255)                                                      cGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAg                                                                                            <  1:193712/147‑1 (MQ=255)                                                         cTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCaa                                                                                       >  6:93412/1‑149 (MQ=255)                                                              ggCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGtt                                                                                  <  4:37653/149‑1 (MQ=255)                                                               gCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGaagaggcc                                                                                                  <  1:114623/132‑7 (MQ=255)                                                               gCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGaagaggcc                                                                                                  >  2:114623/1‑126 (MQ=255)                                                                  ttGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTG‑aaa                                                                             <  8:116026/149‑1 (MQ=255)                                                                      tgcCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTG‑a                                                                               <  8:290336/143‑1 (MQ=255)                                                                      tgcCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTG‑a                                                                               >  7:290336/1‑143 (MQ=255)                                                                       gcCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTG‑AAATcgcg                                                                        >  4:169671/1‑149 (MQ=255)                                                                        ccTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTG‑AAATcgcg                                                                        <  3:209825/148‑1 (MQ=255)                                                                        ccTACCTGAGCGAGAACGAGGGCAAGCTGTAATTGGTGCGTGCATGAAAAAACCGGCCTTCGTGGCCTG‑TTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTAAGCAGATCAATCGTAG‑AAATcgcg                                                                        <  4:215720/148‑1 (MQ=255)                                                                             cTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTG‑AAATCGCGCCAgcc                                                                  >  8:146590/1‑149 (MQ=255)                                                                                  cGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTG‑AAATCGCGCCAGCCGCCCa                                                             >  2:272031/1‑149 (MQ=255)                                                                                   gagaACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTG‑AAATCGCGCCAGCCGCCCAt                                                            >  1:268121/1‑149 (MQ=255)                                                                                   gagaACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTG‑AAATCGCGCCAACCGCCCAt                                                            <  7:156972/149‑1 (MQ=255)                                                                                        cGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGgcgcgc                                                                                                               >  4:14930/1‑94 (MQ=255)                                                                                        cGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGgcgcgc                                                                                                               <  3:14930/94‑1 (MQ=255)                                                                                           tGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGAGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCCCGAATGCCTTTAGCACATCAATCGTTG‑AAAACGAGCCAGCCGCCCAatactttc                                                     >  3:248145/1‑148 (MQ=255)                                                                                                  cTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTG‑AAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGtt                                             <  5:93412/149‑1 (MQ=255)                                                                                                  cTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTG‑AAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGtt                                             <  7:146590/149‑1 (MQ=255)                                                                                                   tGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTG‑AAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGtt                                             >  2:108204/1‑148 (MQ=255)                                                                                                    gTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTAg                                                                                               <  2:491/98‑1 (MQ=255)                                                                                                    gTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTAg                                                                                               >  1:491/1‑98 (MQ=255)                                                                                                     tAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGa                                           <  1:53988/149‑1 (MQ=255)                                                                                                     tAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTG‑AAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGtt                                             <  1:272031/146‑1 (MQ=255)                                                                                                             gcGTGCATGAAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGAt                                                                                          >  3:57025/1‑94 (MQ=255)                                                                                                             gcGTGCATGAAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGAt                                                                                          <  4:57025/94‑1 (MQ=255)                                                                                                               gTGCATGAAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCTGGCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTG‑AAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCg                                   >  1:84051/1‑148 (MQ=255)                                                                                                                      aaaaaaCCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTG‑AAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGCt                         >  1:85213/1‑149 (MQ=255)                                                                                                                      aaaaaaCCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCCTTG‑AAATCGAGCCAGCCGCCCATTTCTTTACAGAGGTTGAAGATGCAGCAGAACAGc                          >  4:265999/1‑148 (MQ=255)                                                                                                                       aaaaaCCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTAGc                                                                                              >  8:14475/1‑80 (MQ=255)                                                                                                                       aaaaaCCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTAGc                                                                                              <  7:14475/80‑1 (MQ=255)                                                                                                                       aaaaaCCGGCCACGAAGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTG‑AAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGc                                     <  3:49053/136‑1 (MQ=255)                                                                                                                       aaaaaCCGGCCACGAAGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTG‑AAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGc                                     >  4:49053/1‑136 (MQ=255)                                                                                                                        aaaaCCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTG‑AAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGCTCa                       <  7:95374/149‑1 (MQ=255)                                                                                                                           aCCGTCCTTCGTGTCCGG‑TTTTTTTTTGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGAGAATGCCTTTAGGAGATCAATCGTTG‑AAATCGCGCCAGCCGTCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGCTCAggg                    <  8:198636/149‑3 (MQ=255)                                                                                                                              ggCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTG‑AAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGCTCAGCGGtt                 <  3:169671/149‑1 (MQ=255)                                                                                                                               gCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTG‑AAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGCTCAGCGGttt                >  4:109479/1‑149 (MQ=255)                                                                                                                                ccTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTG‑AAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGCTCAGCGGtttt               <  5:290077/149‑1 (MQ=255)                                                                                                                                    cGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTG‑AAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGCTCAGCGGTTTtctc            <  6:281556/148‑1 (MQ=255)                                                                                                                                     gTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTG‑AAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGCTCAGCGGTTTTCTCg           <  2:268121/148‑1 (MQ=255)                                                                                                                                     gTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CATCGCGGGCGCTCCCGCGAATGCCTTTAGCAGATCAATGGTTG‑AAATCGCGGCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGCTCAGCGGTTTTCTCg           <  4:243066/148‑1 (MQ=255)                                                                                                                                        gCCGGTTTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTTGCAGATGAATCGTTG‑AAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGCTCAGCGGTTTTCTCGgtg        >  4:121341/1‑149 (MQ=255)                                                                                                                                        gCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTG‑AAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGCTCAGCGGTTTTCTCGgtg        <  2:266611/148‑1 (MQ=255)                                                                                                                                              tttttttATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTG‑AAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGCTCAGCGGTTTTCTCGGTGTCATaa  <  7:145633/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                  CGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGG‑TTTTTTTATGCCTTCCTGTGCCGGCCT‑‑CTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTG‑AAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGCTCAGCGGTTTTCTCGGTGTCATAA  >  NC_002947/1243609‑1243897

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