Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,283,465	G→A	P480S (CCG→TCG)	actP ←	acetate/glycolate:cation symporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,283,465	0	G	A	90.2%	103.2 / ‑1.3	44	P480S (CCG→TCG)	actP	acetate/glycolate:cation symporter
Reads supporting (aligned to +/- strand):  ref base G (3/1);  new base A (29/8);  total (34/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.43e-01

GTCAGACGACCGGGGAAGGGGAGAGATTAATGCGCGCGGCCTTGCTCAACGCCAAAGCCGGTCTGGGAGCGGATAAACTGCGCGCGGAACAGTTCACGCTCACGCGCGCCTTCCGCTGAGTTATCGGTTGCCGAGAAGAACCAGATGCCGAGGAATGCCACGGTGATAGAGAACAGCGCCGGGTATTCATACGGGAAGATGGCTTTTTCGTGACCAAGGATCTGTACC‑CAAATCGTCGGGCCGAGGATCATCAGTACTACTGCGGTAATCAGCCCCAGCCAGCCACCCATCATCGCGCCACGCGTGGTCAGTTTCGACCAGTACATAGAAAGCAGAATGATCGGGAAGTTACAGCTCGCCGCGATGGCAAACGCCAGCCCCACCATAAAGGCGATGTTCTGATTCTCA‑AACAGCACGCCGAGGATAATCGCAATCACGCCGAGGATCAGTACGGTGATTTTTGATACCCGCAGCTCTTCACGT  >  NC_000913/4283227‑4283709                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                      gTCAGACGACCGGGGAAGGGGAGAGATTAATGCGCGCGGCCTTGCTCAACGCCAAAGCCGGTCTGGGAGCGGATAAACTGCGCGCGGAACAGTTCACGCTCACGCGCGCCTTCCGCTGAGTTATCGGTTGCCGAGAAGAACCAGATGCCGAGGAATGCCACAGTGATAGAGAACAGCGCCGGGTATTCATACGGGAAGATGGCTTTTTCGTGACCAAGGATCTGTACC‑CAAATCGTCGAGCCGAGGatcat                                                                                                                                                                                                                                           >  1:605251/1‑251 (MQ=255)   cAGACGACCGGGGAAGGGGAGAGATTAATGCTCGCGGACTTGCTCAACGCCAAAGCCGGTCTGGGAGCGGATAAACTGCGCGCGGCACAGTTCACGCTCACGCGCGCCTTCTGCTGAGTTATCGGATACCGTAAAGAACCAGATGCCGAGGAGTGCCCCGGTGATAGAGAAAAGCGCCGGGAATTCATACGGGAAGATGGCTTTTGCGTGCCCAAGGATCTGTAAC‑CAAATCGTCGAGcctggaatcaac                                                                                                                                                                                                                                          >  2:459069/1‑240 (MQ=255)           cGGGGAAGGGGAGAGATTAATGCGCGCGGCCTTGCTCAACGCCAAAGCCGGTCTGGGAGCGGATAAACTGCGCGCGGAACAGTTCACGCTCACGCGCGCCTTCCGCTGAGTTATCGGTTGCCGCGAAGAACCAGATGCCGAGGAATGCCACGGTGATAGAGAACAGCGCCGGGTATTCATACGGGAAGATGGCTTTTTCGTGACAAAGGATCTGTACC‑CAAATCGTCGCGCCGGGGttcatccatcctact                                                                                                                                                                                                                                 >  2:193492/1‑242 (MQ=255)                         gaTTAATGCGCGCGGCCTTGCTCAACGCCAAAGCCTGTCTGGGAGCGGATAAACTGCGCGCGGAACAGTTCACGCTCACGCGCGCCTTCCGCTGAGTTATCGGTTGCCGAGAAGAACCCGATGCCGCGGAATGCCCAGGTGATAGAGAACAGCGCCGGGTATTCATACGGAAAGATGGCTTTTTCGTGACCAAGGATCTGTACC‑CAACTCGTCGGGCCGAGAATCATCAGTCCTACTCCGGTAATCAGccc                                                                                                                                                                                                                   >  2:89045/1‑251 (MQ=255)                         gaTTAATGCGCGCGGCCTTGCTCAACGCCAAAGCCGGTCTGGGAGCGGATAAACTGCGCGCGGAACAGTTCACGCTCACGCGCGCATTCCGCTGAGTTATCGGTTGCCGAGCAGAACCAGATGCCGCGGACTGCCACGGTGATAGAGAACAGCGCCGGGTATTCATACGGGAAGATGGCTTTTTCGTGACCAAGGATCTGTACC‑CAAATCGTCGAGCCGAGGATCATCAGTACTACTGCGGTAATCAGccc                                                                                                                                                                                                                   >  1:83753/1‑251 (MQ=255)                              atgagcgcgGCCTTGCTCAACGCCAAGGCCGGTCTGGGGGCGGATAAACTGCGCGCGGAACAGTTCACGCTCACGCGCGCCTTCCGCTGAGTTATCGGTTGCCGAGAAGAACCAGATGCCGAGGAATGCCACGGTGATAGAGAACAGCGCCGGGTATTCATACGGGCAGATGGCTTTTTCGTGACCAAGGATCTG‑ACCACAAATCGTCGGGCCGAGGATCATCAGTACTACTGCGGACTTCAGCCccagcc                                                                                                                                                                                                              >  2:105643/5‑251 (MQ=255)                                                       aaGCCGGTCTGGGATCGGTTAAACTGCGCGCGGAACAGTTCACGCTCACGCGCGCCTTCCGCTGAGTTATCGGTTGCCGAGAAGAACCAGATGCCGAGGTATGCCACGGTGATAGAGCACAGCGCCGGGTCTTCATACGGGAAGATGGCTTTTTCGTGTCCAAGGATCTGTACC‑CTAATCATCGAGCCGAGGTTCATCTGTACTCCTTCGTTAATCAGCCCCAGCCCGCCACCCATCATCGCGCCCCGCGt                                                                                                                                                                                     >  2:503811/1‑251 (MQ=255)                                                       aaGCCGGTCTGGGAGCGGATAAACTGCGCGCGGAACAGTTCACGCTCACGCGCGCCTTCCGCTGAGTTATCGGTTGCCGAGAAGAACCAGATGCCGAGGAATGCCACGGTGATAGAGAACAGCGCCGGGTATTCATACGTGAAGATGGCTTTTTCGTGACCAAGGATCTGTACC‑CAAATCGTCGAGCCGAGGATCATCAGTACTACTGCGGTAATCAGCCCCAGCCAGCCACCCATCATCGCGCCACGCGt                                                                                                                                                                                     <  1:83583/251‑1 (MQ=255)                                                              tCTGGGAGCGGATAAACTGCGCGCGGAACAGTTCACGCTCAAGCGCGCCTTCCGCTGAGTTATCGGATGCCGAGAAGAACCAGATGCCGAGGAATGCCACGGTGATAGAGAACAGCGCCGGGTATTCATACGGGAAGATGGATTTTTAGTGACCACGGATCTGTACC‑CAAATCGTCGAGCCGAGGATCATCAGTACCACTGCGGTAATCAGCCCCAGCCAGCCACCCATCATCGCGCCAAGCGGGGTCAGt                                                                                                                                                                              >  2:113582/1‑251 (MQ=255)                                                                                   cgcgGAACAGTTCACGCTCACGCGCGCCTTCCGCTGCGTTATCGGTTGCCGAGAAGAACCAGATGCCGAGGAATGCCACGGTGATAGAGAACAGCGCCGGGTATTCATACGGGAAGATGGCTTTTTCGTGACCAAGGATCTGTACC‑CAAATCGTCGAGCCGAGGATCATCAGTACTACTGCGGTAATCAGCCCCAGCCAGCCACCCATCATCGCGCCACGCGTGGTCAGTTTCGACCAGTACATAGAAAGc                                                                                                                                                         >  1:611297/1‑251 (MQ=255)                                                                                    gcgGAACAGTTCACGCTCACGCGCGCCTTCCGCTGAGTTATCGGTTGCCGAGATGTACCCGATTCCGAGGACTGCCACGGTGCTAGAGACCAGCGCCGGGTATTCCTACGGGCAGCTGGCTTTTTCGTGACCAAGAATCTGTACC‑CACATCGTCGAGCCGAGGCTCATCAGTACTCCTGCGCTAATCAGCCCCAGCCAGCCACCCATCATCGCGCCACGCGTGGTCTGTTTCGACCAGTACAtcgcacccc                                                                                                                                                        >  1:448803/1‑243 (MQ=255)                                                                                        aaCAGTTCACGCTCACGCGCGCCTTCCGCTGAGTTATCGGTTGCCGAGAAGAACCAGATGCCGAGGAATGCCACGGTGATAGAGAACAGCGCCGGGTATTCATACGGGAAGATGGCTTTTTCGTGACCAAGGATCTGTACC‑CAAATCGTCGAGCCGAGGATCATCAGTACTACTGCGGTAATCAGCCCCAGCCAGCACCCCATCATCGCGCCACGCGTGGTCAGTTTCGACCAGCACATAGAAAGCAGa                                                                                                                                                      >  2:11066/1‑249 (MQ=255)                                                                                         acagtacaaggtcacgcGCGCCTTCCGCTGAGTTATCGGTTGACGAGAAGAACCAGATGCCGAGGCATGCCACGGTGATAGAGAAGAGCGCCGGGTATGCATAAGGGAAGATGGCTTTTTCGTGACCAAGGATCAGTACC‑CAAATCGACGAGACGAGGATCATCAGTACTACAACGGTAATCAGACCAAGCCAGCCACCCATCATCACGCCACGCGAGGTCAGGTTAGACCCGTACATAGAAAgacaactg                                                                                                                                                   >  2:584872/12‑244 (MQ=255)                                                                                           aGTTCACGCTCACGCGCTCCTTCCGCTGAGTTATCGGTTGCCGAGAAGAACCAGATGCCGAGGAATGCCACGGTGATAGAGAACAGCGCCGGGTATTCATACGGGAAGATGGCTTTTTCGTGACCAAGGATCTGTACC‑CAAATCGTCGAGCCGAGGATCATCAGTACTACTGCGGTAATCAGCCCCAGCCAGCCACCCATCATCGCGCCACGCGTGGTCAGTTTCGACCAGTACATAGAAAGCAGAATGAt                                                                                                                                                 <  1:536680/251‑1 (MQ=255)                                                                                                       cgcgcgCCTTCCGCTGAGTTATCGGTTGCCGAGAAGAACCAGATGCCGAGGAATGCCACGGTGATAGAGAACAGCGCCGGGAATTCATACGGGAAGATGGCTTTTTCGTGACCAAGGATCTGTACC‑CAAATCGTCGAGCCGAGGATCATCAGTACTACTGCGGTAATCAGCCCCAGCCAGCCACCCATCATCGCGCCACGCGTGGTCAGTTTCGACCAGGACATAGAAAGCAGAATGATCGGGAAGTTAc                                                                                                                                      >  1:437671/1‑250 (MQ=255)                                                                                                                           aTCGGTTGCCGAGAAGAACCAGATGCCGAGGAATGCCACGGTGATAGAGAACAGCGCCGGGTATTCATACGGGAAGATGGCTTTTTCGTGACCAAGGATCTGTACC‑CAAATCGTCGAGCCGAGGATCATCAGTACTACTGCGGTAATCAGCCCCAGCCAGCCACCCATCATCGCGCCACGCGTGGTCAGTTTCGACCAGTACATAGAAAGCAGAATGATCGGGAAGTTACAGCTCGCCGCGATGGCAAACg                                                                                                                 <  1:334013/251‑1 (MQ=255)                                                                                                                             cGGTTGCCGAGCAGAACCAGATGCCGAGGAATGACACGGTGATAGAGAACAGCGCAGGGTATTCATACGGGAAGCTGGCTTTTTCGTGCCCAAGGATCTGTACC‑CAACTCGTCGAGCCGAGGATCATCAGAACTACTGCGGAAGTCAGCCCCAGCCAGCCACCCATCATCGCGACACGCGTGGTCAGTTTCGACCCGTACATAGAAAGCAGCATGAACGGGACGTTACCGCTCGCCGCAGTGGAAAACGcc                                                                                                               >  2:101560/1‑251 (MQ=255)                                                                                                                              ggTTGCCGAGAAGAACCAGATGCCGAGGAATGCCACGGTGATAGAGAACAGCGCCGGGTATTCATACGGGAAGATGGCTTTTTCGTGACCAAGGATCTGTACC‑CAAATCGTCGAGCCGAGGATCATCAGTACTACTGCGGTAATCAGCCCCAGCCAGCCACCCATCATCGCGCCACGCGTGGTCAGTTTCGACCAGTACATAGAAAGCAGAATGATCGGGAAGTTACAGCTCGCCGCGATGGCAAACGcc                                                                                                               >  1:393589/1‑250 (MQ=255)                                                                                                                              ggTTGCCGAGAAGAACCAGATGCCGAGGAATGCCACGGTGATAGAGAACAGCGCCGGGTATTCATACGGGAAGAAGGCTATTTCGTGACCAAGGATCTGAACC‑CAAATCGTCGAGCCGAGGATCATCAGTACTACTGCGGTAATCAGACCCCGCCAGACACCCATACTCGAGCCAAGCGTAGTCAGGTTCGACCAGCACATAGACAGCAGAATGATCAGGATGTTACAGCTCGCCGCGATGGCAACCGcc                                                                                                               >  2:678977/1‑250 (MQ=255)                                                                                                                                     gagaAGAACCAGATGCCGAGGAATGCCACGGTGATAGAGAACAGCGCCGGGTATTCATACGGGAAGATGGCTTTTTCGTGACCAAGGATCTGTACC‑CAAATCGTCGAGCCGAGGATCATCAGTCCTACTGCGGTAATCAGCCCCAGCCAGCCACCCATCATCGCGCCCAGCGTGGTCAGTTTCGTCCAGAACATAGAAAGCAGAATGATCGGGAAGTTACCGCTCGCCGCGATGGCAAACGCACGCCccac                                                                                                       >  2:205459/1‑251 (MQ=255)                                                                                                                                      agaagaACCAGATGACGAGGCATGCCACGGTGATAGAGAACAGCGCCGGGTATTCATACGGGAAGATGGCTTTTTCGTGACCAAGGATCTGTACC‑CAAATCGTCAAGCCGAGGATCATCAGTACTACTGCGGTAATCAGCCCCAGCCAGCAACCCATCATCGCGCCTCGCGTGGTCAGTTTCGACCAGTACATAGAAAACAGAATGATCGGGAAGTTACAGATCGCCGCGATGGCAAAAGCCAGCCccaac                                                                                                      >  2:656594/1‑249 (MQ=255)                                                                                                                                       gaagaaCCAGATGCCGAGGAATGCAAAGGTGATAGAGAACAGCGCCGGGTATTCATACGGGAAGATGGCTTTTTCGTGACCAAGGATCTGTACC‑CAAATCGTCGAGCCGAGGATCATCAGTACTACTGCGGTAATCAGCCCCAGCCAGCCACCCATCATCGCGCCACGCGTGGTCAGTTTCGACCAGTACATAGAAAGCAGAATGATCGGGAAGTTACAGCTCGCCGCGATGGCAAACGCCAGCCccacca                                                                                                     <  1:89045/251‑1 (MQ=255)                                                                                                                                          ccaCCAGAGGCCGAGGAATGCCACGGTGAGAGAGAACAGCGCCGGGTATTCATACTGGAAGATGGCTTTTTCGTGACCAAGGACCTGTACC‑CAAAGCGGCGAGCCGAGGATCATCAGTACTACTGCGGTAATCAGCCCCAGCCAGCCACCCATCATCGCGCCACGCGTGGTCAGTTTCGACCAGTACATAGAAAGCAGAATGATCGGGAAGTTACAGCTCGCCGCGATGGCAAACGCCAGCCCCTCCATaa                                                                                                  <  2:501690/249‑1 (MQ=255)                                                                                                                                            aCCAGGTGCCGAGGCATGCCCCGGGGAGAGAGAACAGCGGCGGGTATTCATTCGGGAAGATGGCTTGTGCGTGACCAATGATCGGTTCC‑CAAATCGTCGAGCCGAGGATCAACAGTACTACTGCGGTAATCAGCCCCAGCCAGCCACCCATCATCGCGCCACGCGTGGTCAGTTTCGACCAGTACATAGAAAGCAGAATGATCTGGAAGTTACAGCTCGCCGCGATGGCAAACGCCAGCCCCACCATAAAg                                                                                                <  2:162271/251‑1 (MQ=255)                                                                                                                                                   gCCGAGGAATGCCACGGTGATAGAGAACAGCGCCGGGTATTCATACGGGAAGATGGCTTTTTCGTGACCAAGGATCTGTACC‑CAAATCGTCGAGCCGAGGATCATCAGTACTACTGCGGAAATCAGCCCCAGCCAGCCACCCATCATCGCGCCACGCGTGGTCAGTTTCGACCAGTACATAGAAAGCAGAATGATCGGGAAGTTACAGCTCGCCGCGATGGCAAACGCCAGCCCCCCCAAAAAGGCGATGt                                                                                         >  2:169874/1‑251 (MQ=255)                                                                                                                                                             gCCACGGTGATAGAGAACAGCGCCGGGTATTCATACGGGAAGATGGCTTTTTCGTGACCAAGGATCTGTACC‑CAAATCGTCGAGCCGAGGATCATCAGTACTACTGCGGTAATCAGCCCCAGCCAGCCACCCATCATCGCGCCACGCGTGGTCAGTTTCGACCAGTACATAGAAAGCAGAATGATCGGGAAGTTACAGCTCGCCGCGATGGCAAACGCCAGCCCCACCATAAAGGCGATGTTCTGATtctc                                                                               >  1:99781/1‑251 (MQ=255)                                                                                                                                                               taCGGTGATAGAGAACAGCGCCGGGTATTCATACGGGAAGATGGCTTTTTCGTGACCAAGGATCTGTACC‑CAAATCGTCGGGCCGAGGATCATCAGTACTACTGCGGTAATCAGCCCCAGCCAGCCACCCATCATCGCGCCACGCGTGGTCAGTTTc                                                                                                                                                                           <  2:193153/156‑1 (MQ=255)                                                                                                                                                               taCGGTGATAGAGAACAGCGCCGGGTATTCATACGGGAAGATGGCTTTTTCGTGACCAAGGATCTGTACC‑CAAATCGTCGGGCCGAGGATCATCAGTACTACTGCGGTAATCAGCCCCAGCCAGCCACCCATCATCGCGCCACGCGTGGTCAGTTTc                                                                                                                                                                           >  1:193153/2‑157 (MQ=255)                                                                                                                                                                      aaaggaaaCAGCGCATGGTTTTAATACGGAAAGATGGCTTTTTCTGGACCAAGGATCTGTACC‑AAAATTGTCGCTCCGAGGCTCATCAGTACTACTGCGGTAATCCGCCCCAGCCAGCCACCCATCATCGCGCCACGCGTGTTCAGTTTCGACCAGTACATAGAAAGCAGAATGATCGGGAAGTTACAGCTCGCCGCGATGGCAACCGCCAGCCCCACCATAAAGGCTATGTTCTGATTCTCA‑AACAGCAc                                                                     <  2:604843/245‑1 (MQ=255)                                                                                                                                                                       tAGAGAACAGCGCCGGGTATTCATACGGGAAGATGGCTTTTTCGTGACCAAGGATCTGTACA‑CAAATCGTCGAGCCGAGGATCATCAGTACTACTGCGGTAATCAGCCCCAGCCAGCCACCCATCATCGCGCCACGCGTGGTCAGTTTCGACCAGTACATAGAAAGCAGAATGATCGGGAAGTTACAGCTCGCCGCGATGGCAAACGCCAGCCCCACCATAAAGGCGATGTTCTGATTCTCA‑AACAGCACg                                                                    >  1:307297/1‑251 (MQ=255)                                                                                                                                                                           gaACAGCGCCGGGTATTCATACGGGAAGATGGCTTTTTCGTGACCAAGGATCTGTACC‑CAAATCGTCGCGCCGAGGATCATCAGTACTACTGCGGTAATCAGCCCCAGCCAGCCACCCATCATCGCGCCACGCGTGGTCAGTTTCGACCAGTACATAGAAAGCAGAATGATCGGGAAGTTACAGCTCGCCGCGATGGCAAACGCCAGCCCCACCATAAAGGCGATGTTCTGATTCTCA‑AACAGCACGCCg                                                                 >  1:312547/1‑250 (MQ=255)                                                                                                                                                                            aaCAGCGCCGGGTATTCATACGGGAAGATGGCTTTTTCGTGACCAAGGATCTGTACC‑CAAATCGTCGAGCCGAGGATCATCAGTACTACTGCGGTAATCAGCCCCAGCCAGCCACCCATCATCGCGCCACGCGTGGTCAGTTTCGACCAGTACATAGAAAGCAGAATGATCGGGAAGTTACAGCTCGCCGCGATGGCAAACGCCAGCCCCACCATAAAGGCGATGTTCTGATTCTCA‑AACAGCACGCCg                                                                 >  1:435148/1‑249 (MQ=255)                                                                                                                                                                               aGCGCCGGGTATTCATACGGGAAGATGGCTTTTTCGTGACCAAGGATCTGTACC‑CAAATCGTCGAGCCGAGGATCATCAGAACTACTGCGGTAATCAGCCCCAGCCAGCCCCCCATCATCGCGCCACGCGTGGTCAGTTTCGACCAGTACATAGAAAGCAGAATGATCGGGAAGTTACAGCTCGCCGCGATGGCAAACGCCAGCCCCACCATAAAGGCGATGATCTGATTCTCA‑AACAGCACGCCGAGGAt                                                            >  1:475062/1‑251 (MQ=255)                                                                                                                                                                                     gggTATTCATACGGGAAGATGGCTTTTTCGTGACCAAGGATCTGTACC‑CAAATCGTCGAGCCGAGGATCATCAGTACTACTGCGGTAATCAGCCCCAGCCAGCCACCCATCATCGCGCCACGCGTGGTCAGTTTCGACCAGTACATAGAAAGCAGAATGATCGGGAAGTTACAGCTCGCCGCGATGGCAAACGCCAGCCCCACCATAAAGGCGATGTTCTGATTCTCA‑AAAAGCACGCCGAGGATAATCGc                                                      >  1:21623/1‑251 (MQ=255)                                                                                                                                                                                             aTACGGGAAGATGGCTTTTTCGTGACCAAGGATCTGTACC‑CAAATCGTCGAGCCGAGGATCATCAGTACTACTGCGGTAATCAGCCCCAGCCAGCCACCCATCATCGCGCCACGCGTGGTCAGTTTCGACCAGTACATAGAAAGCAGAATGATCGGGAAGTTACAGCTCGCCGCGATGGCAAACGCCAGCCCCACCATAAAGGCGATGTTCTGATTCTCA‑AACAGCACGCCGAGGATAATCGCAATCACGc                                              >  1:128745/1‑251 (MQ=255)                                                                                                                                                                                             aTACGGGAAGATGGCTTTTTCGTGACCAAGGATCTGTACC‑CAAATCGTCGAGCCGAGGATCATCAGTACTACTGCGGTAATCAGCCCCAGCCAGCCACCCATCATCGCGCCACGCGTGGTCAGTTTCGACCAGTACATAGAAAGCAGAATGATCGGGAAGTTACAGCTCGCCGCGATGGCAAACGCCAGCCCCACCATAAAGGCGATGTTCTGATTCTCA‑AACAGCACGCCGAGGATAATCGCAATCACGc                                              >  1:181781/1‑251 (MQ=255)                                                                                                                                                                                                        tGGCTTTTTCGTGACCAAGGATCTGTACC‑CAAATCGTCGAGCCGAGGATCATCAGTACTACTGCGGTAATCAGCCCCAGCCAGCCACCCATCATCGCGCCACGCGTGGTCAGTTTCGACCAGTACATAGAAAGCAGAATGATCGGGAAGTTACAGCTCGCCGCGATGGCAAACGCCAGCCCCACCATAAAGGCGATGTTCTGACTCTCA‑AACAGCACGCCGAGGATAATCGCAATCACGCCGAGGATCAGt                                   >  1:24725/1‑251 (MQ=255)                                                                                                                                                                                                        tGGCTTTTTCGTGACCAAGGATCTGTACC‑CAAATCGTCGAGCCGAGGATCATCAGAACTACTGCGGTAATCAGCCCCAGCCAGCCACCCATCATCGCGCCACGCGTGGACAGTTTCGACCAGTACATAGAAAGCAGAATGATCGGGAAGTTACAGCTCGCCGCGATGGCAAGCGCCAGCCACACCCTAAAGGCGATGTTCTGATTCTCA‑AACAGAACGCCGAGGATCACAGCAATCCCGCCGGGGATCAGt                                   >  2:681264/1‑251 (MQ=255)                                                                                                                                                                                                           cTTTTTCGTGACCAAGGATCTGTACC‑CAAATCGTCGAGCCGAGGATCATCAGTACTACTGCGGTAATCAGCCCCAGCCAGCCACCCATCATCGCGCCACGCGTGGTCAGTTTCGACCAGTACATAGAAAGCAGAATGCTCGGGAAGTTACAGCTCGCCGCGATGGCAAACGCCAGCCCCACCATAAAGGCGATGTTCTGATTCTCA‑AACATCACGCCGAGGATAATCGCAATCACGCCGAGGATCAGTACg                                >  1:448821/1‑251 (MQ=255)                                                                                                                                                                                                                gcgggcactaggttcttaccc‑cAAAGCGCCGAGCCGAGGATCATCAGTACTACTGCGTTAATCAGCCCCAGCCAGCCACCCATCATCGCGCCACGCGTGGTCAGTTTCGACCAGTACATAGAAAGCAGAATGATCGGGAAGTTACAGCTCGCCGCGATGGCAAACGCCAGCACCACCATAAAGGCGATGTTCTGATTCTCA‑AACAGCACGCCGAGGATAATCGCAATCACGCCGAGGATCAGTACGGTGAt                           <  2:605251/232‑1 (MQ=255)                                                                                                                                                                                                                    gcccAAGGATCTGTACA‑CAAATCGTCGAGCCGAGGATCAGCAGAACTACTGCGGGAATCCGCACCAGCCAGCCACCCACCATCGCGCCACGCGTGGTCAGTTTCGACCAGTACATAGAAAGCAGAATGATCGGGAAGTTACAGCTCGCCGCGATGGCAAAAGCCAGCCCCACCATAAAGGCGATGTTATGATTC‑CACAACAGCACGCCGAGGATAATAGCAATCACGCCGGGGATCCGTCAGGTGAttttt                       >  2:633020/3‑251 (MQ=255)                                                                                                                                                                                                                                  aCC‑CAAATCGTCGAGCCGAGGATCATCAGTACTACTGCGGTAATCAGCCCCAGCCAGCCACCCATCATCGCGCCACGCGTGGTCAGTTTCGACCAGTACATAGAAAGCAGAATGATCGGGAAGTTACAGCTCGCCGCGATGGCAAACGCCAGCCCCACCATAAAGGCGATGTTCTGATTCCTA‑AACAGCACGCCGAGGATAATCGCATTCACGCCGAGGATCAATAAGGTGATTTTTGTTCCACGCAGctc         >  2:21655/1‑251 (MQ=255)                                                                                                                                                                                                                                         aTCGTCGAGCCGAGGATCATCAGTACTACTGCGGTAATCAGCCCCAGCCAGCCACCCATCATCGCGCCACGCGTGGTCAGTTTCGACCAGTACATAGAAAGCAGAATGATCGGGAAGTTACAGCTCGCCGCGATGGCAAACACCAGCCCCACCATAAAGGCGACGTTCTGATTCTCA‑ACCAGCACGCCGAGGATAATCGCAATCACGCCGAGGCTCAGTACGGAGCTTTTTGATACCCGCCGCTCttaccg   >  2:345276/1‑247 (MQ=255)                                                                                                                                                                                                                                           cgtcgAGCCGAGGATCATCAGTCCTACTGCGGTAATCAGCCCCAGCCAGCCACCCATCATCGCGCCACGCGTGGTCAGTTTCGACCAGTACATAGAAAGCAGAATGATCGGGAAGTTACAGCTCGCCGCGATGGCAAACGCCAGCCCCACCATAAAGGCGATGTTCTGATTCTCA‑AACAGCACGCCGAGGATAATCGCAATCACGCCGAGGATCAGTACGGTGATTTTTGATACCCGCAGCTCTTCACGt  <  1:155144/250‑1 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                      GTCAGACGACCGGGGAAGGGGAGAGATTAATGCGCGCGGCCTTGCTCAACGCCAAAGCCGGTCTGGGAGCGGATAAACTGCGCGCGGAACAGTTCACGCTCACGCGCGCCTTCCGCTGAGTTATCGGTTGCCGAGAAGAACCAGATGCCGAGGAATGCCACGGTGATAGAGAACAGCGCCGGGTATTCATACGGGAAGATGGCTTTTTCGTGACCAAGGATCTGTACC‑CAAATCGTCGGGCCGAGGATCATCAGTACTACTGCGGTAATCAGCCCCAGCCAGCCACCCATCATCGCGCCACGCGTGGTCAGTTTCGACCAGTACATAGAAAGCAGAATGATCGGGAAGTTACAGCTCGCCGCGATGGCAAACGCCAGCCCCACCATAAAGGCGATGTTCTGATTCTCA‑AACAGCACGCCGAGGATAATCGCAATCACGCCGAGGATCAGTACGGTGATTTTTGATACCCGCAGCTCTTCACGT  >  NC_000913/4283227‑4283709

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