Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	987,152	G→T	intergenic (‑170/+433)	ompF ← / ← asnS	outer membrane porin F/asparagine‑‑tRNA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	987,152	0	G	T	100.0%	53.7 / NA	17	intergenic (‑170/+433)	ompF/asnS	outer membrane porin F/asparagine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (9/8);  total (9/8)

TTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATATTTACAGCAATTTTAATCTATTTATATGATTTCCTTATATTTAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTA  >  NC_000913/987001‑987392                                                                                                                                                        |                                                                                                                                                                                                                                                 ttttttttATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAATATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAAAGCCGTGATGATATTTACAGc                                                                                                                                               >  1:132273/1‑251 (MQ=255)        tATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAATATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATATTTACAGCAAtttt                                                                                                                                         <  1:225159/250‑1 (MQ=255)           gACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAACTATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATATTTACAGCACTTTTAATc                                                                                                                                     >  1:142945/1‑251 (MQ=255)           gACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAATATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATATTTACAGCAATTTTAATc                                                                                                                                     <  2:59226/251‑1 (MQ=255)                            tCAATAAGTCCTGTCAATAAAAATTTCCGGAAATATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAATATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATATTTACAGCAATTTTAATCTATTTATATGATTTCCt                                                                                                                    >  2:117095/1‑251 (MQ=255)                                   gTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAATATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGatat                                                                                                                                                      <  1:102350/210‑1 (MQ=255)                                   gTTCTGTCAATAAAAATTTACGGAAATATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAATATTTGATTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGatat                                                                                                                                                      >  2:102350/1‑210 (MQ=255)                                              aaaaaTTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAATATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATATTTACAGCAATTTTAATCTATTTATATGATTTCCTTATATTTAAATTAACTaa                                                                                                  <  1:7351/251‑1 (MQ=255)                                                aaaTTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAATATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATATTTACAGCAATTTTAATCTATTTATATGATTTCCTTATATTTAAATTAACTAAAc                                                                                                >  1:25758/1‑251 (MQ=255)                                                         gAACTATTGATGTGAGTTTAGTTGCTTTATGTGTGTGCAGACATCTTTCCATTCAAACTAACTTATCACTACTGTGAAAGTGCTACAAAAATAAATATTTGGTTTCAAAAAGAAAAAATATTTAAGCAAAAGTAAAATTTAAGTATCTTTGTTTACAACTTCATATTTCCATAGAGACGTGATGATATTTACATCAATTTTAATCTATTTATATTATTTCCTTATATCTAAACGAACTAATCGAAAAtttt                                                                                       <  2:41569/251‑1 (MQ=255)                                                                                                                 aCTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAATATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATATTTACAGCAATTTTAATCTATTTATATGATTTCCTTATATTTAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGa                                >  2:250595/1‑250 (MQ=255)                                                                                                                       gTTTCGCTACCGTGAAAGTGCTACAAAGATAAATATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATATTTACAGCAATTTTAATCTATTTATATGATTTCCTTATATTTAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCaa                         <  1:192940/251‑1 (MQ=255)                                                                                                                           cGCTACCGTGAAAGTGCTACAAAGATAAATATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATATTTACAGCAATTTTAATCTATTTATATGATTTCCTTATATTTAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCaaaaa                      >  1:276269/1‑250 (MQ=255)                                                                                                                                     aaaGTGCTACAAAGATAAATATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATATTTACAGCAATTTTAATCTATTTATATGATTTCCTTATATTTAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCAAAAAAATAGTATGaa            >  1:94086/1‑250 (MQ=255)                                                                                                                                           cTACAAAGATAAATATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATATTTACAGCAATTTTAATCTATTTATATGATTTCCTTATATTTAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTAt               <  1:96135/241‑1 (MQ=255)                                                                                                                                           cTACAAAGATAAATATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATATTTACAGCAATTTTAATCTATTTATATGATTTCCTTATATTTAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTAt               >  2:96135/1‑241 (MQ=255)                                                                                                                                              cAAAGATAAATATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATATTTACAGCAATTTTAATCTATTTATATGATTTCCTTATATTTAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTa  <  1:240411/251‑1 (MQ=255)                                                                                                                                                        |                                                                                                                                                                                                                                                 TTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATATTTACAGCAATTTTAATCTATTTATATGATTTCCTTATATTTAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTA  >  NC_000913/987001‑987392

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑