Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	700,038	G→A	H97Y (CAT→TAT)	umpH ←	UMP phosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	700,038	0	G	A	92.3%	102.0 / 1.2	39	H97Y (CAT→TAT)	umpH	UMP phosphatase
Reads supporting (aligned to +/- strand):  ref base G (1/2);  new base A (22/14);  total (23/16)
Fisher's exact test for biased strand distribution p-value = 5.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.61e-01

ACGGGCTGGGCTTACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGC‑ACACAACGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATGAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGA‑AATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGATCTTGCCCAGTCTGCGAAGGATAGTTGGTCAGCAACACCAGCGGCAGGCCTTTATCCATAATCCCGTGCAAAAATTCCGCTGCACCCGGTACGGCGACGTTAT  >  NC_000913/699789‑700277                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                 aCGGGCTGGGCTTACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGC‑ACACAACGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATGa                                                                                                                                                                                                                                                 <  1:119428/251‑1 (MQ=255)        gggtttatctacatgggaCGGTTTGCGCCCTGAGATTTTCTCAATCACTGC‑ACACAACGCTCTGTAAGCGGGATAAAAACCGCGCCCGTGGTTGTCCGGATTGGTGGCGATAAATCGTGCATCGTTAGCTACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGTCATCAATCACGAAATCAGGGTTTACATCGGTAAGAGTCAATCCGGCTTTGTTCAGTTCATAAAGCAGtt                                                                                                                                                                                                                                          <  2:197401/234‑2 (MQ=255)          gctTACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGC‑ACACAACGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTgcg                                                                                                                                                                                                                                        <  1:288171/251‑1 (MQ=255)           ctTACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGC‑ACACAACGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGAAGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCGTTGTACAGATCATAAATCAGTgcgc                                                                                                                                                                                                                                       >  2:227862/1‑251 (MQ=255)              aCCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGC‑ACACAACGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCtt                                                                                                                                                                                                                                    >  1:10018/1‑251 (MQ=255)              aCCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGC‑ACACAACGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCtt                                                                                                                                                                                                                                    >  1:271943/1‑251 (MQ=255)                        aaCGGTTTGCGCCCGGAGATTTTCGCAATCCCTGC‑ACACAACGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCcaccac                                                                                                                                                                                                                          <  2:60422/251‑1 (MQ=255)                              ttGCGCCCGAAGATTTTCTCAATCCCTGC‑AAAAAAAGAGACACAAGCGGGATAAAAACCGCGACAGGGGGTGTCCGGCTTGGTGAAGATAAAACGTGCACAGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCAAGTTGTAGGAACGGGTTTCGCCAAAAAACACGAAATAAGGGTTCACAACGGAAAGAAAGAAACCGGCTTAGTACAGGTCAAAAATCAgagctccaccaacaacatcataagg                                                                                                                                                                                                                    >  2:238347/1‑227 (MQ=255)                                gccaccGGAGATATTCCCAATCCC‑GCTACACAGCGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGCTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGATCAAATAGGCTGCTGTATGCATCATTTCCCATTTGTAGGATAGCGTGTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATaagcc                                                                                                                                                                                                                   <  2:50010/246‑5 (MQ=255)                                      ggAGATTTTCTCAATCCCTGC‑ACACAACGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCACCCACCACATACGCTTTCTTGc                                                                                                                                                                                                            >  1:201360/1‑251 (MQ=255)                                                    tcacTGC‑ACACAACGCGCCACAAGCGGGATAAAAACAGCGCCCGTGGGTGTCCGGATTGATGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCCTGTCCCAGTTGTAGGTACGAGTTTCGCCAACAATCACGAATTCCGGGTTCACATCGGTAATAGTGAAACCGGCTTTTTACAGTTAATAAATCAGTGCGCCATCGCCCACCAAATAAGTTTTCTTGCCTTCCTgcggaca                                                                                                                                                                                               >  2:29295/4‑244 (MQ=255)                                                       cTGC‑ACACAACGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACACCATACGCTTTCTTGCCTTCCTGGCGACGCAgg                                                                                                                                                                                           >  2:89465/1‑251 (MQ=255)                                                                aaCGCGCCACAAGCGGGATAAAAACCGCGACCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGACTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGA‑AATCGGc                                                                                                                                                                                  >  2:45859/1‑251 (MQ=255)                                                                          aaacGGGATCAAAACCGCGCCCGTGGGTGTTCGGATTGGTGGCCATAAAACGTGCACCGTTAGCGACGAAATAGTCTGCTTTATGCATAATGTCCCAGGTGTAGGAACGAGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGGGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGA‑AATCGGCAGTCGCCATc                                                                                                                                                                        <  2:252541/248‑1 (MQ=255)                                                                                    aaaaCCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCGCATACGCTTTCTTGCCTTCCTGGCGACGCAGGA‑AATCGGCAGTCGCCATCGCAGAGGtat                                                                                                                                                              <  1:177970/251‑1 (MQ=255)                                                                                                  gggTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGA‑AATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTc                                                                                                                                                 >  1:85467/1‑250 (MQ=255)                                                                                                     tttCCGGATTGGTGGCGATAAAACGTGCACCGTTACCGACGAAATATGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTCCAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGC‑ACTCGGCAGTCGCCATCGCAGAGGTATACACCACGCTGTCAGGt                                                                                                                                             >  1:58763/3‑251 (MQ=255)                                                                                                      gtCCGGATTGGTGGCGATAAAACGTGCTCCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGTTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTGCTTGCCTTCCTGGCGACGCAGGA‑AATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTa                                                                                                                                            <  2:285509/251‑1 (MQ=255)                                                                                                       tCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGA‑AATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTAc                                                                                                                                           >  1:203395/1‑251 (MQ=255)                                                                                                        ccGGATTGGTGGCGATAAAACGTGCACCGTGAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCATGA‑AATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTAc                                                                                                                                           <  1:12779/250‑1 (MQ=255)                                                                                                        ccGGATTGGGGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTGCTTGCCTGCCTGGCGACGCAGGA‑AATCGGCAGTCGCCATCGCTGAGGTATAAAACACGCTGTCAGGTAc                                                                                                                                           <  1:203489/250‑1 (MQ=255)                                                                                                          ggATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGA‑AATCGGCAGTCGCCATCGCAGAGGAATAAAACACGCTGTCAGGTACATc                                                                                                                                        >  1:218995/1‑251 (MQ=255)                                                                                                           gATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGA‑AATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCg                                                                                                                                       <  2:180516/251‑1 (MQ=255)                                                                                                               ggtggCGATAAAACATGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGAAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTATCGTGCCTTCCTGGCGACG‑ACGACAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAAGTACATAGacac                                                                                                                                   >  2:25385/1‑251 (MQ=255)                                                                                                                   gCGATAAAACGTGCACCGTTAGAGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGAAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGA‑AATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTATCAGGTACATCGACACCTGc                                                                                                                               >  1:224679/1‑251 (MQ=255)                                                                                                                                         cgacgaAATAGGCTGCTTTATGCATCATGTCCCAGTTTTAGGAACGCGTTTCGCCAACAATCACGAAAGCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGAA‑AATCGACAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGGACATCGACACCTGCGGTGGCAAAGCGGttatcc                                                                                                            >  1:12047/1‑244 (MQ=255)                                                                                                                                                             tGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGA‑AATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAACGCGGTTCGCCAGATCTTGCCCAGTCTGCGAAgg                                                                                     >  1:288847/1‑251 (MQ=255)                                                                                                                                                                         cAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGA‑AATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGATCTTGCCCAGTCTGCGAAGGATAGTTGGTCAg                                                                         <  1:45859/251‑1 (MQ=255)                                                                                                                                                                                       cgTTTCGCCAACAATCACGAATACAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCGTGCCTTCCTGGCGACGCAGGA‑AATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGATCTTGCCCAGTCTGCGAAGGATAGTTGGTCAGCAACACCAGCGGCa                                                           <  2:130129/251‑1 (MQ=255)                                                                                                                                                                                          ttCGTCAACAATCTCGAAAGCATTGTTCACATCGGTAATAGTGAAACCGGCTTTGTATAGTTCATAAATCAGTGCGCCTTCTCCCTCCACATACGCTTTCTTGCCTTCCTGGCGTCGCAGGA‑AAGCGGCAGGCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCTCCAGATCTTGCCCAGTCTGCGAAGGATAGTTTTTCAGCAACACCAGCGGCAgg                                                         <  2:289225/250‑1 (MQ=255)                                                                                                                                                                                                              cAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGA‑AATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGGTCGCCAGATCTTGCCCAGTCTGCGAAGGATAGTTGGTCAGCAACACCAGCGGCAGGCATTTACCCATAATCCCGTGc                                    >  2:55604/1‑251 (MQ=255)                                                                                                                                                                                                                gggTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGA‑AATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGATCTTGCCCAGTCTGCGAAGGATAGTTGGTCAGCAACACCAGCGGCAGGCCTTTATCCATAATCCCGTGCa                                   >  1:10155/1‑250 (MQ=255)                                                                                                                                                                                                                     cacaTCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGA‑AATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGATCTTGCCCAGTCTGCGAAGGATAGTTGGTCAGCAACACCAGCGGCAGGCCTTTATCCATAATCCCGTGCAAAAAtt                             >  1:69036/1‑251 (MQ=255)                                                                                                                                                                                                                        aTCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGA‑AATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGATCTTGCCCAGTCTGCGAAGGATAGTTGGTCAGCAACACCAGCGGCAGGCCTTTATCCATAATCCCGTGCAAAAATTCCg                          <  1:89465/251‑1 (MQ=255)                                                                                                                                                                                                                             tAATAGTGAAACCGGCTTTGTACAGTTTATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGA‑AATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGATCTTGCCCAGTCTGCGAAGGATAGTTGGTCAGCAACACCAGCGGCAGGCCTTTATCCATAATCCCGTGCAAAAATTCCGCTGc                      >  1:86255/1‑250 (MQ=255)                                                                                                                                                                                                                             tAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGCCGCAGGA‑AATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGATCTTGCCCAGTCTGCGAAGGATAGTTGGTCAGCAACACCAGCGGCAGGCCTTTATCCATAATCCCGTGCAAAAATTCCGCTGc                      >  1:274539/1‑250 (MQ=255)                                                                                                                                                                                                                               aTAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGA‑AATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGATCTTGCCCAGTCTGCGAAGGATAGTTGGTCAGCAACACCAGCGGCAGGCCTTTATCCATAATCCCGTGCCAAAATTCCGCTGCAcc                   >  1:139079/1‑251 (MQ=255)                                                                                                                                                                                                                                                       tCATGAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGA‑AATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGATCTTGCCCAGTCTGCGAAGGATAGTTGGTCAGCAACACCAGCGGCAGGCCTTTATCCATAATCCCGTGCAAAAATTCCGCTGCACCCGGTACGGCGACGTTAt  >  1:37745/1‑244 (MQ=255)                                                                                                                                                                                                                                                       tCATGAATCAGGGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGA‑AATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGATCTTGCCCAGTCTGCGAAGGATAGTTGGTCAGCAACACCAGCGGCAGGCCTTTATCCATAATCCCGTGCAAAAATTCCGCTGCACCCGGTACGGCGACGTTAt  <  2:37745/244‑1 (MQ=255)                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                 ACGGGCTGGGCTTACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGC‑ACACAACGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATGAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGA‑AATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGATCTTGCCCAGTCTGCGAAGGATAGTTGGTCAGCAACACCAGCGGCAGGCCTTTATCCATAATCCCGTGCAAAAATTCCGCTGCACCCGGTACGGCGACGTTAT  >  NC_000913/699789‑700277

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