Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 83,296 T→C E78E (GAA→GAG)  leuA ← 2‑isopropylmalate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091383,2960TC96.2% 78.4 / ‑3.6 26E78E (GAA→GAG) leuA2‑isopropylmalate synthase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (16/9);  total (16/10)
Fisher's exact test for biased strand distribution p-value = 3.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.95e-01

CGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTCACCATATAGATAGCGCGTTCGATCACCTCGTCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTTTCCACGCAGCGAGCTAACGCACATACGCGGCTGTTTTTAACCTGGCGGGCGATGGTTTGCACCGATTCAAAATCGCCCGGCGAAGAGACGGGG‑AAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTT  >  NC_000913/83056‑83533
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                              
cgccgcTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTCACCATATAGATAGCGCGTTCGATCACTTCATCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGc                                                                                                                                                                                                                                      <  1:119428/251‑1 (MQ=255)
           ccACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTCACCATATAGATAGCGCGTTCGATCACCTCGTCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAATGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGTCCGCCACGTCGATATCTTTTTCCACGCAGCGAGCTAACGca                                                                                                                                                                                                                           <  1:269945/251‑1 (MQ=255)
                        aTCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTCACCATATAGATAGCGCGTTCGATCACTTCATCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATCTCTTTCTCCACGCAGCGAGCTAACGCGCATAAGCGGCTGt                                                                                                                                                                                                              >  2:54268/1‑251 (MQ=255)
                           ggCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTCACCATATAGATAGCGCGTTCGATCACTTCATCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGtttt                                                                                                                                                                                                           >  2:69436/1‑251 (MQ=255)
                                       aCGCCCGGCATCTTCGCAAGAAAATTCAACATCAGCGGTGTAATTACGGGCGCGTTTCACCATATAGATAGCGCGTTCGATCACTTCATCCAGCGTGCTGCGCAGGTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGGGCTAACGCGCATACGCGGCTGTTTTTAACCTggcgag                                                                                                                                                                                               <  2:268652/251‑3 (MQ=255)
                                           ccGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTCACCATATAGATAGCGCGTTCGATCACTTCATCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGAt                                                                                                                                                                                           <  2:161496/251‑1 (MQ=255)
                                                              aTTCAACATCATCGGTGTAATTACGGGCGCGTTTCACCATATAGATAGCGCGTTCGATCACTTCATCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCaaa                                                                                                                                                                         >  1:128768/1‑250 (MQ=255)
                                                                              gtAATTACGGGCGCGTTTCACCATATAGATAGCGCGTTCGATCACTTCATCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAagaa                                                                                                                                                         <  1:96336/250‑2 (MQ=255)
                                                                                       ggCGCGTTTCACCATATAGATAGCGCGTTCGATCACTTAATCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATAGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGG‑AAAc                                                                                                                                              >  1:115023/1‑251 (MQ=255)
                                                                                          gcgTTTCCCCATATAGATAGCGCGTTCGATCACTTCATCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCAGACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGG‑AAACCGa                                                                                                                                           <  2:128768/251‑1 (MQ=255)
                                                                                             tttCACCATATAGATAGCGCGTTCGATCACTTCATCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGG‑AAACCGACtt                                                                                                                                        >  1:81209/1‑251 (MQ=255)
                                                                                                cACCATATAGATAGCGCGTTCGATCACTTCATCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGG‑AAACCGACTTcc                                                                                                                                      >  1:95654/1‑250 (MQ=255)
                                                                                                   cATATAGATAGCGCGTTCGATCACTTCATCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGG‑AAACCGACTTCCATc                                                                                                                                   >  1:217235/1‑250 (MQ=255)
                                                                                                            agCGCGTTCGATCACTTCATCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGG‑AAACCGACTTCCATCACGTCAacac                                                                                                                         >  1:186197/1‑251 (MQ=255)
                                                                                                                cgTTCGATCACTTCATCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAGGGTATGAATACGGAAGGCGTCGGCAACTTTCAAGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCATCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTGGCAATGATTCAAAATCACCCGGCGAAGAAACGGGG‑AAACCGACTTCCATCACGTCAACACCCAt                                                                                                                     >  1:172/1‑251 (MQ=255)
                                                                                                                    cGATCACTTCATCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGG‑AAACCGACTTCCATCACGTCAACACCCATACGc                                                                                                                 >  2:230845/1‑251 (MQ=255)
                                                                                                                                           cgcAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGG‑AAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGc                                                                                            >  1:128072/1‑249 (MQ=255)
                                                                                                                                                       gCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAATGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTTATTCAAAATCACCCGGCGAAGAAACGTGG‑AAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTGTTTCTTTcac                                                                              <  1:54268/251‑1 (MQ=255)
                                                                                                                                                                     cGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGG‑AAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCt                                                                >  1:67116/1‑251 (MQ=255)
                                                                                                                                                                      gAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGG‑AAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTg                                                               >  2:20637/1‑251 (MQ=255)
                                                                                                                                                                               aTAAAGGTATGAATACGTAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGG‑AAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCACTTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTg                                                      >  2:144427/1‑251 (MQ=255)
                                                                                                                                                                                                     cTTCGGCGACTCTCAGGGATTCGGCCGCCACGTCGCTATCTTTCTCCACGCAGCGCGATAAAGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGG‑AAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCCCACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGt                                 >  2:267642/1‑250 (MQ=255)
                                                                                                                                                                                                       tCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATTCGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGG‑AAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTa                                <  1:81067/249‑1 (MQ=255)
                                                                                                                                                                                                                     ggATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGG‑AAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACtt                <  1:69436/251‑1 (MQ=255)
                                                                                                                                                                                                                      gATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGG‑AAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTg               <  2:52687/251‑1 (MQ=255)
                                                                                                                                                                                                                                    gTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGGAAACCGACTTCCATCACGTCACCACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGtt  >  1:152111/1‑251 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                              
CGCCGCTTCGACCACTCGCGCCAGATCGGCAATGGGTGTACGCCCGGCATCTTCGCAAGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTCACCATATAGATAGCGCGTTCGATCACCTCGTCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTTTCCACGCAGCGAGCTAACGCACATACGCGGCTGTTTTTAACCTGGCGGGCGATGGTTTGCACCGATTCAAAATCGCCCGGCGAAGAGACGGGG‑AAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTT  >  NC_000913/83056‑83533

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