Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,993,531	T→C	R793R (CGT→CGC)	cyaA →	adenylate cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,993,531	0	T	C	94.4%	110.2 / ‑3.5	38	R793R (CGT→CGC)	cyaA	adenylate cyclase
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base C (18/16);  total (19/19)
Fisher's exact test for biased strand distribution p-value = 4.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.27e-01

ATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCAT‑GACCGTTTTACCTACGGCTCAAGCTTCATCAACTTC‑AACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCT  >  NC_000913/3993287‑3993767                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                              aTAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGGGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCAT‑GACCGCTTTAcc                                                                                                                                                                                                                                         <  2:367616/251‑1 (MQ=255)              cGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGGCGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTTGAAGAAACGCAAGACGAGGAAGGCTTTTATATCTACATTCTCGACGAAAGCAACCCGGGTGAGGGATATTACCACTGCGAAGGCAGCAAAGAGGAGCCGGTAAGTGACGTCAGTCGCTTCTACTCGCCATCGCAT‑TACCGCTTTACCTACGGCTCAAGCt                                                                                                                                                                                                                            >  1:144885/1‑251 (MQ=255)                  cTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACAGCCTTGCCAACGAAGGGATAATCCCGTTGTTTTTCGTAAAAAAGCAAGAAAAGAAACGCCTTAATATTTACATTCTCGGCGAAAACAAACGGGTTGAGGTATATTACCCCCGCGAACGCAACACAACGGAACTGATCCGTGACGTCACTCGCTTTTACTCGTCATCGCCT‑CACCCGTTTCCCCACGGCTCAACCtttc                                                                                                                                                                                                                         >  1:252162/1‑248 (MQ=255)                    gTCAGGGCAGGGTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCAT‑GACCGCTTTACCTACGGCTCAAGCTtcatca                                                                                                                                                                                                                      <  2:251870/251‑1 (MQ=255)                            aGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCTCCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCAT‑GACCGCTTTACCTACGGCTCAAGCTTCATCAACTTC‑AAc                                                                                                                                                                                                             <  2:81835/251‑1 (MQ=255)                            aGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATAACTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCAT‑GACCGTTTTACCTACGGCTCAAGCTTCATCAACTTC‑AAc                                                                                                                                                                                                             <  1:60882/251‑1 (MQ=255)                            aGGTTGAAACCAATCACGTCAAATTACCGGCGGCGGTGGCAGGCTCTGGTCGAGAAGGGGGCATCCAGTTCTTTTTTGATGGAACGCAAAACGAGAGTGGCTTTAATCTATACATTTTTGCAGCAAGCAACCGGGTGGAGGGATATCACCCCTGCGAAGGCAGCAAAGGGGAGGGGGTACCGGACGTCCGCCGCTTCTACACGTCATCGAATAGAC‑GCTTTACCTACgtctacagattcttaagcttcacc                                                                                                                                                                                                              >  2:210093/1‑228 (MQ=255)                              gTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCAT‑GACCGCTTTACCTACGGCTCAAGCTTCATCAACTTC‑AACCt                                                                                                                                                                                                           <  1:183933/251‑1 (MQ=255)                                 gAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCCTCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATTTCTACATTCTCGACGAAAGCAACCGGGTTGTGGGATATCACCACTGCGAAGGCAGCAAAGAGGGGCTGGGACGTGACGTCAGTCGCTTCTACTCGGCATCGCAT‑GACCGCTTTACCTACCGCTCAAGCTTCATCAACTTC‑AACCTGcc                                                                                                                                                                                                        >  1:119633/1‑251 (MQ=255)                                          cACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATACCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCAT‑GACCGCTTTACCTACGGCTCAAGCTTCATCAACTTC‑AACCTGCCGCAGTTCt                                                                                                                                                                                                <  1:187093/250‑1 (MQ=255)                                                aaagtACCGGCGGTGTTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATAGCTGCATGCGCGACGAAAGAAACCGGGTTGAGGTAGATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGGCAGCGCAT‑GACCGCTTTACCTACGGCTCAAGCTTCATCAACCTC‑AACCTGCCGCAGTTCTATCAGAt                                                                                                                                                                                         <  2:48916/247‑1 (MQ=255)                                                      ccGGCGGTGGTGGTCGGCTTTGCCATTGAAGGGATCATCCAGGTCTGTGTCGAAGAAACGCAAGACGAGAATGGCTTTAATAGCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGTTACGTGACGTCAGTCGCTTCTACTCGTCATCGCAT‑GACCGCTTTACCTTCGGCTCAAGCTTAATCAACTTC‑AACCTGCCGCAGTTCTCTCAGATTGTGaa                                                                                                                                                                                   <  2:117769/251‑1 (MQ=255)                                                                     ggCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATAACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCAT‑GACCGCTTTACCTACGGCTCCAGCTTCATCACCATC‑AACCTGCCGCAGTTCTATCAGATTGTGCAGGTTGATGGTcctg                                                                                                                                                                     >  2:7969/1‑247 (MQ=255)                                                                                 gAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCAT‑GCCCGCTTTACCTACGGCTCAAGCTTCATCAACTTC‑AACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGAGAttcc                                                                                                                                                        >  1:335333/1‑251 (MQ=255)                                                                                              aGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCAT‑GACCGCTTTACCTACGGCTCAAGCTTCATCAACTTC‑AACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACaaaa                                                                                                                                           <  2:329283/251‑1 (MQ=255)                                                                                                                cGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCAT‑GACCGCTTTACCTACGGCTCAAGCTTCATCAACTTC‑AACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATgccg                                                                                                                         >  1:66313/1‑251 (MQ=255)                                                                                                                         agcaTGGCTTTCATATCTACATTCTCGACGAAAGCAACCGGGGTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGAAGCTGGTACGTGACGTAAGTCGCTTCTACTCGTCATCGCAT‑GACCGCTTTACCTACGGCTCAAGCTTCATCCACTTC‑AACCTGCCGCAGTTCTATCAGATTGTGAGGGTTGATGGTCGTGAACAGGTGCTTCCGTTCCGCACAACATCTATCGGTAACATGCCGCCTGCCaac                                                                                                                >  2:246780/4‑250 (MQ=255)                                                                                                                                 tttAATATATACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGACGGCAGCAAAGCGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCAA‑GACCGCTTTACCTACGGCTCACGCTTCATCAACTTC‑AACCTGCCGAAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGAGATTCCGTTCGGCACAAAATCTATCGGTAACATGCCGCCTGCCTATCAGGAta                                                                                                         >  2:213104/1‑249 (MQ=255)                                                                                                                                  ttAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCAT‑GACCGCTTTACCTACGGCTCAAGCTTCATCAACTTC‑AACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATc                                                                                                         <  1:312863/249‑1 (MQ=255)                                                                                                                                     atctCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCCCTACGCATGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCAT‑GACCGCTTTACCTACGGCTCAAGCTTCATCAACTTC‑AACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCATTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACgct                                                                                                    >  2:171997/4‑249 (MQ=255)                                                                                                                                           aCATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCAT‑GACCGCTTTACCTACGGCTCAAGCTTCATCAACTTC‑AACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATAc                                                                                                  <  1:38870/247‑1 (MQ=255)                                                                                                                                            cATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCAT‑GACCGCTTTACCTACGGCTCAAGCTTCATCAACTTC‑AACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGTTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATAcgccgc                                                                                             >  1:171946/1‑251 (MQ=255)                                                                                                                                                        aaggaacggggTTGAGGGAGAGCACCAGTGCGAAGGCAGCAAAGAGGAGCTTGTACGTTACGTCAGTCGCTGCTATTCGTCATCGCCT‑GTCCGACCTTCCTACGGCTCAAGCTTAATCAATTGT‑AACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACCCGTTATTCCGTTCCGCACAAAATCTATCGGTACCATGCCTCCTGCCAATCAGGATCACGAGACGCCGCTATTACAGCAat                                                                                 <  2:405031/243‑1 (MQ=255)                                                                                                                                                            aaCCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCAT‑GACCGCTTTACCTACGGCTCAAGCTTCCTCAACTTC‑AACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACCTGCCGCCTGCAAATCAGGCTCACGATAAGCCGCTATTACAGCAATAttt                                                                             >  2:28489/1‑251 (MQ=255)                                                                                                                                                             aCCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCAT‑GACCGCTTTACCTACGGCTCAAGCTTCATCAACTTC‑AACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCCCAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGGTACGCCGCTATTACAGCAATAtttt                                                                            >  1:328442/1‑251 (MQ=255)                                                                                                                                                              ccGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGGACGTGACGTCAGTCGCTTCTACTCGTCATCGCAT‑GACCGCTTTACCTACGGCTCAAGCTTCATCAACTTC‑AACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAACCAGGATCACGATACGCCGATATTACAGCAATAttttt                                                                           >  2:144018/1‑251 (MQ=255)                                                                                                                                                                            atCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCAT‑GACCGCTTTACCTACGGCTCAAGCTTCATCAACTTC‑AACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGc                                                             >  2:292110/1‑251 (MQ=255)                                                                                                                                                                             tCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCAT‑GACCGCTTTACCTACGGCTCAAGCTTCATCAACTTC‑AACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGcc                                                            >  1:84808/1‑251 (MQ=255)                                                                                                                                                                                                 aaGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCAT‑GACCGCTTTACCTACGGCTCAAGCTTCATCAACTTC‑AACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCgg                                        <  1:285884/251‑1 (MQ=255)                                                                                                                                                                                                   gaggagCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCAT‑GACCGCTTTACCTACGGCTCAAGCTTCATCAACTTC‑AACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGc                                       >  1:380120/1‑250 (MQ=255)                                                                                                                                                                                                         cTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCAT‑GACCGCTTTACCTACGGCTCAAGCTTCATCAACTTC‑AACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCACTCAGGATCACGATACGCCGCTATTACAGCACTATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTCTCCGGCATGCAAt                                >  2:180018/1‑251 (MQ=255)                                                                                                                                                                                                                    cGTCAGTCGTTTCTACTCGTACTCGCAT‑GCCCTCTTTCCCTACGGCTCAAGCTTCATCATCTTC‑AACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGaaa                     <  2:119633/251‑1 (MQ=255)                                                                                                                                                                                                                    cGTCAGTCGCTTCTACTCGTCATCGCAT‑GACCGCTTTACCTACGGCTCAAGCTTCATCAACTTC‑AACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGaaa                     <  1:292110/251‑1 (MQ=255)                                                                                                                                                                                                                     gTCAGTCGCTTCTACTCGTCATCGCAT‑GACCGCTTTACCTACGGCTCAAGCTTCATCAACTTC‑AACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAac                    >  1:2181/1‑251 (MQ=255)                                                                                                                                                                                                                                  aCTCGTCACCGCAT‑GACCGCTTTACCTACGGCTCAAGCTTCATAAACTTC‑AACCTGCCGCAGTTCTAGCAGATTGTGAAGGTTGATGGGCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCAc        <  2:15746/250‑1 (MQ=255)                                                                                                                                                                                                                                     cGTCATCGCAT‑GACCGCTTTACCTACGGCTCAAGCTTCATCAACTTC‑AACCTGCCGCAGTTCTATCAGATTGTTAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCg     <  1:235870/250‑1 (MQ=255)                                                                                                                                                                                                                                        cagcGCAT‑GCCCGTTTTACCCACGGCTCAAGCTTCA‑CAAATTCCACCCTGTCGCAGTTTTATTAGTGGGTGAAGGTTGATGGTTGTGAACAGGTGATTCCGTTCCGCAAAAAATCTTTCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACTCGCct  <  2:328442/247‑1 (MQ=255)                                                                                                                                                                                                                                        cATCGCAT‑GACCGCTTTACCTACGGCTCAAGCTTCATCAACTTC‑AACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACAGGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCct  <  1:225964/250‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                              ATAACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCAT‑GACCGTTTTACCTACGGCTCAAGCTTCATCAACTTC‑AACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCT  >  NC_000913/3993287‑3993767

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