Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,608,895	C→A	Q34H (CAG→CAT)	tusA ←	sulfur transfer protein TusA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,608,895	0	C	A	100.0%	156.5 / NA	45	Q34H (CAG→CAT)	tusA	sulfur transfer protein TusA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (29/16);  total (29/16)

CCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGATCGTCGGCGATAATCAGCAACGTTTCGCCAGGCTGCATATTGCGCACGGTTTTGCGCACC‑‑‑‑ATCATCACCGGTTCCGGGCAGCGCAGGCCAAGCGCGTCGAGTGTGTGGTCAGGGCTGGAAAAGAGATCGGTCATTTTCTTCTCATCACTTAAAAAAACGGCGCT  >  NC_000913/3608758‑3609026                                                                                                                                          |                                                                                                                                        ccTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGTCGGATCGTCGGCGATAATCAGCAACGTTTCGCCAGGATg                                                                                                                                       >  2:697682/1‑140 (MQ=255) ccTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGATCGTCGGCGATAATCAGCAACGTTTCGCCAGGATg                                                                                                                                       >  1:695901/1‑140 (MQ=255)  cTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGATCGTCGGCGATAATCAGCAACGTTTCGCCAGGATGc                                                                                                                                      >  1:45388/1‑140 (MQ=255)    ttACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGATCGTCGGCGATAATCAGCAACGTTTCGCCAGGATGCat                                                                                                                                    >  1:946391/1‑140 (MQ=255)           tCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGATCGTCGGCGATAATCAGCAACGTTTCGCCAGGATGCATATTgcgc                                                                                                                             >  1:203201/1‑140 (MQ=255)                    gATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGATCGTCGGCGATAATCAGCAACGTTTCGCCAGGATGCATATTGCGCACGGTTTTg                                                                                                                    >  1:800922/1‑140 (MQ=255)                    gATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGATCGTCGGCGATAATCAGCAACGTTTCGCCAGGATGCATATTGCGCACGGTTTTg                                                                                                                    <  2:800922/140‑1 (MQ=255)                    gATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGATCGTCGGCGATAATCAGCAACGTTTCGCCAGGATGCATATTGCGCACGGTTTTg                                                                                                                    >  2:568011/1‑140 (MQ=255)                    gATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGATCGTCGGCGATAATCAGCAACGTTTCGCCAGGATGCATATTGCGCACGGTTTTg                                                                                                                    >  1:695207/1‑140 (MQ=255)                        aGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGATCGTCGGCGATAATCAGCAACGTTTCGCCAGGATGCATATTGCGCACGGTTTTGCGCa                                                                                                                >  2:786397/1‑140 (MQ=255)                        aGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGATCGTCGGCGATAATCAGCAACGTTTCGCCAGGATGCATATTGCGCACGGTTTTGCGCa                                                                                                                >  1:748713/1‑140 (MQ=255)                        aGGCAGTCCATCCGGCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGACCGTCGGCGATAATCAGCAACGTTTCGCCAGGATGCATATTGCGCACGGTTTTGCGCa                                                                                                                <  2:695156/140‑1 (MQ=255)                                        tctTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGATCGTCGGCGATAATCAGCAACGTTTCGCCAGGATGCATATTGCGCACGGTTTTGCGCACC‑‑‑‑ATCATCACCGGTTc                                                                                            >  2:690425/1‑140 (MQ=255)                                            tAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGATCGTCGGCGATAATCAGCAACGTTTCGCCAGGATGCATATTGCGCACGGTTTTGCGCACC‑‑‑‑ATCATCACCGGTTCCggg                                                                                        <  2:700424/140‑1 (MQ=255)                                                  ccAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGATCGTCGGCGATAATCAGCAACGTTTCGCCAGGATGCATATTGCGCACGGTTTTGCGCACC‑‑‑‑ATCATCACCGGTTCCGGGCAgcg                                                                                   >  1:385775/1‑139 (MQ=255)                                                  ccAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGATCGTCGGCGATAATCAGCAACGTTTCGCCAGGATGCATATTGCGCACGGTTTTGCGCACC‑‑‑‑ATCATCACCGGTTCCGGGCAgcg                                                                                   <  2:385775/139‑1 (MQ=255)                                                    aGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGATCGTCGGCGATAATCAGCAACGTTTCGCCAGGATGCATATTGCGCACGGTTTTGCGCACC‑‑‑‑ATCATCACCGGTTCCGGGCAg                                                                                     >  1:835660/1‑135 (MQ=255)                                                    aGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGATCGTCGGCGATAATCAGCAACGTTTCGCCAGGATGCATATTGCGCACGGTTTTGCGCACC‑‑‑‑ATCATCACCGGTTCCGGGCAg                                                                                     <  2:835660/135‑1 (MQ=255)                                                    aGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGATCGTCGGCGATAATCAGCAACGTTTCGCCAGGATGCATATTGCGCACGGTTTTGCGCACC‑‑‑‑ATCATCACCGGTTCCGGGCAGCGCAg                                                                                <  2:695207/140‑1 (MQ=255)                                                        cGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGATCGTCGGCGATAATCAGCAACGTTTCGCCAGGATGCATATTGCGCACGGTTTTGCGCACC‑‑‑‑ATCATCACCGGTTCCGGGCAGCGCAGGCCa                                                                            <  2:748713/140‑1 (MQ=255)                                                                 tAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGATCGTCGGCGATAATCAGCAACGTTTCGCCAGGATGCATATTGCGCACGGTTTTGCGCACC‑‑‑‑ATCATCACCGGTTCCGGGCAGCGCAGGCCAAGCGCGTCg                                                                   >  1:332803/1‑140 (MQ=255)                                                                    aGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGATCGTCGGCGATAATCAGCAACGTTTCGCCAGGATGCATATTGCGCACGGTTTTGCGCACC‑‑‑‑ATCATCACCGGTTCCGGGCAGCGCAGGCCAAGCGCGTCGAgt                                                                >  2:203832/1‑140 (MQ=255)                                                                         cAAAACCCAGGAATATCGCGGGTAGTGGCCGGATCGTCGGCGATAATCAGCAACGTTTCGCCAGGATGCATATTGCGCACGGTTTTGCGCACC‑‑‑‑ATCATCACCGGTTCCGGGCAGCGCAGGCCAAGCGCGTCGAGtgtgtg                                                           >  1:146681/1‑140 (MQ=255)                                                                         cAAAACCCAGGAATATCGCGGGTAGTGGCCGGATCGTCGGCGATAATCAGCAACGTTTCGCCAGGATGCATATTGCGCACGGTTTTGCGCACC‑‑‑‑ATCATCACCGGTTCCGGGCAGCGCAGGCCAAGCGCGTCGAGtgtgtg                                                           >  1:121791/1‑140 (MQ=255)                                                                         cAAAACCCAGGAATATCGCGGGTAGTGGCCGGATCGTCGGCGATAATCAGCAACGTTTCGCCAGGATGCATATTGCGCACGGTTTTGCGCACC‑‑‑‑ATCATCACCGGTTCCGGGCAGCGCAGGCCAAGCGCGTCGAGtgtgtg                                                           >  1:160134/1‑140 (MQ=255)                                                                            aaCCCAGGAATATCGCGGGTAGTGGCCGGATCGTCGGCGATAATCAGCAACGTTTCGCCAGGATGCATATTGCGCACGGTTTTGCGCACC‑‑‑‑ATCATCACCGGTTCCGGGCAGCGCAGGCCAAGCGCGTCGAGTGTGTGGTc                                                        <  1:203832/140‑1 (MQ=255)                                                                            aaCCCAGGAATATCGCGGGTAGTGGCCGGATCGTCGGCGATAATCAGCAACGTTTCGCCAGGATGCATATTGCGCACGGTTTTGCGCACC‑‑‑‑ATCATCACCGGTTCCGGGCAGCGCAGGCCAAGCGCGTCGAGTGTGTGGTc                                                        >  1:222470/1‑140 (MQ=255)                                                                                cAGGAATATCGCGGGTAGTGGCCGGATCGTCGGCGATAATCAGCAACGTTTCGCCAGGATGCATATTGCGCACGGTTTTGCGCACCACAGATCATCACCGGTTCCGGGCAGCGCAGGCCAAGCGCGTCGAGTGTGTGGTc                                                        >  2:918293/1‑140 (MQ=255)                                                                                     atatCGCGGGTAGTGGCCGGATCGTCGGCGATAATCAGCAACGTTTCGCCAGGATGCATATTGCGCACGGTTTTGCGCACC‑‑‑‑ATCATCACCGGTTCCGGGCAGCGCAGGCCAAGCGCGTCGAGTGTGTGGTCAGGGCTGGa                                               <  2:222470/140‑1 (MQ=255)                                                                                         cgcgGGTAGTGGCCGGATCGTCGGCGATAATCAGCAACGTTTCGCCAGGATGCATATTGCGCACGGTTTTGCGCACC‑‑‑‑ATCATCACCGGTTCCGGGCAGCGCAGGCCAAGCGCGTCGAGTGTGTGGTCAGGGCTGGAAAag                                           >  1:628255/1‑140 (MQ=255)                                                                                                aGTGGCCGGATCGTCGGCGATAATCAGCAACGTTTCGCCAGGATGCATATTGCGCACGGTTTTGCGCACC‑‑‑‑ATCATCACCGGTTCCGGGCAGCGCAGGCCAAGCGCGTCGAGTGTGTGGTCAGGGCTGGAAAAGAGATCgg                                    <  2:946391/140‑1 (MQ=255)                                                                                                aGTGGCCGGATCGTCGGCGATAATCAGCAACGTTTCGCCAGGATGCATATTGCGCACGGTTTTGCGCACC‑‑‑‑ATCATCACCGGTTCCGGGCAGCGCAGGCCAAGCGCGTCGAGTGTGTGGTCAGGGCTGGAAAAGACATCgg                                    <  1:568011/140‑1 (MQ=255)                                                                                                  tGGCCGGATCGTCGGCGATAATCAGCAACGTTTCGCCAGGATGCATATTGCGCACGGTTTTGCGCACC‑‑‑‑ATCATCACCGGTTCCGGGCAGCGCAGGCCAAGCGCGTCGAGTGTGTGGTCAGGGCTGGAAAAGAGATCGGTc                                  >  2:231904/1‑140 (MQ=255)                                                                                                    gCCGGATCGTCGGCGATAATCAGCAACGTTTCGCCAGGATGCATATTGCGCACGGTTTTGCGCACC‑‑‑‑ATCATCACCGGTTCCGGGCAGCGCAGGCCAAGCGCGTCGAGTGTGTGGTCAGGGCTGGAAAAGAGAt                                       >  2:901090/1‑133 (MQ=255)                                                                                                    gCCGGATCGTCGGCGATAATCAGCAACGTTTCGCCAGGATGCATATTGCGCACGGTTTTGCGCACC‑‑‑‑ATCATCACCGGTTCCGGGCAGCGCAGGCCAAGCGCGTCGAGTGTGTGGTCAGGGCTGGAAAAGAGAt                                       <  1:901090/133‑1 (MQ=255)                                                                                                    gCCGGATCGTCGGCGATAATCAGCAACGTTTCGCCAGGATGCATATTGCGCACGGTTTTGCGCACC‑‑‑‑ATCATCACCGGTTCCGGGCAGCGCAGGCCAAGCGCGTCGAGTGTGTGGTCAGGGCTGGAAAAGAGATCGGTCAt                                >  1:545662/1‑140 (MQ=255)                                                                                                        gATCGTCGGCGATAATCAGCAACGTCTCGCCAGGATGCATATTGCGCACGGTTTTGCGCACC‑‑‑‑ATCATCACCGGTTCCGGGCAGCGCAGGCCAAGCGCGTCGAGTGTGTGGTCAGGGCTGGAAAAGAGATCGGTCATTttc                            >  1:624671/1‑140 (MQ=255)                                                                                                              cgGCGATAATCAGCAACGTTTCGCCAGGATGCATATTGCGCACGGTTTTGCGCACC‑‑‑‑ATCATCACCGGTTCCGGGCAGCGCAGGCCAAGCGCGTCGAGTGTGTGGTCAGGGCTGGAAAAGAGAt                                       <  1:339357/123‑1 (MQ=255)                                                                                                              cgGCGATAATCAGCAACGTTTCGCCAGGATGCATATTGCGCACGGTTTTGCGCACC‑‑‑‑ATCATCACCGGTTCCGGGCAGCGCAGGCCAAGCGCGTCGAGTGTGTGGTCAGGGCTGGAAAAGAGAt                                       >  2:339357/1‑123 (MQ=255)                                                                                                               ggCGATAATCAGCAACGTTTCGCCAGGATGCATATTGCGCACGGTTTTGCGCACC‑‑‑‑ATCATCACCGGTTCCGGGCAGCGCAGGCCAAGCGCGTCGAGTGTGTGGTCAGGGCTGGAAAAGAGATCGGTCATTTTCTTCtcat                     <  2:45388/140‑1 (MQ=255)                                                                                                                           cAACGTTTCGCCAGGATGCATATTGCGCACGGTTTTGCGCACC‑‑‑‑ATCATCACCGGTTCCGGGCAGCGCAGGCCAAGCGCGTCGAGTGTGTGGTCAGGGCTGGAAAAGAGATCGGTCATTTTCTTCTCATCACTTaaaaaaa         >  1:564730/1‑140 (MQ=255)                                                                                                                              cGTTTCGCCAGGATGCATATTGCGCACGGTTTTGCGCACC‑‑‑‑ATCATCACCGGTTCCGGGCCGCGCAGGCCAAGCGCGTCGAGTGTGTGGTCAGGGCTGGAAAAGAGATCGGTCATTTTCTTCTCATCACTTAAAAAAACgg      <  1:787610/140‑1 (MQ=255)                                                                                                                              cGTTTCGCCAGGATGCATATTGCGCACGGTTTTGCGCACC‑‑‑‑ATCATCACCGGTTCCGGGCAGCGCAGGCCAAGCGCGTCGAGTGTGTGGTCAGGGCTGGAAAAGAGATCGGTCATTTTCTTCTCATCACTTAAAAAAACgg      >  2:527667/1‑140 (MQ=255)                                                                                                                              cGTTTCGCCAGGATGCATATTGCGCACGGTTTTGCGCACC‑‑‑‑ATCATCACCGGTTCCGGGCAGCGCAGGCCAAGCGCGTCGAGTGTGTGGTCAGGGCTGGAAAAGAGATCGGTCATTTTCTTCTCATCACTTAAAAAAACgg      >  2:383158/1‑140 (MQ=255)                                                                                                                                  tCGCCAGGATGCATATTGCGCACGGTTTTGCGCACC‑‑‑‑ATCATCACCGGTTCCGGGCAGCGCAGGCCAAGCGCGTCGAGTGTGTGGTCAGGGCTGGAAAAGAGATCGGTCATTTTCTTCTCATCACTTAAAAAAACGGCGCt  <  2:695901/140‑1 (MQ=255)                                                                                                                                          |                                                                                                                                        CCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGATCGTCGGCGATAATCAGCAACGTTTCGCCAGGCTGCATATTGCGCACGGTTTTGCGCACC‑‑‑‑ATCATCACCGGTTCCGGGCAGCGCAGGCCAAGCGCGTCGAGTGTGTGGTCAGGGCTGGAAAAGAGATCGGTCATTTTCTTCTCATCACTTAAAAAAACGGCGCT  >  NC_000913/3608758‑3609026

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