Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,620,427	G→A	A275A (GCG→GCA)	deoB →	phosphopentomutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,620,427	0	G	A	98.0%	149.7 / ‑4.2	51	A275A (GCG→GCA)	deoB	phosphopentomutase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (33/17);  total (34/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.76e-01

GCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGCGGACATCTACGCCAACTGCGGTATCACCAAA‑AAAGTGAAAGCGACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTAC  >  NC_000913/4620191‑4620672                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                       gCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGGCGTGTTATCGCCCGTCCGTTTATCGGTGACACAGCCGGTAACTTCCAGCGTACCGGTAAACGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGACAACTGCGGTATCACCAAG‑AAAGGTAAAGCAACTGGcacggaagc                                                                                                                                                                                                                                          >  2:155146/1‑243 (MQ=255)       gAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTACGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGtt                                                                                                                                                                                                                                    >  1:447211/1‑251 (MQ=255)          gagcggaacaagggaggcTAAAAAATCGGTCGGGTTATCGCGAGGCAGTGTAACGGTGAAAAAGACGGTAACGTCCAGCGTACCGGTAACCGGAACGACCGGGCGGTTGAGCCGCCAGCACCGACCGGGCTGCAGAAACTGGGTGAGGAAAAACACGGTCAGGGGGTTTCTGTCGGTAAAAATGCAGACATCTACGCCAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGa                                                                                                                                                                                                                                 <  2:163428/236‑1 (MQ=255)              tGACCAACGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGcc                                                                                                                                                                                                                             >  1:86622/1‑251 (MQ=255)                   aaCGGCGGCTACAACATCGGTCGTGTTATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGTTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCAt                                                                                                                                                                                                                        <  1:46538/251‑1 (MQ=255)                                                cGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGCGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAac                                                                                                                                                                                            >  1:25325/1‑250 (MQ=255)                                                       ccGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGt                                                                                                                                                                                    >  1:61734/1‑251 (MQ=255)                                                       ccGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGt                                                                                                                                                                                    >  1:33991/1‑251 (MQ=255)                                                                  tgACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGACTTCACCA‑ACt                                                                                                                                                                        >  2:98955/2‑251 (MQ=255)                                                                           cGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACt                                                                                                                                                               >  1:241862/1‑251 (MQ=255)                                                                               aaCTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGTTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCg                                                                                                                                                            <  1:260362/250‑1 (MQ=255)                                                                               aaCTTCCAGCGTACCGGTAACCGGCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCg                                                                                                                                                            <  2:233418/250‑1 (MQ=255)                                                                                aCTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAG‑AAAGTGACAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAACGAGATGAAAGAAGCGGGAGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGAc                                                                                                                                                          >  2:192085/1‑251 (MQ=255)                                                                                    ccAGCGAACCGGTAACCGTCACGACATGGCTGTTGTGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGGAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAG‑CAAGTGAAAGCAGCTGGCCGGGACGCGCTGTTTGAGGCCACCATCAAAGAGATGAACGAAGAGGGAGATAACACCATCGGCTTCACCA‑ACTTCGTTGACTTCGActttt                                                                                                                                                      >  2:183883/1‑248 (MQ=255)                                                                                              ggTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACCTCTACGCCAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGTCGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCC‑ACTTCGTTGACGTCGACTCTTCCTgggccc                                                                                                                                             >  1:268854/1‑247 (MQ=255)                                                                                              ggTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACAATCGTCTTCACCA‑ACTTCGTTGACTTCTACTCTTCCTGGGGcc                                                                                                                                             >  1:264961/1‑250 (MQ=255)                                                                                                                  tGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGtt                                                                                                                        >  2:134740/1‑251 (MQ=255)                                                                                                                              aGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTg                                                                                                                                                                                                                                  <  2:213141/134‑1 (MQ=255)                                                                                                                              aGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTg                                                                                                                                                                                                                                  >  1:213141/1‑134 (MQ=255)                                                                                                                               gCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCt                                                                                                           >  1:130684/1‑251 (MQ=255)                                                                                                                                 accgaccgTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAATCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCAACGTCGCGACGTCGCCGGTTATGCTGCGGGTCTgg                                                                                                         <  1:415435/251‑1 (MQ=255)                                                                                                                                       cgTGCTGCAGAAACTGGTTGTTGACACACAAGGTCACGTGGTTTATGTCGGTAAAATTGCCGGCATCTACGCAAACTGCGGTATCACAAGG‑AAAGCGAAAGCAACTGGCCTGGACGCGTTGTTTGACGCCACCATCAAAGAGATGAAAGAGGCGGGTGCTAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTCTGCTGCGGGTATGGAGCTGt                                                                                                   >  2:301998/1‑251 (MQ=255)                                                                                                                                          gctgaAGAAAATGGTTGAAGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATGGCAGACATCGACGCCAACTGCGGTAACACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTGCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCg                                                                                                <  2:373389/251‑1 (MQ=255)                                                                                                                                            tgcAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGAc                                                                                              >  1:320015/1‑251 (MQ=255)                                                                                                                                               aGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGAccgc                                                                                           >  1:492654/1‑251 (MQ=255)                                                                                                                                               aGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGAccgc                                                                                           >  1:370290/1‑251 (MQ=255)                                                                                                                                                             gaAAAACTCGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAG‑AAAGGGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCt                                                                             <  1:459015/251‑1 (MQ=255)                                                                                                                                                             gaAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGc                                                                              <  2:307761/250‑1 (MQ=255)                                                                                                                                                               aaaaCACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCtt                                                                            >  2:457614/1‑249 (MQ=255)                                                                                                                                                                 aacacgttcAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATg                                                                         <  1:432244/243‑1 (MQ=255)                                                                                                                                                                  acacGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCtgttgt                                                                        >  2:157239/1‑251 (MQ=255)                                                                                                                                                                   cacGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGActc                                                                                                                                                        <  1:14200/170‑1 (MQ=255)                                                                                                                                                                   cacGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGActc                                                                                                                                                        >  2:14200/1‑170 (MQ=255)                                                                                                                                                                           ggtggtTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCtgctgc                                                               >  1:351695/1‑251 (MQ=255)                                                                                                                                                                           ggtggatcatgacGGTAAAATTGCAGACATCTCCGACAACTGCGGTATAACCAAG‑ACAGACAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCCAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTC‑CCACACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTTTGGGGCTGTTCGACCGCCGTCTGCCGGCGCTGATGTCTCtgctgc                                                               >  2:493516/13‑251 (MQ=255)                                                                                                                                                                                    tGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGGGCTGTTCGACCGCCTTCTGCCGTGGCTGATGTCTCTTCTGCGCGagtacg                                                      >  2:348502/1‑246 (MQ=255)                                                                                                                                                                                       cGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTctgct                                                                 >  1:361863/1‑237 (MQ=255)                                                                                                                                                                                       cGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTctgct                                                                 <  2:361863/237‑1 (MQ=255)                                                                                                                                                                                        ggTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACAt                                                  >  1:319837/1‑251 (MQ=255)                                                                                                                                                                                        ggTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACAt                                                  >  2:406045/1‑251 (MQ=255)                                                                                                                                                                                            aaaTTGCAGACATCTACGCCAACTGCGGTATCACCACG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCCAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTAGGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTg                                              >  2:73568/1‑251 (MQ=255)                                                                                                                                                                                                 gcagACATCTACGCCAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCt                                         <  1:87975/248‑1 (MQ=255)                                                                                                                                                                                                    gACAGCTACGCCAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACCCCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCAc                                      <  2:319837/251‑1 (MQ=255)                                                                                                                                                                                                         cTACGCCAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCActgctg                                 >  1:373871/1‑251 (MQ=255)                                                                                                                                                                                                           aCGCCAACTGCGGTATCACCAAA‑AAAGTGAAAGCGACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGAc                               >  1:145314/1‑251 (MQ=255)                                                                                                                                                                                                              ccAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCCCTGCTGAcccc                            >  2:390609/1‑249 (MQ=255)                                                                                                                                                                                                               cAACTGCGGTATCACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACCTCCTGCTCCTCACTGCTGACCACg                           >  2:401039/1‑251 (MQ=255)                                                                                                                                                                                                                  acgtcgttataCCAAAGAAAGTTAAAGCAACGTGCCTGGACGCGCTTTTTGACGCCACCATCAAAGAGAGGAAAGAGGCGGGTGATATCACCATCGTCTTCTCCA‑ACTGCTTTGACTTCGACGCTTCCTGGGGCCACCGGCGCGCCGTCGCCGGTTATGCTGCGGGTCTGGGGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGt                         <  2:268854/241‑1 (MQ=255)                                                                                                                                                                                                                           cACCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCg                <  2:373871/250‑1 (MQ=255)                                                                                                                                                                                                                            aCCAAG‑AAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGAc              <  1:73568/251‑1 (MQ=255)                                                                                                                                                                                                                                         aaaGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTAc  >  1:295513/1‑251 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                       GCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGCGGACATCTACGCCAACTGCGGTATCACCAAA‑AAAGTGAAAGCGACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCGGGTGATAACACCATCGTCTTCACCA‑ACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTAC  >  NC_000913/4620191‑4620672

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