Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,472,426	(A)8→9	coding (108/135 nt)	pyrL ←	pyrBI operon leader peptide

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,472,418	1	.	A	100.0%	127.3 / NA	36	C39F (TGC→TTC)	pyrL	pyrBI operon leader peptide
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base A (19/17);  total (19/17)

GTTGCGGGTTTGCTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGC  >  NC_000913/4472287‑4472546                                                                                                                                     |                                                                                                                                  gTTGCGGGTTTGCTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑‑aaaaaaaa                                                                                                                          >  2:440001/1‑140 (MQ=255)                                         cAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCA‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATc                                                                                         >  2:565657/1‑134 (MQ=255)                                         cAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCA‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATc                                                                                         <  1:565657/134‑1 (MQ=255)                                                     aTCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCA‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGc                                                                       >  2:297306/1‑140 (MQ=255)                                                     aTCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCA‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGc                                                                       >  1:715421/1‑140 (MQ=255)                                                      tCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCA‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGTGAAGAAAAACGGCa                                                                      >  1:21209/1‑140 (MQ=255)                                                       cgcgACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCA‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAg                                                                     >  1:576957/1‑140 (MQ=255)                                                            cTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCA‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAg                                                                >  2:110176/1‑140 (MQ=255)                                                               aGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCA‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGt                                                             <  1:530553/140‑1 (MQ=255)                                                                              aTGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCA‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGAcgcg                                              >  1:11501/1‑140 (MQ=255)                                                                               tGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCA‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCgg                                             <  1:23252/140‑1 (MQ=255)                                                                                gATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCA‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGt                                            <  1:440001/140‑1 (MQ=255)                                                                                        tCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCA‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACaaa                                    <  2:715421/140‑1 (MQ=255)                                                                                        tCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCA‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACaaa                                    >  1:723544/1‑140 (MQ=255)                                                                                         cTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCA‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACaaaa                                   >  1:106268/1‑140 (MQ=255)                                                                                           gATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCA‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATg                                 >  2:510614/1‑140 (MQ=255)                                                                                                aGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCA‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGaa                            <  2:688739/140‑1 (MQ=255)                                                                                                  cGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCA‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAaca                          <  1:507235/140‑1 (MQ=255)                                                                                                  cGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCA‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAaca                          >  2:190717/1‑140 (MQ=255)                                                                                                      ttAGCCATCTTTTATCTCCTGACGCCTGGGCA‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCGGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACt                      >  1:215914/1‑140 (MQ=255)                                                                                                      ttAGCCATCTTTTATCTCCTGACGCCTGGGCA‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAaca                          >  1:241598/1‑136 (MQ=255)                                                                                                      ttAGCCATCTTTTATCTCCTGACGCCTGGGCA‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAaca                          <  2:241598/136‑1 (MQ=255)                                                                                                      ttAGCCATCTTTTATCTCCTGACGCCTGGGCA‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACt                      <  2:11501/140‑1 (MQ=255)                                                                                                      ttAGCCATCTTTTATCTCCTGACGCCTGGGCA‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACt                      <  1:297306/140‑1 (MQ=255)                                                                                                      ttAGCCATCTTTTATCTCCTGACGCCTGGGCA‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACt                      <  2:710434/140‑1 (MQ=255)                                                                                                       tAGCCATCTTTTATCTCCTGACGCCTGGGCA‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACt                      <  2:200642/139‑1 (MQ=255)                                                                                                       tAGCCATCTTTTATCTCCTGACGCCTGGGCA‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACt                      >  1:200642/1‑139 (MQ=255)                                                                                                       tAGCCATCTTTTATCTCCTGACGCCTGGGCA‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACt                      <  2:229036/139‑1 (MQ=255)                                                                                                       tAGCCATCTTTTATCTCCTGACGCCTGGGCA‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACt                      >  1:229036/1‑139 (MQ=255)                                                                                                        aGCCATCTTTTATCTCCTGACGCCTGGGCA‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGa                    <  1:510614/140‑1 (MQ=255)                                                                                                              cTTTTATCTCCTGACGCCTGGGCA‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATa              >  1:432131/1‑140 (MQ=255)                                                                                                                   aTCTCCTGACGCCTGGGCA‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGA‑‑ACACTGAACCATACAtcctc       <  2:723544/140‑1 (MQ=255)                                                                                                                        cTGACGCCTGGGCA‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGTGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCg    <  2:21209/140‑1 (MQ=255)                                                                                                                        cTGACGCCTGGGCA‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCg    >  2:61406/1‑140 (MQ=255)                                                                                                                        cTGACGCCTGGGCA‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCg    >  2:14291/1‑140 (MQ=255)                                                                                                                         tGACGCCTGGGCAAAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAAACAAAATGTCGAACACACTGAACCATACATCCTCCCg    >  2:490813/1‑140 (MQ=255)                                                                                                                          gACGCCTGGGCA‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGc  <  1:110176/140‑1 (MQ=255)                                                                                                                                     |                                                                                                                                  GTTGCGGGTTTGCTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGC  >  NC_000913/4472287‑4472546

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