Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,041,451	G→A	A47T (GCA→ACA)	ygfZ →	folate‑binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,041,451	0	G	A	96.1%	148.1 / ‑1.3	52	A47T (GCA→ACA)	ygfZ	folate‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (24/25);  total (26/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 = 6.30e-01

GCCTCGCAGAAGGCACAGATCGCGTAGTGAAAGCACTTGCATTGCCTCATAGCTCTTTTACCATTAGTCATTAATACGCCGTTAAGCAACTCAGGACTCTATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCTTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGA‑TGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGA‑TTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCG  >  NC_000913/3041210‑3041687                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                               gCCTCGCAGAAGGCACAGATCGCGTAGTGAAAGCACTTGCATTGCCTCATAGCTCTTTTACCATTAGTCATTAATACGCCGTTAAGCAACTCAGGACTCTATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCTTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAACAGATGTc                                                                                                                                                                                                                                        >  1:451144/1‑250 (MQ=255)   cTCGCAGAAGGCACAGATCGCGTAGTGAAAGCACTTGCATTGCCTCATAGCTCTTTTACCATTAGTCATTAATACGCCGTTAAGCAACTCAGGACTCTATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCTTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAACAGATGTCAGc                                                                                                                                                                                                                                     >  1:355353/1‑251 (MQ=255)   cTCGCAGAAGGCACAGATAGCGTAGTGAAAGCACTTGCATTGCCTCATAGCTCTTTTACCATTAGTCATTAGTACGCCGTTAAGCAACTCAGGACTCTATTATGGATTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCTTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGc                                                                                                                                                                                                                                     >  1:68267/1‑251 (MQ=255)       cAGAAGTCACAGATCGCGTAGTGAAAGCACTTGCATCGCCTCATAGATCTTTTACCATAAGTCATTAATAAGCCGTTAAGCAACTCAAGACTCTATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCACACGGCTTCTGCCCGTTTACCGCTGGCGATGATACCGCTTGATGAATAAGCGCTTGCCCCCATTTCTGGCTAGGACAGCGAACAATATATGCATGGGCAGGTTACAACAGCTGTCAGcc                                                                                                                                                                                                                                    >  2:246211/1‑248 (MQ=255)        aGAAGGCACAGATCGCGTAGTGAAAGCACTTGCATTGCCTCATAGCTCTTTTACCATTAGTCATTAATACGCCGTTAAGCAACTCAGGACTCTATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCTTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAACAGATGTCAGcc                                                                                                                                                                                                                                    >  1:221803/1‑247 (MQ=255)           aGGCACAGATCGCGTAGTGAAAGCACTTGCATTGCCTCATAGCTCTTTTACCATTAGTCATTAATACGCCGTTAAGCAACTCAGGACTCTATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCTTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAACAGATGTCAGCCAGA‑TGGc                                                                                                                                                                                                                            >  1:130846/1‑251 (MQ=255)                  gATCGCGTAGTGAAAGCACTTGCATTGCCTCATAGCTCTTTTACCATTAGTCATTAATACGCCGTTAAGCAACTCAGGACTCTATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCTTTCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAACAGATGTCAGCCAGA‑TGGCAGAAGAt                                                                                                                                                                                                                     <  2:174541/251‑1 (MQ=255)                                  cACTTGCATTGCCTCATAGCTCTTTTACCATTAGTCATTAATACGCCGTTAAGCAACTCAGGACTCTATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCTTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAACAGATGTCAGCCAGA‑TGGCAGAAGATCAGCACCTGCTcgccg                                                                                                                                                                                                     >  1:135564/1‑251 (MQ=255)                                      tGCATTGCCTCATAGCTCTTTTACCATTAGTCATTAATACGCCGTTAAGCAACTCAGGACTCTATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCTTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAACAGATGTCAGCCAGA‑TGGCAGAAGATCAGCACCTGCTCgccgcc                                                                                                                                                                                                   <  1:111773/249‑1 (MQ=255)                                        cATTGCCTCATAGCTCTTTTACCATTAGTCATTAATACGCCGTTAAGCAACTCAGGACTCTATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCTTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAACAGATGTCAGCCAGA‑TGGCAGAAGATCAGCACCTGCTCGCCGCCCAtt                                                                                                                                                                                               >  1:369831/1‑251 (MQ=255)                                         aTTGCCTCATAGCTCTTTTACCATTAGTCGTTAATACGCCGTTAAGCAACTCAGGACTCTATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCTTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAACAGATGTCAGCCAGA‑TGGCAGAAGATCAGCACCTGCTCGCCGCCCAttt                                                                                                                                                                                              >  2:71299/1‑250 (MQ=255)                                                            aCCATTAGTCATTAATACGCCGTTAAGCAGCTCAGGACTCTATTATGGCTTTTACACTCTTTCCTCCCCGTCAGCCTACTGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCGTGATGACTGGGCTCGTGCCACCAGTACTGGCGAGGACAGCGAAAAATATATGCAGGGTCAGGTGACAACAGATGTCAGCCAGA‑TGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGTCGCCAAAGGTAAAAt                                                                                                                                                                            <  2:279162/250‑1 (MQ=255)                                                              cATTAGTCATTAATACGCCGTTAAGCAACTCAGGACTCTATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGCTGACGCTTGATGACTGGGCGCTTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAACAGATGTCAGCCAGA‑TGGCAGAAGATCAGCACCTGCGCGCCGCCCATTGCGACGCCAAAGGAAAAAtgtg                                                                                                                                                                         >  2:352634/1‑251 (MQ=255)                                                              cATTAGTCATTAATACGCCGTTAAGCAACTCAGGACTCTATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCGTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCTTGCCACCATTACTGTCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAACAGATGTCAGCCAGA‑TGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAAAtgtg                                                                                                                                                                         <  1:401401/251‑1 (MQ=255)                                                              cATTAGTCATTAATACGCCGTTAAGCAACGCAGGACTCTATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTTCCGCTGACGCTGATGACGCTTGATGACTGGGCGCTTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAACAGATGTCAGCCAGA‑TGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAAAtgtg                                                                                                                                                                         <  1:49087/251‑1 (MQ=255)                                                                     cATTAATACGCCGTTAAGCAACTCAGGACTCTATTATGGCTTTTACACCTGTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAACAGATGTCAGCCAGA‑TGGCAGAAGATCAGCACCTGTTCGCCGCCCATTGCGACGCCAAAGGTAAAATGTGGAGCAAt                                                                                                                                                                  <  1:367491/251‑1 (MQ=255)                                                                                   tAAGCAACTCAGGACTCTATGATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCTTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAACAGATGTCAGCCAGA‑TGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGATGCCAAAGGTAAAATGTGGAGCAATTTACGTCTGTTcc                                                                                                                                                     <  1:195754/250‑1 (MQ=255)                                                                                    aaGCAACTCAGGACTCTATTATGGCTTTTACACCTTCTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCTTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAACAGATGTCAGCCAGA‑TGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAAATGTGGAGCAATTTACGTCTGTTCcg                                                                                                                                                    <  1:285241/250‑1 (MQ=255)                                                                                        caCTCAGGACTCTATTATGGCTTTTTCACCTTTTCCTCCCCGTCAGCCTTCGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCTTGCCACCATTTCTGGCGCGGACAGCGAAAAATATATGCGGGGTCGGGTGACAACGTATGGCAGCCAGA‑TGGCAGAAGATCTGCGCCTGCTCGCCGTCCATTGCGTCGCCAATGGTAAAGTGTGGAGCAATTTACGTCTGTTCCGCGACg                                                                                                                                               <  1:202216/250‑1 (MQ=255)                                                                                               gcctctATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCTTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAACAGATGTCAGCCAGA‑TGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATgg                                                                                                                                        <  2:96315/249‑1 (MQ=255)                                                                                                aCTCTATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCTTGCCACCCTTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAACAGATGTCAGCCAGA‑TGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGGAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGc                                                                                                                                       >  2:62908/1‑251 (MQ=255)                                                                                                        atGGATTTTACTCCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCTTGCCACCATTACTGGCGCGGACAGCGAAAAAGATATGCAGGGTCAGGTGACAGCAGATGTCTGCGAGA‑TGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAAATGTGGAGCAAGTTACGTCTGTTCCTCGACGGCGATGGCTTTGCAt                                                                                                                                <  2:68267/250‑1 (MQ=255)                                                                                                            cTTTTACAACTTTTCCTACCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCTTGCCACCATTACTGGCGCGGACCGCGAAAAATATATGCAGGGTCAGGTGACAACAGATGTCAGCCAGA‑TGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAAATGTGGAGCAATTTACGTCTGGTCCGCGACGGCGATGGCTTTGCATGGA‑tt                                                                                                                          >  2:218994/1‑251 (MQ=255)                                                                                                                   acCTTTTACTCCCCGTCAGCCTAAGGCTTCTGACCGTTTACCGCTGACGCTGATGCCCCTTGATGACTGGGCGATTGCCACCATTACTGGCGCAGCCAGCGAAACATATATGCAGGGTCAGGTGACAACAGATGTCAGCCAGA‑TGGCAGAAGATCAGCACCTGCTCGCAGCCCATTGCGACGCCAAAGGCAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGAAGGCCTTTCCTGGA‑TTGAACGgc                                                                                                                   >  2:251601/1‑251 (MQ=255)                                                                                                                       tttCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGTCGCTGATGACGCTTGATGACTGGGCGCGTGCCACCAGTACTGGCGCGGACAGAGAAAAATATATGCAGGGTCAGGTGACAACAGATGTCAGCCAGA‑TGGCAGAAGATCAGCGCCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGTCGAGGGCTTTGCATGGA‑TTGAACGGCGCa                                                                                                                <  2:293740/250‑1 (MQ=255)                                                                                                                        ttCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGATGATGACGCTTGATGACTGGGCGCTTGCCACCATTACTGGCGCGGACAGCGACAAATATATGCAGGGTCAGGTGACAACAGATGTCAGCCAGA‑TGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGA‑TTGAACGGCGCAgc                                                                                                              >  1:434868/1‑251 (MQ=255)                                                                                                                                         aCGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCTTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAACAGATGTCAGCCAGA‑TGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGTCGCCAAAGGGAAAATGTGGAGCAATTTACGTCTGTGCCGCGACGGCGATGGCTTTGCATGGA‑TTGAACGGCGCAGCGTGCGTGAACCGCAg                                                                                               <  1:283382/249‑1 (MQ=255)                                                                                                                                            gCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTTATGACTGGGCGCGTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAACAGATGTCAGCCAGA‑TGGCAGAAGATCAGCACCTGCTCGCCTCCCATTGCGACGCCAAAGGTAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATTGCTTTGCATGGA‑TTGAACGGCGCAGCGTGCGTGAACCGCAGctgac                                                                                          <  1:201825/251‑1 (MQ=255)                                                                                                                                              ttCTGCCCGTTTACAGCTGACGCTGATGACGCTTGATGACTGGGCGCTTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAACAGATGTCAGCCAGA‑TGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCACAGGTAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGA‑TTGATCGGCGCAGCGGGCGTGAACCGCAGctgactg                                                                                        >  2:314794/1‑251 (MQ=255)                                                                                                                                                cggCCCGTTTACCGCTGAGGCTGATGCTGCTTGATTACTTGGCGCATGCCACCATTACTGGCGCGGACAGCGAAAAATATACGCAGGGTCAGGCGACAACAGCGGTCAGCCAGACTGG‑AGAAGATCTTCACCTTGTCGCCGCCCATTGCGTCGCCAAAGGTAAAAGTTGGAGCAATTTACGGCGGTTCCGCTATTGCGATGTTTTTCCAGGGA‑TTTAACGGCTCATCGTGCGTGCTCCGCAGCTGactgaa                                                                                      <  2:385917/249‑1 (MQ=255)                                                                                                                                                  gCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCTTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAACAGATGTCAGCCAGA‑TGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGA‑TTGAACGGCGCAGCGTGCGTGAACCGCAGCTGactgaa                                                                                      <  1:153391/249‑1 (MQ=255)                                                                                                                                                   cccGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCTTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAACAGATGTCAGCCAGA‑TGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGA‑TTGAACGGCGCAGCGTGCGTGAACCGCAGCTGactgaactg                                                                                   <  1:352634/251‑1 (MQ=255)                                                                                                                                                   cccGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCTTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAACAGATGTCAGCCAGA‑TGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGA‑TTGAACGGCGCAGCGTGCGTGAACCGCAGCTGactgaactg                                                                                   <  1:282158/251‑1 (MQ=255)                                                                                                                                                                     tgatgaCGCTTGATGACTGGGCGCTTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAACAGATGTCAGCCAGA‑TGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGA‑TTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATtc                                                                 >  1:140945/1‑251 (MQ=255)                                                                                                                                                                        tgaCGATTGATGAATGGGAGCTTGCAACCATTACTGGCGCGGACAGCGAAAAAAATATGCGGGGTCAGGGGACAACACATGTCAGCCAGA‑TGGCAGAAGATCAGCACATGCTCGCCGCCCATTGCGACGCCAACGGTAAAATGGGGAGCAATTTAAGTCTGTTCCGAGACGGCGATGGATTTGCCTGGA‑TTGAACGGCGCAGGGTGCGTGAACCGCAGCTGACTGAACTGAAAAAAAATGAGGTATTctct                                                              >  2:293119/1‑251 (MQ=255)                                                                                                                                                                         gaCGCTTGATGACTGGGCGCTTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAACAGATGTCAGCCAGA‑TGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGA‑TTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTctct                                                              >  1:26271/1‑250 (MQ=255)                                                                                                                                                                         gaCGCTTGATGACTGGGCGCTTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAACAGATGTCAGCCAGA‑TGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGA‑TTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTa                                                             <  2:198751/251‑1 (MQ=255)                                                                                                                                                                            gCTTGATGACTGGGCGCTTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAACAGATGTCAGCCAGA‑TGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGA‑TTGAACGGCGCAGCGGGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTaa                                                            >  2:89916/1‑249 (MQ=255)                                                                                                                                                                                gatgaCTGGGCGCTTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAACAGATGTCAGCCAGA‑TGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGA‑TTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGAc                                                      <  1:363237/251‑1 (MQ=255)                                                                                                                                                                                   gaCTGGGCGCTTGCCACCATTACTGGCGCGGACAGCGAACAATATATGCAGGGTCAGGTGACAACAGATGTCAGCCAGA‑TGGCAGAAGATCAGCACATGCTCGCCGCCCATTGCGACGCCAAAGGTAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGA‑TTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGACCAt                                                   >  1:28247/1‑251 (MQ=255)                                                                                                                                                                                           gcTTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAACAGATGTCAGCCAGA‑TGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGAAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGA‑TTGAACGGCGCAGCGTGCGTGAACCGCACCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGACCATCGCGcc                                             >  2:416395/1‑249 (MQ=255)                                                                                                                                                                                                    caTTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAACAGATGTCAGCCAGA‑TGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGA‑TTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGACCATCGCGCCAGACGACGAGc                                  >  1:187968/1‑251 (MQ=255)                                                                                                                                                                                                              gcgGACAGCGAAAAATATATGCAGGGTCAGGTGACAACAGATGTCAGCCAGA‑TGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGA‑TTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGACCATCGCGCCAGACGACGAGCGTGTgctgct                        <  1:50610/251‑1 (MQ=255)                                                                                                                                                                                                               cgGACAGCGAAAAATATATGCAGGGTCAGGTGACCACAGATGTCAGCCAGA‑TGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGA‑TTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGACCATCGCACCAGACGACGAGCGGGtgctgcgt                       >  1:294198/1‑249 (MQ=255)                                                                                                                                                                                                               cgGACAGAGAAAAATCTCTGCAGGGTCCGATCACACCAGATGTCAGCAAGA‑TGGCAGAAGACCCACACCTGCTCGCCGCGCATTGCGACGCCAAAGGTAAAATGTGAAGCAATTTACGTCTGGTCCGCGACGGCGATGGCTTTGCATGAC‑TTGAACGGCGCAGCGTGCGCGAACCGCAGCTGACTGCACTGAAAAAATATGCGGAATTAGCTAACGGGACCATCGCGCCAGACGACGAGCGTGTGCTGCtt                       >  2:210939/1‑251 (MQ=255)                                                                                                                                                                                                                ggACAGTGACAAATATATGCAGGGTCAGGTGACCACAGATGTCAGCCAGA‑TGGCAGAAGATCAGCACCTGCTCGACGCCCACTGCGACGCCAAAGGTAAAATGTGGAGAAATTTCCGTCTGTTCCGCGACGGCGATGGCTTTGCATGGA‑TTGAACGGCGCAGCGTGCGTGAACAGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGACCATCGCTCCAGACGACGAGCGTGTGCTGTTTg                      >  2:301341/1‑251 (MQ=255)                                                                                                                                                                                                                       gAAAAATATATGCAGGGTCAGGTGACAACAGATGTCAGCCAGA‑TGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGA‑TTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGc               >  1:335/1‑251 (MQ=255)                                                                                                                                                                                                                        aaaaaTATATGCAGGGTCAGGTGACAACAGATGTCAGCCAGA‑TGGCAGAAGATCAGCACCTGCTCGCCGCCAATTGCGACGCCAAAGGTAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGA‑TTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGcc              <  1:446066/251‑1 (MQ=255)                                                                                                                                                                                                                            atatatGCAGGGTCAGGTGACAACAGATGTCAGCCAGA‑TGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGA‑TTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGtt          <  2:221803/251‑1 (MQ=255)                                                                                                                                                                                                                             aatatGCAGGGTCAGGTGACAACAGAGGTCAGCCAGA‑TGGCATAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGA‑TTGAACGGCGCAGCTTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGACCATCGCGCCAGACGACGATCGTGTGCTGCTTGGTGTTGCCGGttt         <  2:274017/250‑1 (MQ=255)                                                                                                                                                                                                                              atatGCAGGGTCAGGTGACAACAGATGTCAGCCAGA‑TGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGA‑TTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGttt         <  1:123998/250‑1 (MQ=255)                                                                                                                                                                                                                                    cgagttaGGTGACAACAGATGTCAGCCAGC‑TGGCAGAATGTGAGCGCCTGCTCGCCTTCCGTTTCCACGCCAAAGGTTAAATGTGTAGCAAGTTACGTCGGTTCCGCGACGGCGATTGCGTTTCTTGGAGTTGA‑CTGTGCTGCGTGCGTGAAGCGCAGCTTTCTGGATGGAAAAAATATGCTGTTGTCGCTATAGGGTCTGTCGCGGCGTAGGACGTGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGgcg  <  1:407606/245‑1 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                               GCCTCGCAGAAGGCACAGATCGCGTAGTGAAAGCACTTGCATTGCCTCATAGCTCTTTTACCATTAGTCATTAATACGCCGTTAAGCAACTCAGGACTCTATTATGGCTTTTACACCTTTTCCTCCCCGTCAGCCTACGGCTTCTGCCCGTTTACCGCTGACGCTGATGACGCTTGATGACTGGGCGCTTGCCACCATTACTGGCGCGGACAGCGAAAAATATATGCAGGGTCAGGTGACAGCAGATGTCAGCCAGA‑TGGCAGAAGATCAGCACCTGCTCGCCGCCCATTGCGACGCCAAAGGTAAAATGTGGAGCAATTTACGTCTGTTCCGCGACGGCGATGGCTTTGCATGGA‑TTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCG  >  NC_000913/3041210‑3041687

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: ‑