Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,598,557	A→G	*368Q (TAA→CAA)	livJ ←	branched‑chain amino acid ABC transporter periplasmic binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,598,557	0	A	G	100.0%	141.4 / NA	40	*368Q (TAA→CAA)	livJ	branched‑chain amino acid ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (23/17);  total (23/17)

TACAATGTGTGTTGCAATTTACTGATTTCTTTGGATCTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAAGCAACGTACTCGATATTAGCAATTTGGCGGCAACCCAAAGTTGCCGATTAATGATTACTTCGCATCGGTCGCCGTGCCGTTGGCGTGCCAGTCAAATACGCCGAACTCAAAGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCACGGAG  >  NC_000913/3598423‑3598669                                                                                                                                       |                                                                                                                 tACAATGTGTGTTGCAATTTACTGATTTCTTTGGATCTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAAGCAACGTACTCGATATTAGCAATTTGGCGGCAACCCAAAGTTGCCGATTAATGATTGCTTc                                                                                                              >  2:563052/1‑139 (MQ=255)                              tttGGATCTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAAGCAACGTACTCGATATTAGCAATTTGGCGGCAACCCAAAGTTGCCGATTAATGATTGCTTCGCATCGGTCGCCGTGCCGTTGGCGTGCCa                                                                                 >  1:1083829/1‑139 (MQ=255)                              tttGGATCTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAAGCAACGTACTCGATATTAGCAATTTGGCGGCAACCCAAAGTTGCCGATTAATGATTGCTTCGCATCGGTCGCCGTGCCGTTGGCGTGCCa                                                                                 >  2:1016785/1‑139 (MQ=255)                                        gTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAAGCAACGTACTCGATATTAGCAATTTGGCGGCAACCCAAAGTTGCCGATTAATGATTGCTTCGCATCGGTCGCCGTGCCGTTGGCGTGCCAGTCAAATACg                                                                       <  1:563052/139‑1 (MQ=255)                                              cGGATAAGGCGTTTACGCCGCATCCGGCATGAAGCAACGTACTCGATATTAGCAATTTGGCGGCAACCCAAAGTTGCCGATTAATGATTGCTTCGCATCgg                                                                                                       <  1:1053182/101‑1 (MQ=255)                                              cGGATAAGGCGTTTACGCCGCATCCGGCATGAAGCAACGTACTCGATATTAGCAATTTGGCGGCAACCCAAAGTTGCCGATTAATGATTGCTTCGCATCgg                                                                                                       >  2:1053182/1‑101 (MQ=255)                                                gATAAGGCGTTTACGCCGCATCCGGCATGAAGCAACGTACTCGATATTAGCAATTTGGCGGCAACCCAAAGTTGCCGATTAATGATTGCTTCGCATCGGTCGCCGTGCCGTTGGCGTGCCAGTCAAATACGCCGAACTc                                                               >  1:72061/1‑139 (MQ=255)                                                 aTAAGGCGTTTACGCCGCATCCGGCATGAAGCAACGTACTCGATATTAGCAATTTGGCGGCAACCCAAAGTTGCCGATTAATGATTGCTTCGCATCGGTCGCCGTGCCGTTGGCGTGCCAGTCAAATACGCCGAACTCa                                                              >  2:480050/1‑139 (MQ=255)                                                 aTAAGGCGTTTACGCCGCATCCGGCATGAAGCAACGTACTCGATATTAGCAATTTGGCGGCAACCCAAAGTTGCCGATTAATGATTGCTTCGCATCGGTCGCCGTGCCGTTGGCGTGCCAGTCAAATACGCCGAACTCa                                                              <  1:1016785/139‑1 (MQ=255)                                                          ttACGCCGCATCCGGCATGAAGCAACGTACTCGATATTAGCAATTTGGCGGCAACCCAAAGTTGCCGATTAATGATTGCTTCGCATCGGTCGCCGTGCCGTTGGCGTGCCAGTCAAATACGCCGAACTCAAAGCCTTTc                                                     >  2:528306/1‑139 (MQ=255)                                                           tACGCCGCATCCGGCATGAAGCAACGTACTCGATATTAGCAATTTGGCGGCAACCCAAAGTTGCCGATTAATGATTGCTTCGCATCGGTCGCCGTGCCGTTGGCGTGCCAGTCAAATACGCCGAACTCaaa                                                            >  2:222828/1‑131 (MQ=255)                                                           tACGCCGCATCCGGCATGAAGCAACGTACTCGATATTAGCAATTTGGCGGCAACCCAAAGTTGCCGATTAATGATTGCTTCGCATCGGTCGCCGTGCCGTTGGCGTGCCAGTCAAATACGCCGAACTCaaa                                                            <  1:222828/131‑1 (MQ=255)                                                              gccgcATCCGGCATGAAGCAACGTACTCGATATTAGCAATTTGGCGGCAACCCAAAGTTGCCGATTAATGATTGCTTCGCATCGGTCGCCGTGCCGTTGGCGTGCCAGTCAAATACGCCGAACTCAAAGCCTTTCAGAt                                                 >  2:31657/1‑139 (MQ=255)                                                                cgcATCCGGCATGAAGCAACGTACTCGATATTAGCAATTTGGCGGCAACCCAAAGTTGCCGATTAATGATTGCTTCGCATCGGTCGCCGTGCCGTTGGCGTGCCAGTCAAATACGCCGAACTCa                                                              <  1:568489/124‑1 (MQ=255)                                                                cgcATCCGGCATGAAGCAACGTACTCGATATTAGCAATTTGGCGGCAACCCAAAGTTGCCGATTAATGATTGCTTCGCATCGGTCGCCGTGCCGTTGGCGTGCCAGTCAAATACGCCGAACTCa                                                              >  2:568489/1‑124 (MQ=255)                                                                  cATCCGGCATGAAGCAACGTACTCGATATTAGCAATTTGGCGGCAACCCAAAGTTGCCGATTAATGATTGCTTCGCATCGGTCGCCGTGCCGTTGGCGTGCCAGTCAAATACGCCGAACTCAAAGCCTTTCAGATCGcc                                             <  1:480050/139‑1 (MQ=255)                                                                        gCATGAAGCAACGTACTCGATATTAGCAATTTGGCGGCAACCCAAAGTTGCCGATTAATGATTGCTTCGCATCGGTCGCCGTGCCGTTGGCGTGCCAGTCAAATACGCCg                                                                    >  2:120986/1‑110 (MQ=255)                                                                        gCATGAAGCAACGTACTCGATATTAGCAATTTGGCGGCAACCCAAAGTTGCCGATTAATGATTGCTTCGCATCGGTCGCCGTGCCGTTGGCGTGCCAGTCAAATACGCCg                                                                    <  1:120986/110‑1 (MQ=255)                                                                              aGCAACGTACTCGATATTAGCAATTTGGCGGCAACCCAAAGTTGCCGATTAATGATTGCTTCGCATCGGTCGCCGTGCCGTTGGCGTGCCAGTCAAATACGCCGAACTCAAAGCCTTTCAGATCGCCTTTCTCATCCCa                                 >  1:367173/1‑139 (MQ=255)                                                                                  aCGTACTCGATATTAGCAATTTGGCGGCAACCCAAAGTTGCCGATTAATGATTGCTTCGCATCGGTCGCCGTGCCGTTGGCGTGCCAGTCAAATACGCCGAACTCAAAGCCTTTCAGATCGCCTTTCTCATCCCAGGTc                             >  1:820844/1‑139 (MQ=255)                                                                                       cTCGATATTAGCAATTTGGCGGCAACCCAAAGTTGCCGATTAATGATTGCTTCGCATCGGTCGCCGTGCCGTTGGCGTGCCAGTCAAATACGCCGAACTCAAAGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCgg                        <  1:528306/139‑1 (MQ=255)                                                                                        tCGATATTAGCAATTTGGCGGCAACCCAAAGTTGCCGATTAATGATTGCTTCGCATCGGTCGCCGTGCCGTTGGCGTGCCAGTCAAATACGCCGAACTCAAAGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGt                       >  1:920579/1‑139 (MQ=255)                                                                                              ttAGCAATTTGGCGGCAACCCAAAGTTGCCGATTAATGATTGCTTCGCATCGGTCGCCGTGCCGTTGGCGTGCCAGTCAAATACGCCGAACTCAAAGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCAtt                 <  2:820844/139‑1 (MQ=255)                                                                                              ttAGCAATTTGGCGGCAACCCAAAGTTGCCGATTAATGATTGCTTCGCATCGGTCGCCGTGCCGTTGGCGTGCCAGTCAAATACGCCGAACTCAAAGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCAtt                 <  2:367173/139‑1 (MQ=255)                                                                                                aGCAATTTGGCGGCAACCCAAAGTTGCCGATTAATGATTGCTTCGCATCGGTCGCCGTGCCGTTGGCGTGCCAGTCAAATACGCCGAACTCAAAGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTAc               >  1:347110/1‑139 (MQ=255)                                                                                                  cAATTTGGCGGCAACCCAAAGTTGCCGATTAATGATTGCTTCGCATc                                                                                                         >  1:860387/1‑47 (MQ=255)                                                                                                  cAATTTGGCGGCAACCCAAAGTTGCCGATTAATGATTGCTTCGCATc                                                                                                         <  2:860387/47‑1 (MQ=255)                                                                                                  cAATTTGGCGGCAACCCAAAGTTGCCGATTAATGATTGCTTCGCATCGGTCGCCGTGCCGTTGGCGTGCCAGTCAAATACGCCGAACTCAAAGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACgg             <  2:920579/139‑1 (MQ=255)                                                                                                    aTTTGGCGGCAACCCAAAGTTGCCGATTAATGATTGCTTCGCATCgg                                                                                                       >  2:50348/1‑47 (MQ=255)                                                                                                    aTTTGGCGGCAACCCAAAGTTGCCGATTAATGATTGCTTCGCATCgg                                                                                                       <  1:50348/47‑1 (MQ=255)                                                                                                      ttGGCGGCAACCCAAAGTTGCCGATTAATGATTGCTTCGCATCGGTCGCCGTGCCGTTGGCGTGCCAGTCAAATACGCCGAACt                                                                >  2:631724/1‑84 (MQ=255)                                                                                                      ttGGCGGCAACCCAAAGTTGCCGATTAATGATTGCTTCGCATCGGTCGCCGTGCCGTTGGCGTGCCAGTCAAATACGCCGAACt                                                                <  1:631724/84‑1 (MQ=255)                                                                                                      ttGGCGGCAACCCAAAGTTGCCGATTAATGATTGCTTCGCATCGGTCGCCGTGCCGTTGGCGTGCCAGTCAAATACGCCGAACTCaaa                                                            >  2:1091010/1‑88 (MQ=255)                                                                                                      ttGGCGGCAACCCAAAGTTGCCGATTAATGATTGCTTCGCATCGGTCGCCGTGCCGTTGGCGTGCCAGTCAAATACGCCGAACTCaaa                                                            <  1:1091010/88‑1 (MQ=255)                                                                                                      ttGGCGGCAACCCAAAGTTGCCGATTAATGATTGCTTCGCATCGGTCGCCGTGCCGTTGGCGTGCCAGTCAAATACGCCGAACTCAAAGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATc         >  2:1044699/1‑139 (MQ=255)                                                                                                          cggcAACCCAAAGTTGCCGATTAATGATTGCTTCGCATCGGTCGCCGTGCCGTTGGCGTg                                                                                    >  1:818470/1‑60 (MQ=255)                                                                                                          cggcAACCCAAAGTTGCCGATTAATGATTGCTTCGCATCGGTCGCCGTGCCGTTGGCGTg                                                                                    <  2:818470/60‑1 (MQ=255)                                                                                                          cggcAACCCAAAGTTGCCGATTAATGATTGCTTCGCATCGGTCGCCGTGCCGTTGGCGTGCCAGTCAAATACGCCGAACTCAAAGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCACg     >  2:522000/1‑139 (MQ=255)                                                                                                          cggcAACCCAAAGTTGCCGATTAATGATTGCTTCGCATCGGTCGCCGTGCCGTTGGCGTGCCAGTCAAATACGCCGAACTCAAAGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCACg     <  1:522000/139‑1 (MQ=255)                                                                                                             caagctaaaGTTGCCGATTAATGCTTGCTTCGTATCCCTCCCCGTGCCGTTGCCGTCCCTGTCACGTTCGCCAAACTCATAGCCTTTCCGATCGCCTTTCTTATCCCAGATCAGCGGTTCCGTTACGGTATCCACGGAg  >  2:498379/7‑139 (MQ=255)                                                                                                                                       |                                                                                                                 TACAATGTGTGTTGCAATTTACTGATTTCTTTGGATCTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAAGCAACGTACTCGATATTAGCAATTTGGCGGCAACCCAAAGTTGCCGATTAATGATTACTTCGCATCGGTCGCCGTGCCGTTGGCGTGCCAGTCAAATACGCCGAACTCAAAGCCTTTCAGATCGCCTTTCTCATCCCAGGTCAGCGGTCCCATTACGGTATCCACGGAG  >  NC_000913/3598423‑3598669

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