Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,102	T→C	G123G (GGT→GGC)	thrC →	threonine synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,102	0	T	C	98.1%	151.2 / NA	54	G123G (GGT→GGC)	thrC	threonine synthase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  major base C (23/29);  minor base A (1/0);  total (25/29)
Fisher's exact test for biased strand distribution p-value = 4.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

GACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCACAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGTGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTT‑ACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAACTGAAAGTGGCG  >  NC_000913/3853‑4348                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                        gactgaAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGTg                                                                                                                                                                                                                                                        >  1:69375/1‑251 (MQ=255)              gagaTGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGa                                                                                                                                                                                                                                            <  1:40056/250‑1 (MQ=255)                     tGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTg                                                                                                                                                                                                                                    >  1:290957/1‑251 (MQ=255)                                       cccGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGGCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGAt                                                                                                                                                                                                                  <  2:320473/251‑1 (MQ=255)                                          gCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCGTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGGCCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCGCCGGTGATAcc                                                                                                                                                                                                               <  2:375503/251‑1 (MQ=255)                                            aGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCTCCAATGTTGACAGCGATGTCGGTTGACTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTTCGGGCGATAAGCCAGTGACCATTCTAACCGCGACCTCCGGTGATACCgg                                                                                                                                                                                                             >  1:457561/1‑251 (MQ=255)                                                 gAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTGCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCgg                                                                                                                                                                                                        <  2:67096/251‑1 (MQ=255)                                                  aaGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGACGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGc                                                                                                                                                                                                       >  1:585389/1‑251 (MQ=255)                                                             ggtgtgtATTGGTGATGAAGTCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTGTGCCTTCCCGGTGCCTGTCGCCAATGTGGAAAGCGATGTCGGTTGTCTTGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATg                                                                                                                                                                                            <  2:204157/245‑1 (MQ=255)                                                                  ttATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTTCGCTCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCTGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGTCCTCCGGTGATACCTGAGCGGCAGTGGCTCATGCttt                                                                                                                                                                                        <  1:551559/250‑1 (MQ=255)                                                                       ggTGATGAAAACCCGCAGCAAAGACTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGTCAATGTTGAAAGCGATGTCGGTTGTCTGGTATTGTTCCACTGGCCAACGCTGGCACTTAAAGATTTCGGCGGTCGCTTTATGGCACAAAGGCTGACCCATAGTGCGGGCGAGAAGCCAGTGACCATTCTGACCTCGACCTCCGGTGATACCGGAGCGCCAGTGGCGCAGGCTTTCTACg                                                                                                                                                                                   <  2:145025/250‑1 (MQ=255)                                                                          gatgaAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCGCCGGTGATTCCGGAGCGGCAGTGGCTCATGCTTTCTACGGttt                                                                                                                                                                               <  2:62878/251‑1 (MQ=255)                                                                           atgaAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTTCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGtttt                                                                                                                                                                              >  1:277547/1‑250 (MQ=255)                                                                                    cgcaGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCTCGGCTCCGGTCGCCAATGGTGAAAGCGATGTCGGTTGTCTGGAATTGGTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCGTGATGGCACAAATGCTGACCCAGATGGCGGGCGATAAGTCAGTGACCATTCTGACCGCGGCCTGCGGTGATACCTGAGCGGGATTGGCTGATGCTTTTGCCGGTTT‑TCCGAAtgtg                                                                                                                                                                    <  2:583889/249‑1 (MQ=255)                                                                                      caGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCCTGCTTTCTACGGTTT‑ACCGAATGTGa                                                                                                                                                                   >  1:236895/1‑250 (MQ=255)                                                                                              ccTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGTCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCGTTCTACGGTTT‑ACCGAAAATGAAAGTGGttt                                                                                                                                                          >  1:7764/1‑250 (MQ=255)                                                                                              ccTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGGCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTGCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCTGTTACCATTCTGACCGCTTCCTCCGGTGATACCGGAGCGGCAGTGGCTCAGGCTTTCTACGGTTT‑ACCGAATGTGAAAGTGGTTa                                                                                                                                                          <  1:259619/251‑1 (MQ=255)                                                                                                 ggAAGAGCGCGTGCGCGCGGAGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTTCGGGCGATAAGCCAGTGACCATTCTGACCACGACCTCCGGTGATACCGGAGCGGCAGTGGCTCCTGCTTTCTACGGTTT‑TCCGGATGTGGAAGTGGTTAttc                                                                                                                                                       >  2:171426/1‑249 (MQ=255)                                                                                                   aaGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTT‑ACCGAATGTGAAAGTGGTTATCctc                                                                                                                                                     <  1:435501/251‑1 (MQ=255)                                                                                                     gagCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTT‑ACCGAATGTGAAAGTGGTTATCCTCTa                                                                                                                                                   <  1:83134/251‑1 (MQ=255)                                                                                                                                ccGGCGCCGGTCGCTAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGGTCCCCGGGCCAACTCTGGCATTTTAAGATTTCGTCGGTCGCTTTATGGCACAAATGCTGGCCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGGCCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTGTCTACGGGTC‑ACCGAATGTTAAAGTGGTTATACTCTATCCACTAGGCAAAATCGGTCCACTGCa                                                                                                                        <  2:585665/251‑1 (MQ=255)                                                                                                                                  ggCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTT‑ACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAg                                                                                                                      <  1:20175/251‑1 (MQ=255)                                                                                                                                        ggTCGCCAATGTTGAACGCGATGGCGGTTGTCTGGAATTGTTACACGAGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAACTTCTGACCCATATTGTGGGCGTTTAGCCCGTGACCATTCTGACCGCGACCTCCGGTTGTCCCAGAGCGGAAGTGGCTCATGCTTTCTACGGTTTTACTGA‑TGTGACAGTGGTTTTCCTCTCTCAACGAGGCACAATCCGGCCACTGCTAGAAAAAc                                                                                                                >  2:104435/1‑251 (MQ=255)                                                                                                                                        ggTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTT‑ACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAAc                                                                                                                <  1:314670/251‑1 (MQ=255)                                                                                                                                                  gTTGAAAGCGATGTCGGTTGGCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGTGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTT‑ACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTAc                                                                                                      <  2:297698/251‑1 (MQ=255)                                                                                                                                                   ttGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGATGACCCATATTGCGGGCGATAAGCCAGTGAACATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAATGGCTCATGCGTTCTAAGGGTT‑ACCGAATGTGAAAATGGTTATCCTCTATCCACGAGGAAAAATCAGTCCACTTCAAGAAAAAATGTTCTGTAc                                                                                                      >  2:529948/1‑250 (MQ=255)                                                                                                                                                     gAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTAGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGAACTCCGGTGATACCGGAGCGGCAGTGGCTAATGCTTTCTACGGTTT‑ACCGAATGTGAAAGTGGTTATCCTCTATCCACGCGGCAAAATCAGTCCACTGCCAGAAAAACTGTTCTGTACAtt                                                                                                   >  1:604445/1‑251 (MQ=255)                                                                                                                                                      aaaGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCTCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTT‑ACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGGCCACTGCAAGAAAAACTCTTCTGTACATTg                                                                                                  <  2:602003/251‑1 (MQ=255)                                                                                                                                                        aGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTT‑ACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTggg                                                                                                <  1:474347/251‑1 (MQ=255)                                                                                                                                                           gATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTT‑ACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAACAACTGTTCTGTACATTGggcgg                                                                                             >  1:623721/1‑251 (MQ=255)                                                                                                                                                              gTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTT‑ACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCa                                                                                           >  1:491622/1‑250 (MQ=255)                                                                                                                                                                           aaTTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGATTTCTACGGTTT‑ACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGtt                                                                             >  1:526402/1‑251 (MQ=255)                                                                                                                                                                              tggtCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATGCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTT‑ACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGc                                                                           <  1:230302/247‑1 (MQ=255)                                                                                                                                                                              tGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTT‑ACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGc                                                                           <  2:216169/250‑1 (MQ=255)                                                                                                                                                                                 tCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTT‑ACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGAGGCAATATCGAAACTGTTGCCAt                                                                        <  1:352818/250‑1 (MQ=255)                                                                                                                                                                                      gggCCCACGCTGGAATTTAAAGATTTCGGCGGTCGCTTTATGGCACAACTGGTGACCCATATTGCGGGCGATAAGCCAGTGACCCTTCTGAACCCGACCTCCGGTGCTACCGGAGCGGCAGCGGCCCACGCTTTCTACGGTTT‑AACGAATGCGAAAGGGGTTATCCACTATCCACGATGCAAACTCAGGCAACGGACAGAACAACTGCTATGTTCATTGGGCGGCAATAACGacaatgtccccctcggcgg                                                                  >  2:371528/1‑233 (MQ=255)                                                                                                                                                                                       ggCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTT‑ACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCCTCGACGGc                                                                 >  2:232986/1‑251 (MQ=255)                                                                                                                                                                                           aaCGCTGGCATTTAAAGATTTCGGCGGGCGCTTTATGGCACAAAAGCTGACACATATTGCGGGCGATAAGCAAGTGACCATTCTGACCGCAACCTCCGGTGATACCGGAGCGGCAGCGGCTCATGCTTTCTACGGTTT‑ACCGAATGTGCAAGTGGTTATCCTCTATCCAAGAGGCAAAATAAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCAGCAATATCGCAACGATTGCCATCGACGGCGAtt                                                             >  1:363932/1‑251 (MQ=255)                                                                                                                                                                                               cTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTT‑ACCGAATGTGAAAGTGGTTATCCTCTCTCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCg                                                                >  1:388155/1‑244 (MQ=255)                                                                                                                                                                                               cTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTT‑ACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCg                                                                <  2:388155/244‑1 (MQ=255)                                                                                                                                                                                                    cgttAAAGATTTCGGTGGTGGCGTTATGGGACAAATGCTGACGCAGATTGCGGGCGATAAGCCTGTGACCATTCTGACCTCCACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTG‑ACCGAATGTGAAAGTGGTTACGCTCTATTCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATAGAAACTGTTGCCATCGACGGCGATTTCGATGCCt                                                    <  2:356223/249‑1 (MQ=255)                                                                                                                                                                                                    aTTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTT‑ACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCt                                                    <  1:436067/251‑1 (MQ=255)                                                                                                                                                                                                          aGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGAGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTT‑ACCGAAAGTGAAAGTGGGTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGACAAACTGTTCTGTACATTGGGCGGCAATATCGACACTGTTGCCATCGACGGCGATTTCGAGGCCTGTCAgg                                              >  2:521993/1‑251 (MQ=255)                                                                                                                                                                                                            aTTTCGGCGGTCGCTTTATGGCAAAAATGCTGACCCATATTGCGGGCGGTTAGCCAGTGACCATTTTGACAGCGACCTCCGGTGATACCGGAGCGGCAGGGGGTAATGCTTTCTACGGTTT‑ACCGAATGTGAAAGGGGGTATCCTCTATCCACGAGGCCAAATAAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAAATGTTGACATAGATGGCGATTTCGGTGCCTGTCAGgcg                                            >  2:335938/1‑251 (MQ=255)                                                                                                                                                                                                                cggaggTCGCTTTATGACCAAAATGCAGACCCATATTTCGGGCGATAAACCAGTGACCCTTCTGAACGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTT‑ACCGAAGGTGAAAATGGATATCCTCTATCCAAGAGGAAAAAACAGTCCCATGAACAAAAAACTGTTCTGTACATTGGGCGGCAAGATCGAAACTTTTCCCCTCGACAGCGATTTAGATGCCTGTAAGGGGCTgg                                        >  2:436027/5‑251 (MQ=255)                                                                                                                                                                                                                       cGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTT‑ACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGc                                  >  2:414830/1‑250 (MQ=255)                                                                                                                                                                                                                           ttATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTT‑ACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGc                              <  1:593507/250‑1 (MQ=255)                                                                                                                                                                                                                            tATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTT‑ACCGAATGTGAAAGTGGTTATACTATATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGt                            >  1:277319/1‑251 (MQ=255)                                                                                                                                                                                                                               ggCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTT‑ACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCTGCAATATCGAAACTGTTGCCATCGACGGCGATGTCGATGCCTGGCTGGCGCTGGTGAAGCAGGCGttt                          <  2:316203/250‑1 (MQ=255)                                                                                                                                                                                                                                gCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTT‑ACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTtga                        <  1:288301/251‑1 (MQ=255)                                                                                                                                                                                                                                  acaAATGCTGACCCATATTGCGGGCGATAAGCAAGTGACCATTCTGACCGCCACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTT‑ACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCCAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTtgctg                      >  1:161924/1‑248 (MQ=255)                                                                                                                                                                                                                                                  atTGCGGGCGATAAGCGAGTGACCATTCGGACCGCGACGTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTT‑ACCGAATGTGAAAGTGGTTATCCGCGATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAGCTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCGGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGt      <  2:291432/251‑1 (MQ=255)                                                                                                                                                                                                                                                  atTGCGGGCGATAAGCCTGTGACCATTCTTACCGCGACCTCCGGTGATATCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTT‑ACCGAATGTGAAAGTTGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGAGTTCTATGCCTATCAGGCGCTGGTGAAGCAGGCGTTtgttgactaataattaatat       <  2:377353/250‑15 (MQ=255)                                                                                                                                                                                                                                                      cGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTT‑ACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAAATGTTCTGTACATTGGGCGGCAATATCGACACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGgcg  >  2:457014/1‑251 (MQ=255)                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                        GACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCACAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGTGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTT‑ACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAACTGAAAGTGGCG  >  NC_000913/3853‑4348

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