Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	67,759	G→A	G193G (GGC→GGT)	araA ←	L‑arabinose isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	67,759	0	G	A	84.6%	106.9 / ‑0.2	41	G193G (GGC→GGT)	araA	L‑arabinose isomerase
Reads supporting (aligned to +/- strand):  ref base G (1/1);  major base A (30/3);  minor base C (0/6);  total (31/10)
Fisher's exact test for biased strand distribution p-value = 2.57e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTTGCCGTGGATTTGTGTG‑GCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCGCCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAATCGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAGACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACCACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATAC  >  NC_000913/67528‑67994                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                            aCCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTT‑GCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTGACCGCCACtt                                                                                                                                                                                                                          >  1:198436/1‑251 (MQ=255)  ccTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTT‑GCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTGACCGCCACTTc                                                                                                                                                                                                                         >  1:179127/1‑251 (MQ=255)             gAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTT‑GCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTGACCGCCACTTCACGCATGTTAt                                                                                                                                                                                                              >  1:18990/1‑251 (MQ=255)             gAAACGCTTCATCCACCGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTT‑GCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTGACCGCCACTTCACGCATGTTAt                                                                                                                                                                                                              >  1:112996/1‑251 (MQ=255)                  gCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTT‑GCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTGACCGCCACTTCACGCATGTTATCGcc                                                                                                                                                                                                          >  2:81523/1‑250 (MQ=255)                      cATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTT‑GCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTGACCGCCACTTCACGCATGTTATCGCCaaaac                                                                                                                                                                                                     >  2:132372/1‑249 (MQ=255)                         ccccAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTT‑GCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGc                                                                                                                                                                                                  >  1:161347/1‑251 (MQ=255)                                      cgcgCCGCTTCCAGCACGTTCTGTCGTTTTTCGCCGTGGATTTGTGTT‑GCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATg                                                                                                                                                                                     >  1:375494/1‑251 (MQ=255)                                                 cggccctttCTGTCGTTTTTCGCCGTGG‑TTTTTGTTTGCAGGCGTCATGGTGTAGCATCTTTCGTACTCATCGACCAGCCCGTAAACATCGCCGTCGCTGATGGAGTTCCCCCCCTGCCCAAGTTCCCCACCCGCCCGGGTATTGCCGGAAAACCCGAACTTGATCTGTGCGGCAACTTTTTCCCCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCt                                                                                                                                                                          <  1:56088/244‑1 (MQ=255)                                                     aCGTTCTGTCGTTTTTCGCCGTGGATTTGTGTT‑GCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGttt                                                                                                                                                                      >  1:123889/1‑251 (MQ=255)                                                         tCTGTGGTTTTTCGCCGGGGGCTTGTGTT‑GCAGACGTCATGGTGTAGCGGCTTTCGTACTCATCGACCAGCGCGTTAACACCGCCGTCGCTGGTGGAGTTCCCCACCTGCACCAGATCGCCAACCGCCCGGGTATTGACGGAGAACCCGCACTTGATCTGTCCGGCAACTTTATCACCATCGGTGACCGCCACTTCACGCATGTTATCGCCACAACGGCATACTTTCAGATGACGGGTATCCTGTTTaga                                                                                                                                                                   >  2:16595/1‑250 (MQ=255)                                                              cgattgtcgCCGCGGATTTGTGTT‑GCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCAAACCGCCCAGGCATTGCCGGAGAAACCGAACTTGATCTGCGCGGCAACTTTATCACCATCGGGGACCGCCACTTCACGCATGTTATCGCCAAAACGGCCGACTTTCCGATGACGGGTCTCCTGTTTCGAAACCGc                                                                                                                                                             >  2:331528/9‑251 (MQ=255)                                                                         gTGGATTTGTGTG‑GCAGGCGTCTGGGTGTCGCAGTTTTCGTCCTCATCCACCAGCGCTTTAACACCGCCGGCGCTGATGGAGTCCCCCCCCTCCACCAGATCGCCAACCGCCCAGGTATTGCCGGAAAAACCGAACTTGCTCTGTGCGCCACCTTTACCACCATCGTTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATcc                                                                                                                                                  <  2:259673/251‑1 (MQ=255)                                                                          tGGATTTGTGTT‑GCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATcc                                                                                                                                                  >  1:325181/1‑250 (MQ=255)                                                                                 tgggttgCAGGCTTCATGGTGTAGCAGCTTTCGTCCCCATCGCCCCGCGCTTTACCATCGCGTTCTGGGAGGGAGTCCACCACCGGCTCCAAAGGCCCACCCCCCCAGGTATGGACGGAGAAACCGATTTTGATCTGGGCGGCAACTTTACCACCATCGGGGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGc                                                                                                                                            <  1:238432/244‑1 (MQ=255)                                                                                    gttgCAGGCGTCGTGGTGTAGCATCTTCCGTACCCATCGCCCAGCGCTTTAACATCCCCTTCGCTGATGGATTTCCCCACCGGCCCCAGATCGCCAACCGCCCAGGTATGGAGGGAAAAACCGAACTTGTTCTGGGCGGCAACTTTATCCCCACCGGGGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGa                                                                                                                                         <  1:220480/247‑1 (MQ=255)                                                                                       gCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACg                                                                                                                                     >  1:222815/1‑251 (MQ=255)                                                                                             gTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATg                                                                                                                               >  2:79140/1‑251 (MQ=255)                                                                                                  ggTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCtt                                                                                                                          >  1:44535/1‑251 (MQ=255)                                                                                                  ggTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCtt                                                                                                                          >  1:175584/1‑251 (MQ=255)                                                                                                     gtAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCGCCATCGGTGACCGCCACTTCACGCATGTTATCCCCAAATCGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAGACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGtt                                                                                                                       >  2:14564/1‑251 (MQ=255)                                                                                                           gcTTTCGTAATCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCt                                                                                                                 >  1:249840/1‑251 (MQ=255)                                                                                                              ttCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGcc                                                                                                              >  1:224380/1‑251 (MQ=255)                                                                                                                 gTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGCCGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGt                                                                                                           >  2:244009/1‑251 (MQ=255)                                                                                                                                                   cTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCTTAcg                                                                         >  2:295283/1‑251 (MQ=255)                                                                                                                                                   cTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATAcg                                                                         >  1:293938/1‑251 (MQ=255)                                                                                                                                                            ttCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAAt                                                                >  2:335058/1‑251 (MQ=255)                                                                                                                                                                ccaccTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGa                                                              >  1:98392/1‑249 (MQ=255)                                                                                                                                                                        aCCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTg                                                                                                          >  2:195056/1‑197 (MQ=255)                                                                                                                                                                        aCCAGACCGCCACCCGCCCAGGTATTGACGGAGAAACCGATCTTGATCTGGGCGGCACTTTTATCCCCAGGGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTg                                                                                                          <  1:195056/197‑1 (MQ=255)                                                                                                                                                                                          cAGGTATTGACGGAAAAACCGAACTTAATCTGTGCGACAATTTTACCCCCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGt                                  <  1:228358/251‑1 (MQ=255)                                                                                                                                                                                             gTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTg                               >  2:92483/1‑251 (MQ=255)                                                                                                                                                                                                 tGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTggttc                           >  1:196914/1‑251 (MQ=255)                                                                                                                                                                                                  gagggaaaaaCCGACTTTGAGCTGGGCGGCAACTTTACCACCACCGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTggttca                          <  1:270628/244‑1 (MQ=255)                                                                                                                                                                                                      gagaAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGCGCCCGACCCTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGAAACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTggttcaggtt                      >  1:367616/1‑251 (MQ=255)                                                                                                                                                                                                       agaAACCGAACTTGATCTGTACGGCAACTTTATCACCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTggttcaggttc                     >  2:366733/1‑251 (MQ=255)                                                                                                                                                                                                         aaaCCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCAt                   >  1:210499/1‑251 (MQ=255)                                                                                                                                                                                                                   tGATCTGGGCGGCAACTTTATCCCCACGGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCat          <  1:3728/250‑1 (MQ=255)                                                                                                                                                                                                                    gATCTGTGCGGCAACTTTACCCCCAGCGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCata         <  1:130338/250‑1 (MQ=255)                                                                                                                                                                                                                         gtgCGGCATCTTAACCGCCACCGGTGACCGCCTCTTCACGCATGTTATCCCCAAAGCGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAGACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACCACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCg      <  1:14564/248‑1 (MQ=255)                                                                                                                                                                                                                          tgCGGCAACTTTATCACCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGATATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATAc  >  2:185572/1‑251 (MQ=255)                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                            ACCTTGTTCCAGGAAACGCTTCATCCCCAGCTCAATACGCGCCGCTTCCAGCACGTTCTGTCGTTTTTTGCCGTGGATTTGTGTG‑GCAGGCGTCATGGTGTAGCAGCTTTCGTACTCATCGACCAGCGCGTTAACATCGCCGTCGCTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCGCCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAATCGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAGACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACCACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATAC  >  NC_000913/67528‑67994

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑