Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,337,394	A→G	S522G (AGC→GGC)	acnA →	aconitate hydratase 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,337,394	0	A	G	98.4%	204.1 / NA	61	S522G (AGC→GGC)	acnA	aconitate hydratase 1
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (28/32);  minor base T (0/1);  total (28/33)
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.98e-01

AACCCAAGTGTGCTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAAGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCATCGGCACAAGAAATTGCCCGTGCGGTAGAACAAGTCTCCA  >  NC_000913/1337148‑1337637                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    aaCCCAAGTGTGCTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCg                                                                                                                                                                                                                                                  >  1:259575/1‑250 (MQ=255)    ccAAGTGTGCTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGAttt                                                                                                                                                                                                                                              >  1:350343/1‑251 (MQ=255)        gtgtGCTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAAcc                                                                                                                                                                                                                                          >  1:361887/1‑251 (MQ=255)                  ggCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAGGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTcggtgcg                                                                                                                                                                                                                                 <  2:14998/250‑1 (MQ=255)                         ggCTTGCTGGCGAAAAAAGCGGTAACTCTGGGCCTCAAGCGGCAACCGTGGGTCAAATCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGTGGTTTAACCTTGTGGGATACGGTTGGACCACCTGTATGGGTAACTCTGGGCCGCTGCCCGATCCTATCGAATCGGCAATCAAAAAAGGCGATTTAACCGCCGGTGTGGTGCTGTc                                                                                                                                                                                                                         <  2:197926/251‑1 (MQ=255)                           cTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCg                                                                                                                                                                                                                       >  1:348982/1‑251 (MQ=255)                           cTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTAGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGCGCTGTCCg                                                                                                                                                                                                                       >  1:173482/1‑251 (MQ=255)                                 ggCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAAcc                                                                                                                                                                                                                 >  1:196010/1‑251 (MQ=255)                                 ggCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAAcc                                                                                                                                                                                                                 <  2:361887/251‑1 (MQ=255)                                 ggCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCGTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAAcc                                                                                                                                                                                                                 <  1:349765/251‑1 (MQ=255)                                 ggCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGGCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAAcc                                                                                                                                                                                                                 <  2:453908/251‑1 (MQ=255)                                   cGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGc                                                                                                                                                                                                                     >  2:300124/1‑245 (MQ=255)                                   cGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGc                                                                                                                                                                                                                     <  1:300124/245‑1 (MQ=255)                                    gAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGCTTCACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCTGTGCTGTCCGGCAACCGt                                                                                                                                                                                                               >  2:406635/1‑250 (MQ=255)                                           gCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGCTTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTACCCATGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCACTCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTg                                                                                                                                                                                                        >  1:29724/1‑250 (MQ=255)                                                     tGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTatc                                                                                                                                                                                             >  1:56899/1‑251 (MQ=255)                                                             aaGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCt                                                                                                                                                                                     <  2:56899/251‑1 (MQ=255)                                                               gCGGCAACCATGGGTCAAAGCGGCGCTGGCACCGGTTGCGAAAGTCGTTTCTGATTAGCTGGCAAAAGCGAAACTGACACCGTATCTCGCCGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAGCGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTgg                                                                                                                                                                                   <  2:300873/251‑1 (MQ=255)                                                                 ggCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGtt                                                                                                                                                                                 >  1:63764/1‑251 (MQ=255)                                                                    aaCCATTGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTGTCTGATTATCTGGCAAAAGCGAAACTGACGCCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAGCGGCAATCAAAAAAGGCGATTTAACCTTCGGTGCGGTGATGTCCGTCAACCGTAACTTTGAAGGCCGTATCCCTCCGCTGGTTaaa                                                                                                                                                                              <  2:39386/251‑1 (MQ=255)                                                                              cAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGCCACCGTATCTCGACGAACTGGGGTTTAACCGTGGGGGATACGGTGGTACCACCTGTATTGGTAACTCTGGGCCGCTCCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAActgg                                                                                                                                                                     <  2:151659/250‑1 (MQ=255)                                                                                aaGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTGTGAAGGCCGTATCCATCCGCTGGTTAAAACTAActggctg                                                                                                                                                                  <  2:259575/251‑1 (MQ=255)                                                                                         tGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGCTACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAAATTTGAAGGCCGTATCCATCNGCTGGTTAAAACTAACTGGCTGGCCTcgccg                                                                                                                                                         >  2:211846/1‑251 (MQ=255)                                                                                         tGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGAAACTTTGAAGGCCGTATCCATCCGCTGGTTAACACTAACTGGCTGGCCTcgccg                                                                                                                                                         >  2:317176/1‑251 (MQ=255)                                                                                                     cGAAAGTCGTTTCTTGTTATCTGGCAAAAGCGAAAGTGACACAGGATCTCGCCGAACTGGGGTTTAACCGTGTGGGATACGGTTGTACCCCCTGTATTGTTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAGCCGGCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAATGCCGTATCCATCCGCTGGTTAAAACTAATTTGCTGGCCTCGCCGTCTCTGGTGGtt                                                                                                                                             <  2:264368/251‑1 (MQ=255)                                                                                                           tCGTTTATGAATATCTGACAAAAGCGAAACTGACACCGTATCTCTACGAACTGGGGATTAACCATGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCACTGCCCGATCCTATCGAAACGGCAATCAAAAACGGCGATTTAACCGTCGGTGCGGTGATGTCCGGCGACATAAAATTAGAAGACCGAATGCATCCGCAGGTTAAAACTACCtgggtggcacttccggcgctggtgtttgccact                                                                                                                                       >  2:440220/1‑226 (MQ=255)                                                                                                            cGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCTATCCTATCGAAACGGCAATCAAAAAATGCGATTTAACCGTCGGTGCGGTGCTGACGGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTTGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGGGGTTGCCTAt                                                                                                                                       <  2:419076/250‑1 (MQ=255)                                                                                                              tttCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAACAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGCAACTTTGAAGGCCGTATCCATCCGATGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGCGGTGGCATATgcg                                                                                                                                    >  1:291013/1‑251 (MQ=255)                                                                                                               tacggattacttGGCAAAAGCGAAACTGACACCGTATCGCGACGAACTGGGGTTTAACCTTGTGGGAGACGGTTGTACCACCGGTATTGGTATCTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGTCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCGTCCGCTGGTTAAATCTAACTGGCTGGCCTCGCCGCTGCTGGTGGTTGCCTATgcg                                                                                                                                    <  2:310032/239‑1 (MQ=255)                                                                                                                   gATTATCTGGCAAAAGCGAAGCTGACACCGTATCTCGACGAACGGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTAGCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCt                                                                                                                                  <  2:334243/248‑1 (MQ=255)                                                                                                                       aTCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGCCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCggg                                                                                                                            >  1:456770/1‑250 (MQ=255)                                                                                                                                 aGCGAAACTGACGCCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCTATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGATATATGAata                                                                                                                 <  2:245913/251‑1 (MQ=255)                                                                                                                                 aGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAata                                                                                                                 <  1:25659/251‑1 (MQ=255)                                                                                                                                 aGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAat                                                                                                                  >  2:359305/1‑250 (MQ=255)                                                                                                                                          gACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTgg                                                                                                        <  1:42102/251‑1 (MQ=255)                                                                                                                                           acacCGTATCTCGTCGTACTGGGTTTTAACCTTGGTGGCTACTGTTGTACCACCTGTATTGGTATCTCGGGGCCGCTGCCCTATCCTATCGAAACGGAAATCAAAAAAGGCGATTGAACCGTCGGTGCGTTGCGGCCCGGCACCCGTAACTTTGAAGGCCGTATCCATCCGCCGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGCGGTTGCCTATGCGCTGGCGTGAAATATGAATGTCAACCTGGc                                                                                                       <  2:12346/251‑1 (MQ=255)                                                                                                                                           acacCGTATCTCGACTAACTGGTGTTTAACCNTTTGTGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTACCCGTCGGTGCGGTTCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTTGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGc                                                                                                       <  2:237933/251‑1 (MQ=255)                                                                                                                                             acCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCATCCTGGCtt                                                                                                     <  2:356381/251‑1 (MQ=255)                                                                                                                                                 ttcttcGACGAACTGGGGTTTAACCTTGTGGGTTCCGGTTGAACCACCTGTATTGGTATCTCTGGGCCGCCGCCCGATCCTATAGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGAAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATTAACCTGGCTTCt                                                                                                   <  2:312189/245‑1 (MQ=255)                                                                                                                                                                     tAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGc                                                                              <  1:20657/250‑1 (MQ=255)                                                                                                                                                                         cTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTc                                                                                                    <  2:359164/224‑1 (MQ=255)                                                                                                                                                                         cTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTc                                                                                                    >  1:359164/1‑224 (MQ=255)                                                                                                                                                                             ggggATACGTGTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGAt                                                                     <  2:196010/250‑1 (MQ=255)                                                                                                                                                                               ggATACGGTTGTACCACCTGTATTGGTAACTNTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATcc                                                                   <  2:225979/251‑1 (MQ=255)                                                                                                                                                                                      gTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAAACTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTAt                                                            >  1:113511/1‑251 (MQ=255)                                                                                                                                                                                             caccTGGATGGGTAACTCGTGGCCGCGTCCCGATCCTATCGAAACGTCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGCCCGGCAACCGTAACTTTGAAGGCCGTATCCAGCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAg                                                     <  2:437987/251‑1 (MQ=255)                                                                                                                                                                                                  gTATTGGTAACTCTGAGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAACAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGACATATGAATATCAACCTGGCTTCTGAGCCTCTCGGCCAGGATCGCCAAGGCGACCCGGTTTTTTTGAAAGaattc                                                >  2:26002/1‑247 (MQ=255)                                                                                                                                                                                                        gTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCCATCGGCCATGATCGCAAAGNCGATCCGGTTTATTTGAAAGATATCTGGcc                                           >  2:226446/1‑250 (MQ=255)                                                                                                                                                                                                        gTAACTCGGGGCCGCGGCCCGATCCTAGCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCa                                          <  2:348982/251‑1 (MQ=255)                                                                                                                                                                                                           aCTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGGGATCCGGGTTATCTGAAAGATATCTGGCCATCg                                       >  1:55221/1‑251 (MQ=255)                                                                                                                                                                                                            ctctGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCATCgg                                      >  1:155646/1‑251 (MQ=255)                                                                                                                                                                                                                   ccGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGCAATATGAATATCAACCCGGCTTCTGAACCTATCGGCCCTGCTCGCAAAGGCGTCCCGGGTTATCTGAAAGATATCTGGCCCTCGGcccc                                  >  2:35146/1‑245 (MQ=255)                                                                                                                                                                                                                         cccGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCATCGGCACAAGAAATTGc                         <  1:226446/251‑1 (MQ=255)                                                                                                                                                                                                                            gATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTTAGCCTATCGGCCCTGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATTTGGGCATCGGCACAAGAAAATGCCCg                      >  2:358092/1‑251 (MQ=255)                                                                                                                                                                                                                            gATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCAATCCGGTTTATCTGAAAGATATCTGGCCATCGGCACAAGAAATTGCCCg                      <  1:261234/251‑1 (MQ=255)                                                                                                                                                                                                                                       aaCGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAatat                                                                                                                >  1:38592/1‑150 (MQ=255)                                                                                                                                                                                                                                       aaCGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAatat                                                                                                                <  2:38592/150‑1 (MQ=255)                                                                                                                                                                                                                                         cGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCCCAAAGGCGATCCGGTTTATCTGAAAGGTAACTTGCCATCCGCACAAGAAATTGCCCGTGCGGGAGAAcc          >  2:149898/1‑249 (MQ=255)                                                                                                                                                                                                                                         cGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGGAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGGGGTTGCCTATGCGCTGGCGGGAAATATGAATCTCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGGTTATCTGAAAGATATCTGGCCATCGGCACAAGAAATTGCCCGTGCGGTAGAACa          >  1:119875/1‑250 (MQ=255)                                                                                                                                                                                                                                           gCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCATCGGCACAAGAAATTGCCCGTGCGGTAGAACAAGt       <  1:4656/251‑1 (MQ=255)                                                                                                                                                                                                                                                cAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGTGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCATCGGCACAAGAAATTGCCCGTGCGGTAGAACAAGTCTcc   <  1:423339/250‑1 (MQ=255)                                                                                                                                                                                                                                                 aaaaaagGAGATTTAAGAGTCGGTGCGGTGCTGTCCGGCAAACGTAACTTTGAAGGCCGTTTCCATCAGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGCTGCCTATGCGCTGGAGGGAAATATGAATATCAACCCGGCTTATGAGCCGATAGGCCATTATCGCAAAGGGCATCCGGAGTAATTGAAAGATATTTTGCCCTCGGCACAAGGAATTGccccgggcgaagaaaaagtccccc  >  2:265427/1‑229 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    AACCCAAGTGTGCTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAAGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCATCGGCACAAGAAATTGCCCGTGCGGTAGAACAAGTCTCCA  >  NC_000913/1337148‑1337637

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