Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,955,394	C→T	10.2%	A538A (GCG→GCA)	torZ ←	trimethylamine N‑oxide reductase 2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,955,394	0	C	T	10.2%	113.2 / 3.2	39	A538A (GCG→GCA)	torZ	trimethylamine N‑oxide reductase 2
Reads supporting (aligned to +/- strand):  ref base C (17/18);  new base T (2/2);  total (19/20)
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.77e-01

ATTTTTCACCGGATCGGCGCGGAAATCACCATAACGAAC‑GTACTGTTCGTTCTTCTCGCTGTGGCGCATTTCGATCAGTTTATTTTGCTGCCAGAAGGCATTAAACATTGGCATAGTGACGCGTTGCGCACGGGCACCTTTCTGAGCGGCATCATAGAAAAATTTCAGCCACGCCATTTCATCTTTACCTTCGGTATAGATCTCTTTTCCGCCAGGTTTGAGTAATTCCGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTGCTGTAATCACCGGTCATCGTCAAGTCATTGCGCTCAAACGATGTGGTGATCGGTAATACGATATCTGCGTGTTTGGCTGCTGCGGTCCAGTAGCATTCAGAAACGACGATCATCTCCGGTTTCTGCCATGCTTTAATCAGACGATTAGTATCCTGGTGATGGGTAAAGTTACCGCCACCCGCCCACCAGATCATTTTGATATTCGGATACGT  >  NC_000913/1955133‑1955657                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                        aTTCTTCACCGGATCGGCGCGGAAATCACCATAACGAAC‑GTACTGTTCGTTCTTCTCGCTGTGGCGCATTTCGATCAGTTTATTTTGCTGCCAGAAGGCATTAAACATTGGCATAGTGACGCGTTGCGCACGGGCACCTTTCTGAGCGGCATCATAGAAAAATTTCAGCCACGCCATTTCATCTTTACCTTCGGTATAGATCTCTTTTCCGCCAGGTTTGAGTAATTCCGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAAt                                                                                                                                                                                                                                                                >  1:28598/1‑271 (MQ=255)       ctccGGATCGGCGCGGAAATCACCATAACGATC‑GTACGTTTCGTTCTTCTCGCTGTGGCGCATTTCGATCAGTTTATTTTGCTGCCAGAAGGCATTAAACATTGGCATAGTGACGCGTTGCGAACGGGCACCTTTCTGAGCGGCATCATAGAAAAATTTCAGCCACGCCATTTCATCTTTACCTTCGGTATAGATCTCTTTTCCGCCAGGTTTGAGTAATTCCGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAg                                                                                                                                                                                                                                                         <  1:178524/270‑1 (MQ=255)          cGGATCGGCGCGGAAATCACCATACCGATC‑TTACTTTTCGTTCTTCTCGCTGTGGCGCATTTCGATCAGTTTATTTTGCTGCCAGAAGGCATTAAACATTGGCATAGTGACGCGTTGCGCACGGGCACCTTTCTGAGCGGCATCATAGAAAAATTTCAGCCACGCCATTTCATCTTTACCTTCGGTATAGATCTCTTTTCCGCCAGGTTTGAGTAATTCCGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCaa                                                                                                                                                                                                                                                      <  2:33785/272‑1 (MQ=255)                            aCCCTAACGCACGGTTCTTTTCTTTTTTTTCGCTGTGGCTCATTTCGATCAGTTTATTTTGCTGCCAGAAGGCATTAAACATTGGCATAGTGACGCGTTGCGCACGGGCACCTTTCTGAGCGGCATCATAGAAAAATTTCAGCCACGCCATTTCATCTTTACCTTCGGTATAGATCTCTTTTCCGCCAGGTTTGAGTAATTCCGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGc                                                                                                                                                                                                                                     <  2:28598/272‑1 (MQ=255)                                            cTGTTCGTTCTTCTCGCTGTGGCGCATTTCGATCAGTTTATTTTGCTGCCAGAAGGCATTAAACATTGGCATAGTGACGCGTTGCGCACGGGCACCTTTCTGAGCGGCATCATAGAAAAATTTCAGCCACGCCATTTCATCTTTACCTTCGGTATAGATCTCTTTTCCGCCAGGTTTGAGTAATTCCGCAAGATCGTCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATgaatgttgc                                                                                                                                                                                                                     >  1:218516/1‑264 (MQ=255)                                                  gtttTTCTCGCTTGGGCGCATTTCGATCATTTTATTTTGCTGCCAGAAGGCATTAAACATTGGCATAGTGACGCGTTGCGCACGGGCACCTTTCTGAGCGGCATCATAGAAAAATTTCAGCCACGCCATTTCATCTTTACCTTCGGTATAGATCTCTTTTCCGCCAGGTTTGAGTAATTCCGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTGCTGTAAt                                                                                                                                                                                                               <  2:286978/272‑1 (MQ=255)                                                    tcttctCGCTGTGGCGCATTTCGATCAGTTTATTTTGCTGCCAGAATGCATTAAACATTGGCATAGTGACGCGTTGCGCACGGGCACCTTTCTGAGCGGCATCATAGAAAAATTTCAGCCACGCCATTTCATCTTTACCTTCGGTATAGATCTCTTTTCCGCCAGGTTTGAGTAATTTCGCAAGATCGGCAAACACGTCACAATCGTTACGCGCTTCAAATTACGGAGCAACAGCCTGCTTCATCGGCACAATATGCTAGTTGCTGTaa                                                                                                                                                                                                                >  1:286978/1‑269 (MQ=255)                                                       tctcGCTGTGGCGCATTTCGATCCGTTTATTTTGCTTCCAGAAGGCATTAAACATTGGCATAGTGACGCGTTGCGCACGGGCACCTTTCTGAGCGGCATCATAGAAAAATTTCAGCCACGCCATTTCATCTTTACCTTCTGTATAGATCTCTTTTCCGCCAGGTTTGAGTAATTCCGCAAGATCGGCAAACACGTCAAAATCGTTACgcgc                                                                                                                                                                                                                                                                       >  1:191050/1‑211 (MQ=255)                                                                   gcATTTCGATCAGTTTATTTTGCTGCCAGAAGGCATTAAACATTGGCATCGTGACGCGTTGCGCACGGGCACCTTTCTGAGCGGCATCATAGAAAAATTTCAGCCACGCCATTTCATATTTACCTTCGGTATAGATCTCTTTTCCGCCAGGTTTGAGTAATTCCGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTCCGGAGCAACCGCCTGCTTCATCGGCACAATATGCTGGTTTCTGTAATCACCGGTCATCGTCaaa                                                                                                                                                                                              >  1:376955/1‑271 (MQ=255)                                                                     aTTTCGATCAGTTTATTTTGCTGCCAGAAGGCATTAAACATTGGCATAGTGACGCGTTGCGCACGGGCACCTTTCTGAGCGGCATCATAGAAAAATTTCAGCCACGCCATTTCATCTTTACCTTCGGTATAGATCTCTTTTCCGCCAGGTTTGAGTAATTCCGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTGCTGGAACTAACGGTCAACGTCAAGTc                                                                                                                                                                                            >  2:97571/1‑272 (MQ=255)                                                                          caTCAGTTTTTTTTGCTCCCAGAAGGCATGAAACATTGGCAGAGTGACGCGCTGCGCACGGGAACCTTTCTGAGCGGCATCATAGCAAAATTTCAGCCACGCCATTTCATCTTTACCTTCGGTATAGATCTCTTTTCCGCCAGGTTTGAGTAATTCCGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTGCTGTAATCACCGGTCATCGTCAAGTCATTgc                                                                                                                                                                                       <  2:132259/271‑1 (MQ=255)                                                                                   ttttgGCTGCCGGAAGGCGTTAAACATTGGCATGGTGCAGCGTTGCGAACGGGCACCTTTCTGAGCGGCATCATAGAAAAATGTCAGCCACGCCATTTCATCTTTACCTTCGGTATAGATCTCTTTTCCGCCAGGTTTGAGTAATTCCGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTGCTGTAATCACCGGTCATCGTCAAGTCATTGCGCTCAAAc                                                                                                                                                                               <  2:376955/270‑1 (MQ=255)                                                                                        gctgcCAGAAGGCATTAAACATTGGCATAGTGACGCGTTGCGCACGGGCACCTTTCTGAGCGGCATCATAGAAAAATTTCAGCCACGCCATTTCATCTTTACCTTCGGTATAGATCTCTTTTCCGCCAGGTTTGAGTAATTCCGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTGCTGTAATCACCGGTCATCATCAAGTCATTGCGCTCAAACGAtgtg                                                                                                                                                                         <  1:97571/272‑1 (MQ=255)                                                                                                         aaCATTGGCATAGTGACGCGTTGCGCACGGGCACCTTTCTGAGCGGCATCATAGAAAAATTTCAGCCACGCCATTTCATCTTTACCTTCGGTATAGATCTCTTTTCCGCCAGGTTTGAGTAATTCCGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGa                                                                                                                                                                                                                                                          >  1:119434/1‑174 (MQ=255)                                                                                                          aCATTGGCATAGTGACGCGTTGCGCACTGGCACCTTTCTGAGCGGCATCATAGAAAAATTTCAGCCACGCCATTTCATCTTTACCTTCGGTATAGATCTCTTTTCCGCCAGGTTTGAGTAATTCCGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAg                                                                                                                                                                                                                                                         <  2:119434/174‑1 (MQ=255)                                                                                                                 cATAGTGACGCGTTGCGCACGGGCACCTTTCTGAGCGGCATCATAGAAAAATTTCAGCCACGCCATTTCATCTTTACCTTCGGTATAGATCTCTTTTCCGCCAGGTTTGAGTAATTCCGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTGCTGTAATCACCGGTCATCGTCAAGTCATTGCGCTCAAACGCTGTGGTTATCGGTAACACGAAATcgtcgt                                                                                                                                                >  2:130864/1‑267 (MQ=255)                                                                                                                 cATAGTGACGCGTTGCGCACGGGCACCTTTCTGAGCGGCATCATAGAAAAATTTCAGCCACGCCATTTCATCTTTACCTTCGGTATAGATCTCTTTTCCGCCAGGTTTGAGTAATTCCGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATACGCTGGTTGCTGTAATCACCGGTCCTCGTCAAGTCATTGCTACCAAACGCAGGGGTGACCGGTAAAACGataatcccgc                                                                                                                                                >  2:322769/1‑265 (MQ=255)                                                                                                                  aTAGTGACGCGTTGCGCACGGGCACCTTTCTGAGCGGCATCATAGAAAAATTTCAGCCACGCCATTTCATCTTTACCTTCGGTATAGATCTCTTTTCCGCCAGGTTTGAGTAATTCCGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTGCTGTAATCACCGGTCATCGTCAAGTCATTGCGCTCAAACGATGTGGTGATCGGTAATACGATATCt                                                                                                                                                    <  1:130864/267‑1 (MQ=255)                                                                                                                                   aCGGGCACCTTTCTGAGCGGCATCATAGAAAAATTTCAGCCACGCCATTTCATCTTTACCTTCGGTATAGATCTCTTTTCCGCCAGGTTTGAGTAATTCCGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTGCTGTAATCACCGGTCATCGTCAAGTCATTGCGCTCAAACGATgtgg                                                                                                                                                                        >  1:55630/1‑230 (MQ=255)                                                                                                                                    cGGGCACCTTTCTGAGCGGCATCATAGAAAAATTTCAGCCACGCCATTTCATCTTTACCTTCGGTATAGATCTCTTTTCCGCCAGGTTTGAGTAATTCCGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTGCTGTAATCACCGGTCATCGTCAAGTCATTGCGCTCAAACGATgtggt                                                                                                                                                                       <  2:55630/230‑1 (MQ=255)                                                                                                                                                    cGGCATCATAGAAAAATTTCAGCCACGCCATTTCATCTTTACCTTCGGTATAGATCTCTTTTCCGCCAGGTTTGAGTAATTCCGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTGCTGTAATCACCGGTCATCGTCAAGTCATTGCGCTCAAACGAtgt                                                                                                                                                                          >  1:327516/1‑211 (MQ=255)                                                                                                                                                    cGGCATCATAGAAAAATTTCAGCCACGCCATTTCATCTTTACCTTCGGTATAGATCTCTTTTCCGCCAGGTTTGAGTAATTCCGCAAGATCGGCAAACACGACAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTGCTGTAATCACCGGTCATCGTCACGTCATTGCCCTAAAACGATGTGGTGATCGGTAATACGATATCTTCGTGTTTGgcc                                                                                                                                        >  1:228578/1‑244 (MQ=255)                                                                                                                                                     ggCATCATAGAAAAATTTCAGCCACGCCATTTCATCTTTACCTTCGGTATAGATCTCTTTTCCGCCAGGTTTGAGTAATTCCGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTGCTGTAATCACCGGTCATCGTCAAGTCATTGCGCTCAAACGAtgtg                                                                                                                                                                         <  2:327516/211‑1 (MQ=255)                                                                                                                                                        atcatAGAAAAATTTCAGCCACGCCATTTCATCTTTACCTTCGGTATAGATCTCTTTTCCGCCAGGTTTGAGTAATTCCGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATg                                                                                                                                                                                                                             >  2:270197/1‑156 (MQ=255)                                                                                                                                                         tcatAGAAAAATTTCAGCCACGCCATTTCATCTTTACCTTCGGTATAGATCTCTTTTCCGCCAGGTTTGAGTAATTCCGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGc                                                                                                                                                                                                                            <  1:270197/156‑1 (MQ=255)                                                                                                                                                          catAGAAAAATTTCAGCCACGCCCTTTCATCTTTACCTTCGATATAGATCCCTTTTCTGCCAGGACTTAGTAATTCCGCAAGAACGGCAAACACGTCAAAATCGTTACGTGCATCATATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTGTATGCTGTACTCACCGGTCATCGTCAAGTCATTGCGCTCAAAAGAGGTGGTTATCGGTAATACGATAACAGCGTGGTTGGCTACTTGTGTCCCGTAGAAGTCAGaaaacacaaca                                                                                                       >  2:179762/1‑264 (MQ=255)                                                                                                                                                          catAGAAAAATTTCAGCCACGCCATTTCATCTTTACCTTCGGTATAGATCTCTTTTCCGCCAGGTTTGAGTAATTCCGCAAGATCGGCAAACACGTCAAAATCGTTACGTGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTGCTGTAATCACCGGTCATCGTCAAGTCATTGCGCTCAAACGAGTTGGTTAACGGTAAAACGATATTGGCGTTTTTGGCTGCTGGGGTCCAGTAGACTTCAaaaaacaaaaac                                                                                                       >  2:179755/1‑264 (MQ=255)                                                                                                                                                          catAGAAAAATTTCAGCCACGCCATTTCATCTTTACCTTCGGTATAGATCTCTTTTCCGCCAGGTTTGAGTAATTCCGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTGCTGTAATCACCGGTCATCGTCAAGTCATTGCGCTCAAACGATGTGGTGATCGGTAATACGATATCTGCGTGTTTGGCTGCTGCGGTCCAGTTGCACTCAGAAacgacggct                                                                                                       >  2:28303/1‑269 (MQ=255)                                                                                                                                                                 aaaTTTCAGCCCCGCCATTTCATCTTTACCTTCGGTATAGATCTCTTTTCCGCCAGGTTTGAGTAATTCCGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTGCTGTAATCACCGGTCATCGTCAAGTCATTGCGCTCAAACGATGTGGTGATCGGTAATACGATATCTGCGTGTTTGGCTGCTGCGGTCCAGTAGCATTCAGAAACGACGATCATCTcc                                                                                                 <  1:28303/271‑1 (MQ=255)                                                                                                                                                                                                            tctTTTCCGCCAGGTTTGAGTAATTCCGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTGCTGTAATCACCGGTCATCGTCAAGTCATTGCGCTCAAACGATGTGGTGATCGGTAATACGATATCTGCGTGTTTGGCTGCTGCGGTCCAGTAGCATTCAGAAACGACGATCATCCCCGGTTTCTGCCATGCTTTAATCAGACGATTCGTATCCAGGTTAt                                                      >  2:277230/1‑271 (MQ=255)                                                                                                                                                                                                                      cAGGTTTGAGTAATTCCGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTGCTGTAATCACCGGTCATCGTCAAGTCATTGCGCTCAAACGATGTGGTGATCGGt                                                                                                                                                                >  2:424837/1‑155 (MQ=255)                                                                                                                                                                                                                       aGGTTTGAGTAATTCCGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTGCTGTAATCACCGGTCATCGTCAAGTCATTGCGCTCAAACGATGTGGTGATCGGTa                                                                                                                                                               <  1:424837/155‑1 (MQ=255)                                                                                                                                                                                                                             gAGTAATTCCGCATGTGCTGCAAACACGTCAAATGTGTTACGCTCTTCAAATTGCGGATCAACAGCCTGCTTCATCGGCACAATGTGCTTGTTGCTGTAATCACCTGTCATCGTCAAGTCATTGCGCTCAAACGATGTGGTGATCGTTAATACGATATCTGCGTGTTTGGCTGCTTCGGTCCAGTAGCATTCAGAAACGACGATCATCTCCGGTTTCTGCCATGCTTTAATCAGACGATTAGTATCCTGGTGATGGGTAAAGTTACCGCCAc                                    <  1:348825/272‑1 (MQ=255)                                                                                                                                                                                                                                                       cgtataaaTCGTTACGTGCTTCAAATTGCGGAGCAACTGCCTGCTTCATCGGCACAATATGCTGGTTGCTGTAATCACCGGTCATCGGCAAGTCATTGCGCTCAAACGATGTGGTGATCGGTAATACGATATCTGCGTGTTTGGCTGCTGCGGTCCAGTAGCATTCAGAAACGACGATCATCTCCGGTTTCTGCCATGCTTTAATCAGACGATTAGGATCCTGGTGATGGGTAAAGTTACCGCCACCCGCCCACCAGATCATTTTGATATTc          <  1:179762/267‑1 (MQ=255)                                                                                                                                                                                                                                                       cGTCAAAATCGTTACGTGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTGCTGTAATCACCGGTCATCGTCAAGTCATTGCGCTCAAACGATGTGGTGATCGGTAATACGATATCTGCGTGTTTGGCTGCTGCGGTCCAGTAGCATTCAGAAACGACGATCATCTCCGGTTTCTGCCATGCTTTAATCAGACGATTAGTATCCTGGTGATGGGTAAAGTTACCGCCACCCGCCCACCAGATCATTTTGATATTc          <  1:179755/272‑1 (MQ=255)                                                                                                                                                                                                                                                             aaTCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTGCTGTAATCACCGGTCATCGTCAAGTCATTGCGCTCAAACGATGTGGTGATCGGTAATACGATATCTGCGTGTTTGGCTGCTGCGGTCCAGTAGCATTCAGAAACGACGATCATCTCCGGTTTCTGCCATGCTTTAAt                                                                             >  2:84189/1‑199 (MQ=255)                                                                                                                                                                                                                                                              aTCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTGCTGTAATCACCGGTCATCGTCAAGTCATTGCGCTCAAACGATGTGGTGATCGGTAATACGATATCTGCGTGTTTGGCTGCTGCGGTCCAGTAGCATTCAGAAACGACGATCATCTCCGGTTTCTGCCATGCTTTAATc                                                                            <  1:84189/199‑1 (MQ=255)                                                                                                                                                                                                                                                                cGTTACGCGCTTCAAATTGCTGGGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTGCTGTAATCACCGGTCATCGTCAAGTCATTGCGCTCAAACGATGTGGTGATCGGTAATACGATATCTGCGTGTTTGGCTGCTGCGGTCCAGTAGCATTCAGAAACGACGATCATCTCCGGTTTCTGCCATGCTTTAATCAGACGATTAGTATCCTGGTGATGGGTAAAGTTACCGCCACCCGCCCACCAGATCATTTTGATATTCGGATACGt  <  1:277230/271‑1 (MQ=255)                                                                                                                                                                                                                                                                   tACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTGCTGTAATCACCGGTCATCGTCAAGTCATTGCGCTCAAACGATGTGGTGATCGGTAATACGATATCTGCGTGTTTGGCTGCTGCGGTCCAGTAGCATTCAGAAACGACGATCATCTCCGGTTTCTGCCATGCTTTAATCAGa                                                                         >  2:343109/1‑197 (MQ=255)                                                                                                                                                                                                                                                                    aCGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTGCTGTAATCACCGGTCATCGTCAAGTCATTGCGCTCAAACGATGTGGTGATCGGTAATACGATATCTGCGTGTTTGGCTGCTGCGGTCCAGTAGCATTCAGAAACGACGATCATCTCCGGTTTCTGCCATGCTTTAATCAGAc                                                                        <  1:343109/197‑1 (MQ=255)                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                        ATTTTTCACCGGATCGGCGCGGAAATCACCATAACGAAC‑GTACTGTTCGTTCTTCTCGCTGTGGCGCATTTCGATCAGTTTATTTTGCTGCCAGAAGGCATTAAACATTGGCATAGTGACGCGTTGCGCACGGGCACCTTTCTGAGCGGCATCATAGAAAAATTTCAGCCACGCCATTTCATCTTTACCTTCGGTATAGATCTCTTTTCCGCCAGGTTTGAGTAATTCCGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTGCTGTAATCACCGGTCATCGTCAAGTCATTGCGCTCAAACGATGTGGTGATCGGTAATACGATATCTGCGTGTTTGGCTGCTGCGGTCCAGTAGCATTCAGAAACGACGATCATCTCCGGTTTCTGCCATGCTTTAATCAGACGATTAGTATCCTGGTGATGGGTAAAGTTACCGCCACCCGCCCACCAGATCATTTTGATATTCGGATACGT  >  NC_000913/1955133‑1955657

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑