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%	107.5 / NA	39	L733V (TTA→GTA)	glnD	uridylyltransferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (26/13);  total (26/13)

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