Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,770,383	C→T	G7G (GGG→GGA)	unknown ←	Inner membrane protein yciS

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,770,383	0	C	T	100.0%	93.4 / NA	34	G7G (GGG→GGA)	unknown	Inner membrane protein yciS
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (19/15);  total (19/15)

TTCGTTGATTTGTTGCGTTTGGCGTTTGATTTGGCGATTCAGCGCAATATTTTTTAACTTTAATTTGAGATAGAAGAATCCGGTGATAAACCACCCTAGAATTAATCCCAGGCCGAATAAAATAGCAACTAAGGTAGAAAGCTGGAATTGGCTTTGCGCCACAATATAATTAAAGGTAATGATTTGATCATTGTTGGCACCAATGGTAACGGC‑AACCAGCACAATGGCTAAGGTAATGATCAACCCCAAAATATATTTAATCATTTGTTTCTCCTTTAGATAAATGTGTGAGCATTATGCCAAAAATACATATAAAAAACGACACTTAGAAGTGCCGTTATGAAAAGTTAGTTAGTAAAAAATTATGCTGCATAAACGTCGACACGTTCTCTTAATTCTTTGCCTGCCTTAAAATATGGCACACATTTGGCATCCAATTTTACTTGTTCGCCGGTTTTAGGATTACGTCCGATGCGCGGTTGGCGAC  >  1/1770140‑1770626                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    ttCGTTGATTTGTTGCGTTTGGCGTTTGATTTGGCGATTCAGCGCAATATTTTTTAACTTTAATTTGAGATAGAAGAATCCGGTGATAAACCACCCTAGAATTAATCCCAGGCCGAATAAAATAGCAACTAAGGTAGAAAGCTGGAATTGGTTTTGCGCCACAATATAATTAAAGGTAATGATTTGATCATTGTTGGCACCAATGGTAACGGC‑AACCAGCACAATGGCTAAGGTAATGATCAATCCCaaaa                                                                                                                                                                                                                                              <  1:691821/251‑1 (MQ=255)                                     aTTCAGCGCAATATTTTTTAACTTTAATTTGAGATAGAAGAATCCGGTGATAAACCACCCTAGAATTAATCCCAGGCCGAATAAAATAGCAACTAAGGTAGAAAGCTGGAATTGGCTTTGCGCCACAATATAATTAAAGGTAATGATTTGATCATTGTTGGCACCAATGGTAACGGC‑AACCAGCACAATGGCTAAGGTAATGATCAATCCCAAAatat                                                                                                                                                                                                                                           >  1:668733/1‑218 (MQ=255)                                                         aCTTTAATTTGAGATAGAAGAATCCGGTGATAAACCACCCTAGAATTAATCCCAGGCCGAATAAAATAGCAACTAAGGTAGAAAGCTGGAATTGGCTTTGCGCCACAATATAATTAAAGGTAATGATTTGATCATTGTTGGCACCAATGGTAACGGC‑AACCAGCACAATGGCTAAGGTAATGATCAATCCCAAAATATATTTAATCATTTGTTTCTCCTTTAGATAAATGTGTGAGCATTATGCCAAAAAt                                                                                                                                                                                      >  1:85179/1‑251 (MQ=255)                                                                                   gTGATAAACCACCCTAGAATTAATCCCAGGCCGAATAAAATAGCAACTAAGGTAGAAAGCTGGAATTGGCTTTGCGCCACAATATAATTAAAGGTAATGATTTGATCATTGTTGGCACCAATGGTAACGGC‑AACCAGCACAATGGCTAAGGTAATGATCAATCCCAAAATATATTTAATCATTTGTTTCTCCTTTAGATAAATGTGTGAGCATTATGCCAAAAATAc                                                                                                                                                                                    <  1:339463/227‑1 (MQ=255)                                                                                                   gAATTAATCCCAGGCCGAATAAAATAGCAACTAAGGTAGAAAGCTGGAATTGGCTTTGCGCCACAATATAATTAAAGGTAATGATTTGATCATTGTTGGCACCAATGGTAACGGC‑AACCAGCACAATGGCTAAGGTAATGATCAATCCCAAAATATATTTAATCATTTGTTTCTCCTTTAGATAAATGTGTGAGCATTATGCCAAAAATACATATAAAAAACGACACTTAGAAg                                                                                                                                                             >  1:610088/1‑234 (MQ=255)                                                                                                   gAATTAATCCCAGGCCGAATAAAATAGCAACTAAGGTAGAAAGCTGGAATTGGCTTTGCGCCACAATATAATTAAAGGTAATGATTTGATCATTGTTGGCACCAATGGTAACGGC‑AACCAGCACAATGGCTAAGGTAATGATCAATCCCAAAATATATTTAATCATTTGTTTCTCCTTTAGATAAATGTGTGAGCATTATGCCAAAAATACATATAAAAAACGACACTTAGAAGTGCCGTTATGAAAagtt                                                                                                                                            >  1:483994/1‑251 (MQ=255)                                                                                                         aTCCCAGGCCGAATAAAATAGCAACTAAGGTAGAAAGCTGGAATTGGCTTTGCGCCACAATATAATTAAAGGTAATGATTTGATCATTGTTGGCACCAATGGTAACGGC‑AACCAGCACAATGGCTAAGGTAATGATCAATCCCAAAATATATTTAATCATTTGTTTCTCCTTTAGATAAATGTGTGAGCATTATGCCAAAAATACATATAAAAAACGACACTTAGAAGTGCCGTTATGaaa                                                                                                                                                >  1:710379/1‑241 (MQ=255)                                                                                                                 ccGAATAAAATAGCAACTAAGGTAGAAAGCTGGAATTGGCTTTGCGCCACAATATAATTAAAGGTAATGATTTGATCATTGTTGGCACCAATGGTAACGGC‑AACCAGCACAATGGCTAAGGTAATGATCAATCCCAAAATATATTTAATCATTTGTTTCTCCTTTAGATAAATGTGTGAGCATTATGCCAAAAATACATATAAAAAACGACACTTAGAAGTGCCGTTATGAAAAGTTAGTTAGTAAAAAAt                                                                                                                              <  1:622485/251‑1 (MQ=255)                                                                                                                            aGCAACTAAGGTAGAAAGCTGGAATTGGCTTTGCGCCACAATATAATTAAAGGTAATGATTTGATCATTGTTGGCACCAATGGTAACGGC‑AACCAGCACAATGGCTAAGGTAATGATCAATCCCAAAATATATTTAATCATTTGTTTCTCCTTTAGATAAATGTGTGAGCATTATGCCAAAAATACATATAAAAAACGACACTTAGAAGTGCCGTTATGaaa                                                                                                                                                >  1:398442/1‑222 (MQ=255)                                                                                                                                   aaGGTAGAAAGCTGGAATTGGCTTTGCGCCACAATATAATTAAAGGTAATGATTTGATCATTGTTGGCACCAATGGTAACGGC‑AACCAGCACAATGGCTAAGGTAATGATCAATCCCAAAATATATTTAATCATTTGTTTCTCCTTTAGATAAATGTGTGAGCATTATGCCAAAAATACATATAAAAAACGACACTTAGAAg                                                                                                                                                             <  1:733218/202‑1 (MQ=255)                                                                                                                                                ggAATTGGCTTTGCGCCACAATATAATTAAAGGTAATGATTTGATCATTGTTGGCACCAATGGTAACGGC‑AACCAGCACAATGGCTAAGGTAATGATCAATCCCAAAATATATTTAATCATTTGTTTCTCCTTTAGATAAATGTGTGAGCATTATGCCAAAAATACATATAAAAAACGACACTTAGAAg                                                                                                                                                             >  1:102245/1‑189 (MQ=255)                                                                                                                                                ggAATTGGCTTTGCGCCACAATATAATTAAAGGTAATGATTTGATCATTGTTGGCACCAATGGTAACGGC‑AACCAGCACAATGGCTAAGGTAATGATCAATCCCAAAATATATTTAATCATTTGTTTCTCCTTTAGATAAATGTGTGAGCATTATGCCAAAAATACATATAAAAAACGACACTTAGAAGTGCCGTTATGAAAAGTTAGTTAGTAAAAAATTATGCTGCATAAACGTCGACACGTTCTCTTa                                                                                               >  1:87738/1‑251 (MQ=255)                                                                                                                                                               ccACAATATAATTAAAGGTAATGATTTGATCATTGTTGGCACCAATGGTAACGGC‑AACCAGCACAATGGCTAAGGTAATGATCAATCCCAAAATATATTTAATCATTTGTTTCTCCTTTAGATAAATGTGTGAGCATTATGCCAAAAATACATATAAAAAACGACACTTAGAAGTGCCGTTATGAAAAGTTAGTTAGTAAAAAATTATGCTGCATAAACGTCGACACGTTCTCTTAATTCTTtgcctgc                                                                                  >  1:140778/1‑249 (MQ=255)                                                                                                                                                                 acaATATAATTAAAGGTAATGATTTGATCATTGTTGGCACCAATGGTAACGGC‑AACCAGCACAATGGCTAAGGTAATGATCAATCCCAAAATATATTTAATCATTTGTTTCTCCTTTAGATa                                                                                                                                                                                                               <  1:212992/122‑1 (MQ=255)                                                                                                                                                                 acaATATAATTAAAGGTAATGATTTGATCATTGTTGGCACCAATGGTAACGGC‑AACCAGCACAATGGCTAAGGTAATGATCAATCCCAAAATATATTTAATCATTTGTTTCTCCTTTAGATAAATGTGTGAGCATTATGCCAAAAATACATATAAAAAACGACACTTAGAAGTGCCGTTATGAAAAGTTAGTTAGTAAAAAATTATGCTGCATAAACGTCGACACGTTCTCTTAATTCTTTGCCTGCCTTa                                                                              >  1:37613/1‑251 (MQ=255)                                                                                                                                                                  caATATAATTAAAGGTAATGATTTGATCATTGTTGGCACCAATGGTAACGGCAAACCAGCACAATGGCTAAGGTAATGATCAATCCCAAAATATATTTAATCATTTGTTTCTCCTTTAGATAAATGTGTGAGCATTATGCCAAAAATACATATAAAAAACGACACTTAGAAGTGCCGTTATGAAAAGTTAGTTAGTAAAAAATTATGCTGCATAAACGTCGACACGTTCTCTTAATTCTTTGCCTGCCTTa                                                                              >  1:552340/1‑251 (MQ=255)                                                                                                                                                                  caATATAATTAAAGGTAATGATTTGATCATTGTTGGCACCAATGGTAACGGC‑AACCAGCACAATGGCTAAGGTAATGATCAATCCCAAAATATATTTAATc                                                                                                                                                                                                                                   <  1:801207/101‑1 (MQ=255)                                                                                                                                                                    atatAATTAAAGGTAATGATTTGATCATTGTTGGCACCAATGGTAACGGC‑AACCAGCACAATGGCTAAGGTAATGATCAATCCCAAAATATATTTAATCATTTGTTTCTCCTTTAGATAAATGTGTGAGCATTATGCCAAAAAt                                                                                                                                                                                      >  1:847280/1‑144 (MQ=255)                                                                                                                                                                    atatAATTAAAGGTAATGATTTGATCATTGTTGGCACCAATGGCAACGGC‑AACCAGCACAATGGCTAAGGTAATGATCAATCCCAAAATATATTTAATCATTTGTTTCTCCTTTAGATAAATGTGTGAGCATTATGCCAAAAATACATATAAAAAACGACACTTAGAAGTGCCGTTATGAAAAGTTAGTTAGTAAAAAATTATGCTGCATAAACGTCGACACGTTCTCTTAATTCTTTGCCTGCCTTaaaa                                                                           >  1:625189/1‑251 (MQ=255)                                                                                                                                                                     tataATTAAAGGTAATGATTTGATCATTGTTGGCACCAATGGTAACGGC‑AACCAGCACAATGGCTAAGGTATTGATCAATCCCAAAATATATTTAATCATTTGTTTCTCCTTTAGATAAATGTGTGAGCATTATGCCAAAAATACATATAAAAAACGACACTTAGAAGTGCCGTTATGAAAAGTTAGTTAGTAAAAAATTATGCTGCATAAACGTCGACACGTTCTCTTAATTCTTTGCCTGCCTTAAAat                                                                          <  1:796810/251‑1 (MQ=255)                                                                                                                                                                     tataATTAAAGGTAATGATTTGATCATTGTTGGCACCAATGGTAACGGC‑AACCAGCACAATGGCTAAGGTAATGATCAATCCCAAAATATATTTAATCATTTGTTTCTCCTTTAGATAAATGTGTGAGCATTATGCCAAAAATACATATAAAAAACGACACTTAGAAGTGCCGTTATGAAAAGTTAGTTAGTAAAAAATTATGCTGCATAAACGTCGACACGTTCTCTTAATTCTTTGCCTGCCTTAAAat                                                                          <  1:729538/251‑1 (MQ=255)                                                                                                                                                                       taATTAAAGGTAATGATTTGATCATTGTTGGCACCAATGGTAACGGC‑AACCAGCACAATGGCTAAGGTAATGATCAATCCCAAAATATATTTAATCATTTGTTTCTCCTTTAGATAAATGTGTGAGCATTATGCCAAAAATACATATAAAAAACGACACTTAGAAGTGCCGTTATGAAAAGTTAGTTAGTAAAAAATTATGCTGCATAAACGTCGACACGTTCTCTTAATTCTTTGCCTGCCTTAAAatat                                                                        >  1:17892/1‑251 (MQ=255)                                                                                                                                                                                                 gttgGCACCAATGGTAACGGC‑AACCAGCACAATGGCTAAGGTAATGATCAATCCCAAAATATATTTAATCa                                                                                                                                                                                                                                  <  1:385840/71‑1 (MQ=255)                                                                                                                                                                                                 gttgGCACCAATGGTAACGGC‑AACCAGCACAATGGCTAAGGTAATGATCAATCCCAAAATATATTTAATCATTTGTTTCTCCTTTAGATAAATGTGTGAGCATTATGCCAAAAATACATATAAAAAACGACACTTAGAAGTGCCGTTATGAAAagttag                                                                                                                                          <  1:122028/159‑1 (MQ=255)                                                                                                                                                                                                  ttgGCACCAATGGTAACGGC‑AACCAGCACAATGGCTAAGGTAATGATCAATCCCAAAATATATTTAATCATTTGTTTCTCCTTTAGATAAATGTGTGAGCATTATGCCAAAAATACATATAAAAAACGACACTTAGAAGTGCCGTTATGAAAAGTTAGTTAGTAAAAAATTATGCTGCATAAACGTCGACACGTTCTCTTAATTCTTTGCCTGCCTTAAAATATGGCACACATTTGGCATCCAATTTTACt                                             <  1:419749/251‑1 (MQ=255)                                                                                                                                                                                                     gCACCAATGGTAACGGC‑AACCAGCACAATGGCTAAGGTAATGATCAATCCCAAAATATATTTAATCATTTGTTTCTCCTTTAGATAAATGTGTGAGCATTATGCCAAAAATACATATAAAAAACGACACTTATAAGTGCCGTTATGAAAAGTTAGTTAGTAAAAAATTATGCTGCATAAGCGTCGACACGTTCTCTTAATTc                                                                                           <  1:613550/202‑1 (MQ=255)                                                                                                                                                                                                                    gC‑AACCAGCACAATGGCTAAGGTAATGATCAATCCCAAAATATATTTAATCATTTGTTTCTCCTTTAGATAAATGTGTGAGCATTATGCCAAAAATACATATAAAAAACGACACTTAGAAGTGCCGTTATGAAAAGTTAGTTAGTAAAAAATTATGCTGCATAAACGTCGACACGTTCTCTTAATTCTTTGCCTGCCTTAAAATATGGCACACATTTGGCATCCAATTTTACTTGTTCGCCGGTTTTAGGa                           >  1:92894/1‑251 (MQ=255)                                                                                                                                                                                                                       aaCCAGCACAATGGCTAAGGTAATGATCAATCCCAAAATATATTTAATCATTTGTTTCTCCTTTAGATAAATGTGTGAGCATTATGCCAAAAATACATATAAAAAACGACACTTAGAAGTGCCGTTATGAAAAGTTAGTTAGTAAAAAATTATGCTGCATAAACGTCGACACGTTCTCTTAATTc                                                                                           <  1:298805/185‑1 (MQ=255)                                                                                                                                                                                                                       aaCCAGCACAATGGCTAAGGTAATGATCAATCCCAAAATATATTTAATCATTTGTTTCTCCTTTAGATAAATGTGTGAGCATTATGCCAAAAATACATATAAAAAACGACACTTAGAAGTGCCGTTATGAAAAGTTAGTTAGTAAAAAATTATGCTGCATAAACGTCGACACGTTCTCTTAATTCTTTGCCTGCCTTAAAATATGGCACACATTTGGCATCCAATTTTACTTGTTCGCCGGTTTTAGGAtt                         >  1:56575/1‑251 (MQ=255)                                                                                                                                                                                                                             cacaATGGCTAAGGTAATGATCAATCCCAAAATATATTTAATCATTTGTTTCTCCTTTAGATAAATGTGTGAGCATTATGCCAAAAATACATATAAAAAACGACACTTAGAAGTGCCGTTATGAAAAGTTAGTTAGTAAAAAATTATGCTGCATAAACGTCGACACGTTCTCTTAATTCTTTGCCTGCCTTAAAATATGGCACACATTTGGCATCCAATTTTACTTGTTCGCCGGTTTTAGGATTACGTcc                   <  1:397930/251‑1 (MQ=255)                                                                                                                                                                                                                                   ggCTAAGGTAATGATCAATCCCAAAATATATTTAATCATTTGTTTCTCCTTTAGATAAATGTGTGAGCATTATGCCAAAAATACATATAAAAAACGACACTTAGAAGTGCCGTTATGAAAAGTTAGTTAGTAAAAAATTATGCTGCATAAACGTCGACACGTTCTCTTAATTCTTTGCCTGCCTTAAAATATGGCACACATTTGGCATCCAATTTTACTTGTTCGCCGGTTTTAGGATTACGTCCGATgcg             >  1:125741/1‑251 (MQ=255)                                                                                                                                                                                                                                          gTAATGATCAATCCCAAAATATATTTAATCATTTGTTTCTCCTTTAGATAAATGTGTGAGCATTATGCCAAAAATACATATAAAAAACGACACTTAGAAGTGCCGTTATGAAAAGTTAGTTAGTAAAAAATTATGCTGCATAAACGTCGACACGTTCTCTTAATTCTTTGCCTGCCTTAAAATATGGCACACATTTGGCATCCAATTTTACTTGTTCGCCGGTTTTAGGATTACGTCCGATGCGCGGTTgg      >  1:646424/1‑251 (MQ=255)                                                                                                                                                                                                                                          gTAATGATCAATCCCAAAATATATTTAATCATTTGTTTCTCCTTTAGATAAATGTGTGAGCATTATGCCAAAAATACATATAAAAAACGACACTTAGAAGTGCCGTTATGAAAAGTTAGTTAGTAAAAAATTATGCTGCATAAACGTCGACACGTTCTCTTAATTCTTTGCCTGCCTTAAAATATGGCACACATTTGGCATCCAATTTTACTTGTTCGCCGGTTTTAGGATTACGTCCGATGCGCGGTTgg      >  1:791839/1‑251 (MQ=255)                                                                                                                                                                                                                                              tGATCAATCCCAAAATATATTTAATCATTTGTTTCTCCTTTAGATAAATGTGTGAGCATTATGCCAAAAATACATATAAAAAACGACACTTAGAAGTGCCGTTATGAAAAGTTAGTTAGTAAAAAATTATGCTGCATAAACGTCGACACGTTCTCTTAATTCTTTGCCTGCCTTAAAATATGGCACACATTTGGCATCCAATTTTACTTGTTCGCCGGTTTTAGGATTACGTCCGATGCGCGGTTGGCGAc  <  1:209888/251‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    TTCGTTGATTTGTTGCGTTTGGCGTTTGATTTGGCGATTCAGCGCAATATTTTTTAACTTTAATTTGAGATAGAAGAATCCGGTGATAAACCACCCTAGAATTAATCCCAGGCCGAATAAAATAGCAACTAAGGTAGAAAGCTGGAATTGGCTTTGCGCCACAATATAATTAAAGGTAATGATTTGATCATTGTTGGCACCAATGGTAACGGC‑AACCAGCACAATGGCTAAGGTAATGATCAACCCCAAAATATATTTAATCATTTGTTTCTCCTTTAGATAAATGTGTGAGCATTATGCCAAAAATACATATAAAAAACGACACTTAGAAGTGCCGTTATGAAAAGTTAGTTAGTAAAAAATTATGCTGCATAAACGTCGACACGTTCTCTTAATTCTTTGCCTGCCTTAAAATATGGCACACATTTGGCATCCAATTTTACTTGTTCGCCGGTTTTAGGATTACGTCCGATGCGCGGTTGGCGAC  >  1/1770140‑1770626

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