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	95.3%	120.1 / ‑5.0	44	S522G (AGC→GGC)	acnA	aconitate hydratase 1
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (17/24);  total (18/26)
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

TGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAA‑AGCGATTTAACCGTCGGTGCGGTGCTGTCCGG‑CAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCATCGGCACAAGAAATTGCCCGTGCGGTAGAACAAGTCTCC  >  NC_000913/1337161‑1337636                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                    tGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGGTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAACCGGGAATCAAAAAA‑GGCGATTTAACCGTcggt                                                                                                                                                                                                                                    >  2:407879/1‑251 (MQ=255)            ggctggttGGCGAAAAAAGCCGTAACTCTTGGCCGCAAGCTGCAACCATGGGTCAAAGCGTCGCTGGCAGCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTTGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGTAACGGCAAGCAAAAAA‑GGCGATTGAACCGTCTGTGCGGTGCTGTc                                                                                                                                                                                                                         <  2:480893/244‑1 (MQ=255)                          aaaaGCCGTAACTCTGGTCATCAAGCGGCAAACATGGGTCAAAGCGTCGCCGGCAGCGGGTGCGAAAGTCGGTTCTGATTATCTGGCAAAAGCGAAACTTACACCGTATCTCGACGAACTGGGGTTTAGCCTTGTGGGATACTGTTTTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAA‑GGCGATTTAACCGTCGGTGCGGTGCTGTCCGG‑CAACCGTAACt                                                                                                                                                                                                          <  1:253177/251‑1 (MQ=255)                          aaaaGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAA‑GGCGATTTAACCGTCGGTGCGGTGCTGTCCGG‑CAACCGTAACt                                                                                                                                                                                                          >  1:100354/1‑251 (MQ=255)                                        tGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACACCCTGTATGGGTAACTCTGGGCCGCGGCCCGATCTTTCCGAAACGGCAATCAAAAAAGAG‑GATTAAACCGCCGGTGCGGGGCTGCCCGG‑AAACGGTAACTTTGAAGCCCGAatc                                                                                                                                                                                            >  1:257855/1‑251 (MQ=255)                                         gggcaccAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAA‑GGCGATTTAACCGTCGGTGCGGTGCTGTCCGG‑CAACCGTAACTTTGAAGGCCGTatc                                                                                                                                                                                            <  1:1619/244‑1 (MQ=255)                                         gggCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGTGTTCGAAAGGCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGGCGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAA‑GGCGATTTAACCGTCGGTGCGGTGCTGTCCGG‑CATCCGTAACTTTGAAGGCCGTatcc                                                                                                                                                                                           <  2:502465/251‑1 (MQ=255)                                          ggCCTCAAGCGGCAATCGTGGGTCCAAGCGTCGCTGAGACCGGGTTCGACAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGCGACCCTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGACACGGCAATCAAAAAA‑GGCGATTTAACCGTCGGTGCGGTGCTGTCCGG‑CAACCGTAACTTTGAAGGCCGTatcc                                                                                                                                                                                           >  2:243209/1‑250 (MQ=255)                                           gcctcaggcgtccacggtaagaaaaaGGGTCGCTGGAACCGGGCTCGCAAGTAGCTTCTGATTATATGGCACAAGCGAAACTGCCACCGTATCTCGACGAACTGGGGTTTAACATTGTGGGATACGGTTGTACAACCTGGATTGGTAGCTGTGGTCCGCTGCCGGATCCTATCGAAACGGCAATCGGCAAA‑GGAGAGTTAACCGTCGGCGCGGGGCGGACCAG‑CAACCGTACCTTTGAAGGCCGTatccat                                                                                                                                                                                         >  2:197014/24‑251 (MQ=255)                                                aaGCGGCAACCATGGGTCAAAGCGTCGCTGGCTCCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAA‑GGCGATTTAACCGTCGGTGCGGTGCTGTCCGG‑CAACCGTAACTTTGAAGGCCGGATCCATCCGCt                                                                                                                                                                                    <  2:307427/251‑1 (MQ=255)                                                 aGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAA‑GGCGATTTAACCGTCGGTGCGGTGCTGTCCGG‑CAACCGTAACTTTGAAGGCCGTATCCATCCGCTg                                                                                                                                                                                   <  2:327717/251‑1 (MQ=255)                                                    ggCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAA‑GGCGATTTAACCg                                                                                                                                                                                                                                         <  2:81190/195‑1 (MQ=255)                                                    ggCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAA‑GGCGATTTAACCg                                                                                                                                                                                                                                         >  1:81190/1‑195 (MQ=255)                                                          cATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAA‑GGCGATTTAACCGTCGGTGCGGTGCTGTCCGG‑CAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAaact                                                                                                                                                                          >  1:31439/1‑251 (MQ=255)                                                                 cAAAGCGTCGCTGGCGCCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGGACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAA‑GGCGATTTAACCGTCGGTGCGGTGCTGTCCGG‑CAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAACCTAAccggc                                                                                                                                                                   >  2:23035/1‑251 (MQ=255)                                                                             ggCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAA‑GGCGATTTAACCGTCGGTGCGGTGCTGTCCGG‑CAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTcgcc                                                                                                                                                         <  1:429398/249‑1 (MQ=255)                                                                                  cGGGGTCGAAAGTCGGTTCTGATTCTCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGATGTACCACCTGTATTGGTAACTCTGGGACGCTGCACGATCCTATAGAAACCGCAATCAAAAAA‑GGAGATTTAACCGTCGGAGCGATGATGTCAGGACAAC‑GTAAATTTGACGGACGAATACATACGCTGGTTAAAACTAACTGGCTGGCCTtgcccacgctg                                                                                                                                                  >  1:292766/1‑244 (MQ=255)                                                                                        cGAAAGTCGTTTCTGCTCATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGTACGGGGGTTTATCCTTATGGGGTACGGTTGTACCACCTGTAATGGTAACTCTTGGCCGCTGCTCGATCCTATCGATACGGCAATCAAAAAA‑AGTGATTTAACCGTCGGTGCTGTGCTTTCCTG‑CAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCtggtgggt                                                                                                                                            <  1:425008/251‑3 (MQ=255)                                                                                         gAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAA‑GGCGATTTAACCGTCGGTGCGGTGCTGTCCGG‑CAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTg                                                                                                                                           >  1:126129/1‑251 (MQ=255)                                                                                         gAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAA‑GGCGATTTAACCGTCGGTGCGGTGCTGTCCGG‑CAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTg                                                                                                                                           <  2:109314/251‑1 (MQ=255)                                                                                               cGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCGCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAA‑GGCGATTTAACCGTCGGTGCGGTGCTGTCCGG‑CAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATg                                                                                                                                     <  2:100354/251‑1 (MQ=255)                                                                                                  ttCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAA‑GGCGATTTAACCGTCGGTGCGGTGCTGTCCGG‑CAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATgcg                                                                                                                                   <  2:173923/250‑1 (MQ=255)                                                                                                                  aaaGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAA‑GGCGATTTAACCGTCGGTGCGGTGCTGTCCGG‑CAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGa                                                                                                                   >  1:17196/1‑250 (MQ=255)                                                                                                                      cGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAA‑GGCGATTTAACCGTCGGTGCGGTGCTGTCCGG‑CAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATc                                                                                                              >  1:247303/1‑251 (MQ=255)                                                                                                                         aaCTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAA‑GGCGATTTAACCGTCGGTGCGGTGCTGTCCGG‑CAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAAc                                                                                                           >  1:181389/1‑251 (MQ=255)                                                                                                                                      tctcggagaaatGGGGTTTAAACTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAA‑GGCGATTTAACCGTCGGTGCGGTGCTGTCCGG‑CAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAg                                                                                               <  2:257855/239‑1 (MQ=255)                                                                                                                                        ccgacgaACTGGGGTTTAACCTTGTGGGATACTGTTGTACCTCCGGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAA‑GGCGATTTAGCCTTCGGTGCGGTGCTGTCCGG‑CAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTTGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCt                                                                                            <  1:320327/250‑1 (MQ=255)                                                                                                                                         cgacgaACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAA‑GGCGATTTAACCGTCGGTGCGGTGCTGTCCGG‑CAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGGATATCAACCTGGCTTCTGAg                                                                                               <  1:202434/247‑1 (MQ=255)                                                                                                                                              aaCTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAA‑GGCGATTTAACCGTCGGTGCGGTGCTGTCCGG‑CAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGc                                                                                      <  1:340974/251‑1 (MQ=255)                                                                                                                                                              tgtgGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAA‑GGCGATTTAACCGTCGGTGCGGTGCTGTCCGG‑CAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAgg                                                                        <  2:297134/249‑1 (MQ=255)                                                                                                                                                                      aCGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTAACGAAACGGCAATCAAAAAA‑TGCGATTTAACCGTCGGTGCGGTGCTGTCCGG‑CAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAAAATCAACCTGGCTTCTGATCCTATCGGCCTTGATCGCAAAGGCGATCCGGtt                                                              <  1:405351/251‑1 (MQ=255)                                                                                                                                                                                caccTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGAAATCAAAAAA‑GGCGATTTAACCGTCGGTGCGGTGCTGTCCGG‑CACCCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGACATAAGAATATCAACCTGGCTTCTGTGCCTATCGGCAATGATCGCACAGGAGATCCGGGTTATTTGaa                                                      >  2:492290/1‑249 (MQ=255)                                                                                                                                                                                 accTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAA‑GGCGATTTAACCGTCGGTGCGGTGCTGTCCGG‑CAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGCAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGaa                                                      >  1:212889/1‑248 (MQ=255)                                                                                                                                                                                  ccTGTAGTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAA‑GGCGATGTAACCGTCGGTGCGGTGCTGTCCGG‑CAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGat                                                  <  2:489885/251‑1 (MQ=255)                                                                                                                                                                                       aTTGGTAACTCTGGGCCGCTGCCCGATCCTATCTAAACGGCAATCAAAAAA‑GGCGATTTAACCGTCGGTGCGGTGCTGTCCGG‑CAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCGATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTg                                             <  2:328239/251‑1 (MQ=255)                                                                                                                                                                                         tGGTATCTCTGGGCCGCTGCCCGATCCTAGCGAAACGGCAATCAAAAAA‑GGCGATTTAACCGTCTGTGCGGTGCTGTCCGG‑CAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTGCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCTGTTTATCTGAAAGATATCTgg                                            <  2:273621/250‑1 (MQ=255)                                                                                                                                                                                             aaCTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAA‑GGCGATTTAACCGTCGGTGCGGTGCTGTCCGG‑CAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGGTCGCAAAGGCGATCCGGTTTATCTGAAAGATATTTGGccctc                                       >  2:214394/1‑248 (MQ=255)                                                                                                                                                                                                   gggCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAA‑GGCGATTTAACCGTCGGTGCGGTGCTGTCCGG‑CAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCATCGGcaca                                 <  1:214394/251‑1 (MQ=255)                                                                                                                                                                                                           gCCCGATCCTATCGAAACGGCAATCAAAAAA‑GGCGATTTAACCGTCGGTGCGGTGCTGTCCGG‑CAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTAGGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCa                                                                                    <  2:83188/192‑1 (MQ=255)                                                                                                                                                                                                           gCCCGATCCTATCGAAACGGCAATCAAAAAA‑GGCGATTTAACCGTCGGTGCGGTGCTGTCCGG‑CAACCGTAACTTTGAAGGCCGTATACATCCGCTGGTTAACACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCa                                                                                    >  1:83188/1‑192 (MQ=255)                                                                                                                                                                                                             ccGATCCTATCGAAACGGCAATCAAAAAA‑GGCGATTTAACCGTCGGTGCGGTGCTGTCCGG‑CAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCATCGGCACAAGAAATTGcc                       >  2:429688/1‑251 (MQ=255)                                                                                                                                                                                                              cGATCCTATCGAAACGGCAATCAAAAAA‑GGCGATTTAACCGTCGGTGCGGTGCTGTCCGG‑CAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCATCGGCACAAGAAATTGccc                      <  2:372133/251‑1 (MQ=255)                                                                                                                                                                                                                           aCGGCAATCAAAAAA‑GGCGATTTAACCGTCGGTGCGGTGCTGTCCGG‑CAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCATCGGCACAAGAAATTGCCCGTGCGGTAGAACa         <  2:361218/251‑1 (MQ=255)                                                                                                                                                                                                                                  acAAAAAA‑GGCGATTTAACCGTCGGTGCGGTTCTGTCCGG‑CAACCGGAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCATCGGCACAAGAAATTGCCCGTGCGTTAGAACAAGTCTcc  <  1:393793/250‑1 (MQ=255)                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                    TGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAA‑AGCGATTTAACCGTCGGTGCGGTGCTGTCCGG‑CAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCATCGGCACAAGAAATTGCCCGTGCGGTAGAACAAGTCTCC  >  NC_000913/1337161‑1337636

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