Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,471,278	T→A	D57V (GAT→GTT)	pyrI ←	aspartate carbamoyltransferase, PyrI subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,471,278	0	T	A	100.0%	132.6 / NA	36	D57V (GAT→GTT)	pyrI	aspartate carbamoyltransferase, PyrI subunit
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (16/20);  total (16/20)

CCGTGGCTTGCGGCGCATACAATGCCAGTTGATCTACTTGATCTTCACTCAAAAAGGTATTTTCGATTTTGATCAGATCTTTGCGGCCCATCTCGCCAGAAGGCAGGTTCAGACCAATGGTGATGCGCTGATCCGTTTCGGTCAGCTTGAA  >  NC_000913/4471201‑4471351                                                                              |                                                                          ccGTGGCTTGCGGCGCATACAATGCCAGTTGATCTACTTGATCTTCACTCAAAAAGGTATTTTCGATTTTGATCAGAACtt                                                                        <  2:848126/81‑1 (MQ=255)  cGTGGCTTGCGGCGCATACAATGCCAGTTGATCTACTTGATCTTCACTCAAAAAGGTATTTTCGATTTTGATCAGAACttt                                                                       <  1:848128/81‑1 (MQ=255)  cGTGGCTTGCGGCGCATACAATGCCAGTTGATCTACTTGATCTTCACTCAAAAAGGTATTTTCGATTTTGATCAGAACttt                                                                       <  1:848127/81‑1 (MQ=255)    tGGCTTGCGGCGCATACAATGCCAGTTGATCTACTTGATCTTCACTCAAAAAGGTATTTTCGATTTTGATCAGAACTTTGc                                                                     >  1:848170/1‑81 (MQ=255)      gCTTGCGGCGCATACAATGCCAGTTGATCTACTTGATCTTCACTCAAAAAGGTATTTTCGATTTTGATCAGAACTTTGCgg                                                                   >  2:848164/1‑81 (MQ=255)      gCTTGCGGCGCATACAATGCCAGTTGATCTACTTGATCTTCACTCAAAAAGGTATTTTCGATTTTGATCAGAACTTTGCgg                                                                   <  2:848129/81‑1 (MQ=255)          gcggcgCATACAATGCCAGTTGATCTACTTGATCTTCACTCAAAAAGGTATTTTCGATTTTGATCAGAACTTTGCGGCCCa                                                               >  1:848162/1‑81 (MQ=255)           cggcgCATACAATGCCAGTTGATCTACTTGATCTTCACTCAAAAAGGTATTTTCGATTTTGATCAGAACTTTGCGGCCCa                                                               <  2:848130/80‑1 (MQ=255)           cggcgCATACAATGCCAGTTGATCTACTTGATCTTCACTCAAAAAGGTATTTTCGATTTTGATCAGAACTTTGCGGCCCAt                                                              <  1:848131/81‑1 (MQ=255)             gcgcATACAATGCCAGTTGATCTACTTGATCTTCACTCAAAAAGGTATTTTCGATTTTGATCAGAACTTTGCGGCCCAtct                                                            >  1:848171/1‑81 (MQ=255)                    cAATGCCAGTTGATCTACTTGATCTTCACTCAAAAAGGTATTTTCGATTTTGATCAGAACTTTGCGGCCCATCTCGCCAg                                                      >  1:848175/1‑80 (MQ=255)                      aTGCCAGTTGATCTACTTGATCTTCACTTAAAAAGGCATGTTAGATTCTGATCAGAACTTTGCTGCCCATCTCGCCAGGAg                                                   >  1:891215/1‑81 (MQ=255)                         ccAGTTGATCTACTTGATCTTCACTCAAAAAGGTATTTTCGATTTTGATCAGAACTTTGCGGCCCATCTCGCCAGAAGGc                                                 <  1:848132/80‑1 (MQ=255)                               gATCTACTTGATCTTCACTCAAAAAGGTATTTTCGATTTTGATCAGAACTTTGCGGCCCATCTCGCCAGAAGGCAGGtt                                            >  2:848185/1‑79 (MQ=255)                                     cTTGATCTTCACTCAAAAAGGTATTTTCGATTTTGATCAGAACTTTGCGGCCCATCTCGCCAGAAGGCAGGTTCAg                                         <  2:848133/76‑1 (MQ=255)                                     cTTGATCTTCACTCAAAAAGGTATTTTCGATTTTGATCAGAACTTTGCGGCCCATCTCGCCAGAAGGCAGGTTCAGAcc                                      <  2:848134/79‑1 (MQ=255)                                      ttGATCTTCACTCAAAAAGGTATTTTCGATTTTGATCAGAACTTTGCGGCCCATCTCGCCAGAAGGCAGGTTCAGACCAAt                                   >  2:848169/1‑81 (MQ=255)                                        gATCTTCACTCAAAAAGGTATTTTCGATTTTGATCAGAACTTTGCGGCCCATCTCGCCAGAAGGCAGGTTCAGACCAATg                                  <  2:848135/80‑1 (MQ=255)                                              cACTCAAAAAGGTATTTTCGATTTTGATCAGAACTTTGCGGCCCATCTCGCCAGAAGGCAGGTTCAGACCAATGGTGa                              <  1:848136/78‑1 (MQ=255)                                              cACTCAAAAAGGTATTTTCGATTTTGATCAGAACTTTGCGGCCCATCTCGCCAGAAGGCAGGTTCAGACCAATGGTGATgc                           >  1:848180/1‑81 (MQ=255)                                               aCTCAAAAAGGTATTTTCGATTTTGATCAGAACTTTGCGGCCCATCTCGCCAGAAGGCAGGTTCAGACCAATGGTGATgcg                          <  2:848137/81‑1 (MQ=255)                                                  cAAAAAGGTATTTTCGATTTTGATCAGAACTTTGCGGCCCATCTCGCCAGAAGGCAGGTTCAGACCAATGGTGATGCGCt                        >  1:848174/1‑80 (MQ=255)                                                   aaaaaGGTATTTTCGATTTTGATCAGAACTTTGCGGCCCATCTCGCCAGAAGGCAGGTTCAGACCAATGGTGATGCGCt                        <  2:848138/79‑1 (MQ=255)                                                     aaaGGTATTTTCGATTTTGATCAGAACTTTGCGGCCCATCTCGCCAGAAGGCAGGTTCAGACCAATGGTGATGCGCTGATc                    >  2:848166/1‑81 (MQ=255)                                                       aGGTATTTTCGATTTTGATCAGAACTTTGCGGCCCATCTCGCCAGAAGGCAGGTTCAGACCAATGGTGATGCGCTGATCCg                  <  2:848139/81‑1 (MQ=255)                                                         gTATTTTCGATTTTGATCAGAACTTTGCGGCCCATCTCGCCAGAAGGCAGGTTCAGACCAATGGTGATGCGCTGATCCGtt                <  1:848140/81‑1 (MQ=255)                                                               tCGATTTTGATCAGAACTTTGCGGCCCATCTCGCCAGAAGGCAGGTTCAGACCAATGGTGATGCGCTGATCCGTTTCGGt           <  1:848141/80‑1 (MQ=255)                                                                cGATTTTGATCAGAACTTTGCGGCCCATCTCGCCAGAAGGCAGGTTCAGACCAATGGTGATGCGCTGATCCGTTTCGGt           >  1:848172/1‑79 (MQ=255)                                                                cGATTTTGATCAGAACTTTGCGGCCCATCTCGCCAGAAGGCAGGTTCAGACCAATGGTGATGCGCTGATCCGTTTCGGTc          <  2:848143/80‑1 (MQ=255)                                                                cGATTTTGATCAGAACTTTGCGGCCCATCTCGCCAGAAGGCAGGTTCAGACCAATGGTGATGCGCTGATCCGTTTCGGTc          <  1:848142/80‑1 (MQ=255)                                                                 gATTTTGATCAGAACTTTGCGGCCCATCTCGCCAGAAGGCAGGTTCAGACCAATGGTGATGCGCTGATCCGTTTCGGTCa         <  1:848144/80‑1 (MQ=255)                                                                  aTTTTGATCAGAACTTTGCGGCCCATCTCGCCAGAAGGCAGGTTCAGACCAATGGTGATGCGCTGATCCGTTTCGGTCAGc       >  1:848191/1‑81 (MQ=255)                                                                      tGATCAGAACTTTGCGGCCCATCTCGCCAGAAGGCAGGTTCAGACCAATGGTGATGCGCTGATCCGTTTCGGTCAGCTTg    <  1:848145/80‑1 (MQ=255)                                                                          cAGAACTTTGCGGCCCATCTCGCCAGAAGGCAGGTTCAGACCAATGGTGATGCGCTGATCCGTTTCGGTCAGCTTGaa  >  1:848193/1‑78 (MQ=255)                                                                          cAGAACTTTGCGGCCCATCTCGCCAGAAGGCAGGTTCAGACCAATGGTGATGCGCTGATCCGTTTCGGTCAGCTTGaa  >  1:848190/1‑78 (MQ=255)                                                                          cAGAACTTTGCGGCCCATCTCGCCAGAAGGCAGGTTCAGACCAATGGTGATGCGCTGATCCGTTTCGGTCAGCTTGaa  >  2:848194/1‑78 (MQ=255)                                                                              |                                                                          CCGTGGCTTGCGGCGCATACAATGCCAGTTGATCTACTTGATCTTCACTCAAAAAGGTATTTTCGATTTTGATCAGATCTTTGCGGCCCATCTCGCCAGAAGGCAGGTTCAGACCAATGGTGATGCGCTGATCCGTTTCGGTCAGCTTGAA  >  NC_000913/4471201‑4471351

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