Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	186,454	A→C	L733V (TTA→GTA)	glnD ←	uridylyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	186,454	0	A	C	100.0%	105.3 / NA	38	L733V (TTA→GTA)	glnD	uridylyltransferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (26/12);  total (26/12)

GTCAGTACTTGCTCCAGACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGCGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTCCAGATAAAAATCTCGGTGCCTCCACGCGTAGCCTGCGGGCTAAGCAATACCAGCGGTTTGCTTAAATCATGCTGTAATAAA  >  NC_000913/186315‑186593                                                                                                                                            |                                                                                                                                            gTCAGTACTTGCTCCAGACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATGCCGCGGCGAGTGGTGAAAGTTTGTGCGTCGTGAGCACTTACAtt                                                                                                                                          >  1:83585/1‑143 (MQ=255)     gTACTTGCTCCAGACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGCGTCGTGAACACTTACATTgcgg                                                                                                                                      >  1:103696/1‑143 (MQ=255)             tCCAGACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGCGTCGTGAAGACTTACATTGCGGCGGGCtga                                                                                                                              >  1:148486/1‑141 (MQ=255)              ccAGACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGCGGTGAAAAATTGTGCGTCGTGAACACTTACATTGCGGCGGTCCAAt                                                                                                                             >  2:665593/1‑143 (MQ=255)               cAGACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGCGTCGTGAACACTTGCATTGCGGCGGTCTAAtt                                                                                                                            >  1:513931/1‑143 (MQ=255)                  aCCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGCGTCGTGAACACTTCCCTTGCGGCGGTCTAATTCgg                                                                                                                         >  1:757036/1‑143 (MQ=255)                           aaTAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCCTACCCGCGCGAGTGGTGAAAATTTGTGCGTCGTGAACACTTACATCGCCGCGGTTTAATTCGGCCCAggagg                                                                                                                >  2:184392/1‑139 (MQ=255)                                         gttCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGGGAGTGGTGAAAATTGGTGCGTCGTGAACACTTACATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATaa                                                                                                  <  2:190026/141‑1 (MQ=255)                                         gagcggCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGCGTCGTGAACACTTACATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATaa                                                                                                  <  2:896997/138‑1 (MQ=255)                                          aTCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGCGTCGTGAACACTTACATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAg                                                                                                 >  2:387496/1‑143 (MQ=255)                                           tCTGCGGACAGCGGTTTGCCGTCGGGTGCCGGCGCGATAAAGGTCTCCATCGCCATACCGTCGCGAGGGGTGAAAATTTGTGCGTCGTGAACACGTACATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAgg                                                                                                <  2:708468/143‑1 (MQ=255)                                             tGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCAATAGCCATGCCGTCGCGAGTGGTGACAGTTTGGGCGTCGTGAACACTTACCTAGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACGGATAGGGGc                                                                                              >  2:744466/1‑143 (MQ=255)                                                      cGGGTGGCGATCGGGTGCCAGGTCTATAAGGGGATGCGTCGCCATATCGGCGCGAGTGGTGAGAATTTGGGCGTCGGGAACACTTCCATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGc                                                                                     <  2:107055/143‑1 (MQ=255)                                                         gTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGCGTCGTGAACACTTACATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTcc                                                                                  <  1:248905/143‑1 (MQ=255)                                                                       ccAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGCGTCGTGAACACTTACATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTCCAGATAAAAATCTCg                                                                    <  1:387496/143‑1 (MQ=255)                                                                        cAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGCGTCGTGAACACTTACATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTGCGGGCTCCAGAGAAAAATCTCgg                                                                   >  2:87205/1‑143 (MQ=255)                                                                         aGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGCGTCGTGAACACTTACATTGCGGCGGTCTAATTCGGCACCGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTCCAGATAAAAATCCCGGt                                                                  >  1:89562/1‑143 (MQ=255)                                                                              gATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGCGTCGTGAACACTTACATTGCGGCGGTCTAAGACGGCACGGACGGCGGCGAACGGATAAGGGCGGTCCGGGCTCCAGCTCAAAATCTCGCTGCCTc                                                             >  1:504061/1‑143 (MQ=255)                                                                                   atctccccATCGCCAAACCTGCGGTAGTGGTTCAAATTTGTGCCTCGCGAACATGTACATTGTGGCGGCCTAATGCGGCTCAGACGCCGGCAAATAGATAAGGGCAGTCCGGGCTCCAGATAAAAATCTCGGTGCCTCCAcgc                                                        <  2:330191/137‑1 (MQ=255)                                                                                     gTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGCGTCGTGAACACTTACAGTGCGGCGGGCTAATTCGGCACAGACGGCGGCAAACAGATAAGAGCGGTCCGGGCGCCAGATACAAATGTCGGTGCCTCCACGCGt                                                      >  1:302009/1‑143 (MQ=255)                                                                                      tATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGCGTCGTGAACACTTACATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAAGGCGGTCCGGGCTCCAGATAAAAATCTCGGTGCCTCCACGCGTa                                                     >  1:467489/1‑143 (MQ=255)                                                                                        tCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGCGTCGTGAACACTTACATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTCCAGATAAAAATCTCGGTGCCTCCACGCGTAGc                                                   >  1:578276/1‑143 (MQ=255)                                                                                          cATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGCGTCGTGAACACTTACATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTCCAGATAAAAATCTCGGTGCCTCCACGCGTAGCCt                                                 >  1:560149/1‑143 (MQ=255)                                                                                          cATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGCGTCGTGAACACTTACATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTCCAGATAAAAATCTCGGTGCCTCCACGCGGAGcca                                                 >  2:812599/1‑142 (MQ=255)                                                                                           aTCGCCATACCGTCGCGAGTGGTGAAAATTTGTGCGTCGTGAACACTTACATTGCGGCGGTCTAATTCGGCACAGGCGGCGGCAAACGGATAAGGGCGGTCCGGGCTCCAGATAAAAATCTCGGTGCCGCCACGCGCAGCCTg                                                >  2:377072/1‑143 (MQ=255)                                                                                              gtcgtcccGTGGCGACTGGTCGGAATTGGTGCGTCGCGACCACCTACAGTGCGGCGTGCTAATTCGGCACAGGCGGCGGCAAACAGAGGGGGGCGCTCCGGGCTCCAGATAAAAATCTCGGTGCCTCCACGCGTAGCCTGCgg                                             <  2:148486/137‑1 (MQ=255)                                                                                              gCCATACCGTCGCGAGTGGTGAAAATTTGTGCGTCGTGAACACTTACATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTCCAGATAAAAATCTCGGTGCCTCCACGCGTAGCCTGCgg                                             <  1:959255/143‑1 (MQ=255)                                                                                                   aCCGTCGCGAGTGGTGAAAATTTGTGCGTCGTGAACACTTACATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTCCAGATAAAAATCTCGGTGCCTCCACGCGTAGCCTGCGGGCTaa                                        <  1:801338/143‑1 (MQ=255)                                                                                                    ccGTCCCGGGCGGTGAGAACTTGTGCGTGGTGCCCACGCCCATCGCGGCGGTGTAGTGGGGAACAGACGGCGGCAGACAGATACGGGCGGTCCGGGCTCCAGAGAAAAATCGCGGTGCCTCCACGCGTAGCCTGCGGGCTAAg                                       <  2:851091/143‑1 (MQ=255)                                                                                                      gTCGCGAGTGGTGAAAATTTGTGCGTCGTGAACACTTACATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTCCAGATAAAAATCTCGGTGc                                                                <  2:282362/116‑1 (MQ=255)                                                                                                      gTCGCGAGTGGTGAAAATTTGTGCGTCGTGAACACTTACATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTCCAGATAAAAATCTCGGTGc                                                                >  1:282362/1‑116 (MQ=255)                                                                                                           gAGTGGTGAAAATTTGTGCGTCGTGAACACTTACATTGCGGCGGTCTAATTCGGCCCAGACGGTGGCCAACAGATAAGGGCGGTCCGGGCTCCAGATAAAAATCTCGGTGCCTCCAAGCGTAGCCTGCGGGCTAAGCAATAcc                                >  2:917007/1‑143 (MQ=255)                                                                                                                gtgAAAATTTGTGCGTCGTGAACACTTACATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTCCAGAGAAAAATCTCGGTGCCTCCACGCGTAGCCTGCGGGCTAAGCAAACCCACCgg                           >  2:695180/1‑143 (MQ=255)                                                                                                                 tgAAAATTTGTGCGTCGTGAACACTTACATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAACCAGATCAGGGCGGTCCGGGCTCCAGATAAAAATCTCGGTGCCTCCACGCGTAGCCTGCGGGCTAAGCAATACCAGCGGt                          >  2:805938/1‑143 (MQ=255)                                                                                                                 tgAAAATTTGTGCGTCGTGAACACTTACATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTCCAGATAAAAATCTCGGTGCCTCCACGCGTAGCCTGCGGGCTAAGCAATACCAGCGGt                          >  2:819180/1‑143 (MQ=255)                                                                                                                                      cacTTACATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGCCTCCAGATAAAAATCTAGGTGCCTCCACGCGTAGCCTGCGGGCTAAGCAATACCAGCGGTTTGCTTAAATCATGCAGtact     >  2:248520/1‑141 (MQ=255)                                                                                                                                         ttACATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTCCAGATAAAAATCTCGGTGCCTCCACGCGTAGCCTGCGGGCTAAGCAATACCAGCGGTTTGCTTAAATCATGCTGTGATaaa  >  2:897203/1‑143 (MQ=255)                                                                                                                                         ttACATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTCCAGATAAAAATCTCGGTGCCTCCACGCGTAGCCTGCGGGCTAAGCAATACCAGCGGCTTGCTTAAATCATGCAGCAATaaa  >  2:942089/1‑143 (MQ=255)                                                                                                                                            |                                                                                                                                            GTCAGTACTTGCTCCAGACCAAACCGAATAACCTCATGACGATCTGCGGACAGCGGGTTGCCATCGGGTTCCAGCACGATAAAGGTATCCATCGCCATACCGTCGCGAGTGGTGAAAATTTGTGCGTCGTGAACACTTAAATTGCGGCGGTCTAATTCGGCACAGACGGCGGCAAACAGATAAGGGCGGTCCGGGCTCCAGATAAAAATCTCGGTGCCTCCACGCGTAGCCTGCGGGCTAAGCAATACCAGCGGTTTGCTTAAATCATGCTGTAATAAA  >  NC_000913/186315‑186593

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