Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,148,972 C→T Q72* (CAG→TAG)  fabH → 3‑oxoacyl‑[acyl carrier protein] synthase 3

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,148,9720CT96.0% 67.7 / NA 25Q72* (CAG→TAG) fabH3‑oxoacyl‑[acyl carrier protein] synthase 3
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (11/13);  minor base G (0/1);  total (11/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGAAAAGTGACTGAGCGTACATGTATACGAAGATTATTGGTACTGGCAGCTATCTGCCCGAACAAGTGCGGACAAACGCCGATTTGGAAAAAATGGTGGACACCTCTGACGAGTGGATTGTCACTCGTACCGGTATCCGCGAACGCCACATTGCCGCGCCAAACGAAACCGTTTCAACCATGGGCTTTGAAGCGGCGACACGCGCAATTGAGATGGCGGGCATTGAGAAAGACCAGATTGGCCTGATCGTTGTGGCAACGACTTCTGCTACGCACGCTTTCCCGAGCGCAGCTTGTCAGATTCAAAGCATGTTGGGCATTAAAGGTTGCCCGGCATTTGACGTTGCAGCAGCCTGCGCAGGTTTCACCTATGCATTAAGCGTAGCCGATCAATACGTGAAATCTGGGGCGGTGAAGTATGCTCTGGTCGTCGGTTCCGATGTACTGGCGCGCACCTGCGATCCAACCGATCGTGGGACTA  >  NC_000913/1148739‑1149218
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                       
cGAAAAGTGACTGAGCGTACATGTATACGAAGATTATTGGTACTGGCAGCTATCTGCCCGAACAAGTGCGGACAAACGCCGATTTGGAAAAAATGGTGGACACCTCTGACGAGTGGATTGTCACTCGTACCGGTATCCGCGAACGCCACATTGCCGCGCCAAACGAAACCGTTTCAACCATGGGCTTTGAAGCGGCGACACGCGCAATTGAGATGGCGGGCATTGAGAAAGACTAGATTGGCCTGATCGtt                                                                                                                                                                                                                                        <  2:5908/251‑1 (MQ=255)
                              aGATTATTGGTACTGGCAGCTATCTGCCCGAACAAGTGCGGACAAACGCCGATTTGGAAAAAATGGTGGACACCTCTGACGAGTGGATTGTCACTCGTACCGGTATCCGCGAACGCCACATTGCCGCGCCAAACGAAACCGTTTCAACCATGGGCTTTGAAGCGGCGACACGCGCAATTGAGATGGCGGGCATTGAGAAAGACTAGATTGGCCTGATCGTTGTGGCAACTACTTCTGCTACGCACGCTTTc                                                                                                                                                                                                         <  2:202237/251‑1 (MQ=255)
                                  tattGGTACTGGCAGCTATCTGCCCGAACAAGTGCGGACAAACGCCGATTTGGAAAAAATGGTGGACACCTCTGACGAGTGGATTGTCCCTCGTACCGGTATCCGCGAACGCCACATTGCCGCGCCAAACGAAACCGTTTCAACCATGGGCTTTGAAGCGGCGACACGCGCAATTGAGATGGCGGGCATTGAGAAAGACTAGATTGGCCCGATCGTTGTGGCAACGACTTCTGCTAGCAACGCTTTCCCgg                                                                                                                                                                                                     >  2:271650/1‑250 (MQ=255)
                                             gCAGCTATCTGCCCGAACAAGTGCGGACAAACGCCGATTTGGAAAAAATGGTGGACCCCTCTGACGAGTGGATTGTCACTCGTACCGGTATCCGCGAACGCCACATTGCCGCGCCAAACGAAACCGTTTCAACCATGGGCTTTGAAGCGGCGACACGCGCAATTGAGATGGCGGGCATTGAGAAAGACTAGATTGGCCTGATCGTTGTGGCAACGACTTCTGCTACGCACGCTTTCCAGAGCGCAGCTTGt                                                                                                                                                                                          <  2:225470/251‑1 (MQ=255)
                                                gCTATCTGCCCAAACAAGCGCGGACAAACGCCGATTTGGAAAAAATGGTGGACACCTCTGACGAGTGGATTATCACTCGTACCGGTATCCGCGAACGCCACATTGCCGCGCCAAACGAAACCGATTCACCCCTGGGCTTTGAAGAGGCGACCCGCGCAATTGAAATGGCGGGCATTGAGCCAGACTAGATTGGCCTGATCGTTGTGGCAACGACTTCTGCTACGCACGATTTCCCGAGCGCAGCTTGTc                                                                                                                                                                                         >  2:130801/1‑249 (MQ=255)
                                                   aTCTGCCCGAACAAGTGCGGACAAACGCCGATTTGGAAAAAATGGTGGACACCTCTGACGAGTGGATTGTCACTCGTACCGGTATCCGCGAACGCCACATTGCCGCGCCAAACGAAACCGTTTCAACCATGGGCTTTGAAGCGGCGACACGCGCAATTGAGATGGCGGGCATTGAGAAAGACTAGATTGGCCTGATCGTTGTGGCAACGACTTCTGCTACGCACGCTTTCCCGAGCGCAGCTTGTCAGa                                                                                                                                                                                      <  2:9118/249‑1 (MQ=255)
                                                       gCCCGAACAAGTGCGGACAAACGCCGATTTGGAAAAAATGGTGGACACCTCTGACGAGTGGATTGTCACTCGTACCGGTATCCGCGAACGCCACATTGCCGCGCCAAACGAAACCGTTTCAACCATGGGCTTTGAAGCGGCGACACGCGCAATTGAGATGGCGGGCATTGAGAAAGACTAGATTGGCCTGATCGTTGTGGCAACGACTTCTGCTACGCACGCTTTCCCGAGCGCAGCTTGTCAGAtt                                                                                                                                                                                    <  1:68350/247‑1 (MQ=255)
                                                            aaCAAGTGCGGACAAACGCCGATTTGGAAAAAATGGTGGACACCTCTGACGAGTGGATTGTCACTCGTACCGGTATCCGCGAACGCCACATTGCCGCGCCAAACGAAACCGTTTCAACCATGGGCTTTGAAGCGGCGACACGCGCAATTGAGATGGCGGGCATTGAGAAAGACTAGATTGGCCTGATCGTTGTGGCAACGACTTCTGCTACGCACGCTTTCCCGAGCGCAGCTTGGCAGATTCAAAGCATg                                                                                                                                                                           >  1:118602/1‑251 (MQ=255)
                                                                                          caatggagtacacCTCTGACGAGTGGATTGTCACTCGTACCGATATCCGCGAACGCCACATTGCCGCGCCAAACGAAACCGTTTCAACCATGGGCTTTGAAGCGGCGACACGCGCAATTGAGATGGCGGGCATTGAGAAAGACTAGATTGGCCTGATCGTTGTGGCAACGACTTCTGCTACGCACGCTTTCCCGAGCGCAGCTTGTCAGATTCAAAGCATGTTGGGCATTCAAGGTAGCCCGGCAtttta                                                                                                                                              >  2:58236/10‑248 (MQ=255)
                                                                                          aaaTGGTGGACACCTCTGACGAGTGGATTGTCACTCGTACCGGTATCCGCGAACGCCACATTGCCGCGCCAAACGAAACCGTTTCAACCATGGGCTTTGAAGCGGCGACACGCGCAATTGAGATGGCGGGCATTGAGAAAGACTAGATTGGCCTGATCGTTGTGGCAACGACTTCTGCTACGCACGCTTTCCCGAGCGCAGCTTGTCAGATTCAAAGCATGTTGGGCATTAAAGGTTGCCCGGCATTTGAc                                                                                                                                             >  2:182112/1‑251 (MQ=255)
                                                                                                 ggACACCTCTGACGAGTGGATTGTCACTCGTACCGGTATCCGCGAACGCCACATTGCCGCGCCAAACGAAACCGTTTCAACCATGGGCTTTGAAGCGGCGACACGCGCAATTGAGATGGCGGGCATTGAGAAAGACTAGATTGGCCTGATCGTTGTGGCAACGACTTCTGCTACGCACGCTTTCCCGAGCGCAGCTTGTCAGATTCAAAGCATGTTGGGCATTAAAGGTTGCCCGGCATTTGACGTTgca                                                                                                                                       <  1:268169/250‑1 (MQ=255)
                                                                                                              gAGTGGATTGTCACTCGTACCGGTATCCGCGAACGCCACATTGCCGCGCCAAACGAAACCGTTTCAACCATGGGCTTTGAAGCGGCGACACGCGCAATTGAGAAGGCGGGCATTGAGAAAGACTAGATTGGCCTGATCGTTGTGGCAAAGACTTCTGCTACGCACGCTTTCCCGAGCGCCGCTTGTCAGATTCAAAGCATGTTGGGCATTAAAGGTTGCCCGGCATTTGACGTTGAAGCAGCCTGCGCAgg                                                                                                                         >  1:123304/1‑251 (MQ=255)
                                                                                                                             cGTACCGGTATCCGCGAACGCCACATTGCCGCGCCAAACGAAACCGTTTCAACCATGGGCTTTGAAGCGGCGACACGCGCAATTGAGATGGCGGGCATTGAGAAAGACTAGATTGGCCTGATCGTTGTGGCAACGACTTCTGCTACGCACGCTTTCCCGAGCGCAGCTTGTCAGATTCAAAGCATGTTGGGCATTAAAGGTTGCCCGGCATTTGACGTTGCAGCAGCCTGCGCAGGTTTCACCTATGCAtt                                                                                                          <  1:271650/251‑1 (MQ=255)
                                                                                                                                tccGGTATCCGCGAATGCCACATTGCCGCGCCAAACGAAACCGTTTCAACCATGGGCTTTGAAGCGGCGACACGCGCAATTGAGATGGCGGGCATTGAGAAAGACTAGATTGGCCTGATCGTTGTGGCAACGACGTCTGCTACGCACGCTTTCCCGAGCGCAGCTTGTCAGATTCAAAGCATGTTGGGCATTAAAGGTTGCCCGGCATTTGACGGTGCAGCAGCCTGCTCAGGTTTCACCTATGCAttca                                                                                                        <  2:239674/249‑3 (MQ=255)
                                                                                                                                       tCCGCGAACGCCACATTGCCGCGCCAAACGAAACCGTTTCCACCATGAGCTTTGAAGCGGCGACACGCGCAATTGAGCTGGCGGGCATTGAGAAAGACTAGATTGGCCTGATCGTTGTGGCAACGACTTCTGCTACGCACGCTTTCCCGAGCGCAGCTTGTCAGATTCAACGCATGTTGGGCATTAAAGGGTGCCCGGCATTTGACGTTGCAGCAGCCTGCGCAGGTTTCACCTATGCATTAAGCGTAGcc                                                                                                >  2:146407/1‑251 (MQ=255)
                                                                                                                                                 ccACAGTGCCGCGCCAAACGAAACCGTTTCAACCATGGGCTGTGAAGCGGCGACACGCGCAATTGAGATGGCGGGCATTGAGAAAGACTAGATTGGCCTGATCGTTGTGGCAACGACTTCTGCTACGCACGCTTTCCCGAGCGCAGCTTGTCAGATTCAAAGCATGTTGGGCATTAAAGGTTGCCCGGCATTTGACGTTGCAGCAGCCTGCGCAGGTTTCACCTATGCATTAAGCGTAGCCGATCAATACg                                                                                      <  2:329208/251‑1 (MQ=255)
                                                                                                                                                   acaTTGCCGCGCCAAACGAAACCGTTTCAACCATGGGCTTTGAAGCGGCGACACGCGCAATTGAGATGGCGGGCATTGAGAAAGACTAGATTGGCCTGATCGTTGTGGCAACGACTTCTGCTACGCACGCTTTCCCGAGCGCAGCTTGTCAGATTCAAAGCATGTTGGGCATTAAAGGTTGCCCGGCATTTGACGTTGCAGCAGCCTGCGCAGGTTTCACCTATGCATTAAGCGTAGCCGATCAATACGTg                                                                                    >  1:88549/1‑251 (MQ=255)
                                                                                                                                                                     aaaCCGTTTCAACCATGGGCTTTGAAGCGGCGACACGCGCAATTGAGATGGCGGGCATTGAGAAAGACTAGATTGGCCTGATCGTTGTGGCAACGACTTCTGCTACGCACGCTTTCCCGAGCGCAGCTTGTCAGATTCAAAGCATGTTGGGCATGAAAGTTTGCCCGGCATTTGACGTTGCAGCAGCCTGCGCATGTTTCTCCTATGCATTAAGCGTAGCCGATCAATACGTGAAATCTGGGGCGGTGAAg                                                                  <  2:64246/251‑1 (MQ=255)
                                                                                                                                                                                aCCATGGGCTTTGAAGCGGCGACACGCGCAATTGAGATGGCGGGCATTGAGAAAGACTAGATTGGCCTGATCGTTGTGGCAACGACTTCTGCTACGCTCGCTTTCCCGAGCGCAGCTTGTCAGATTCAAAGCATGTTGGGCATTAAAGGTTGCCCGGCATTTGACGTTGCAGCAGCCTGCGCAGGTTTCACCTATGCATTAAGCGTAGCCGATCAATACGTGAAAACTGGGGCGGTGAAGTATGCTCTGgt                                                       <  2:136588/251‑1 (MQ=255)
                                                                                                                                                                                 ccatggtgtttGATGCGGCGATACTCGCAATTGAGATGGCGAGCATTGAGAAAGACTAGATTTGCCTGATCGTTGTGGCAACGACTTCTGCTACGCACGCTTTCCCGAGCGCAGCTTGTCATATTCAAAGCATGTTGGGCATTAAAGGTTGCCCGGCATTTGACGTTGCAGCAGCCTGCGCAGGTTTCACCTATGCATTAAGCGTAGTCGATCAATACGTGAAATCTGGTGCGTTGAAGTAGTCTCTGgtc                                                      <  2:179898/243‑1 (MQ=255)
                                                                                                                                                                                         tttGAAGCGGCGACACGCGCAATTGAGATGGCGGGCATTGAGAAAGACTAGATTGGCCTGATCGTTGTGGCAACGACTTCTGCTACGCACGCTTTCCCGAGCGCAGCTTGTCAGATTCAAAGCATGTTGGGCATTAAAGGTTGCCCGGCATTTGACGTTGCAGCAGCCTGCGCAGGTTTCACCTATGCATTAAGCGTAGCCGATCAATACGTGAAATCTGGGGCGGTGAAGTATGCTCTGGTCGTCGGTTc                                              >  1:261131/1‑251 (MQ=255)
                                                                                                                                                                                                           gcAATTGAGATGGCGGGCATTGAGAAAGACTAGATTGGCCTGATCGTTGTGGCAACGACTTCTGCTACGCACGCTTTCCCGAGCGCAGCTTGTCAGATTCAAAGCATGTTGGGCATTAAAGGTTGCCCGGCATTTGACGTTGCAGCAGCCTGCGCAGGTTTCACCTATGCATTAAGCGTAGCCGATCAATACGTGAAATCTGGGGCGGTGAAGTATGCTCTGGTCGTCGGTTCCGATGTACTGGCGCGCAc                            >  1:326430/1‑251 (MQ=255)
                                                                                                                                                                                                            cAATTGAGATGGCGGGCATTGAGAAAGCCTCGCTTGGACTGATCGTTGAGGCAACGACTTCTGCTACGCACGCTTTCCCGAGCGCAGCTTGTCAGATTCAAAGCATGTTGGGCATTAAAGGTTGCCCGGCATTTGACGTTGCAGCAGCCTGCGCAGGTTTCACCTATGCATTAAGCGTAGCCGATCAATACGTGAAATCTGGGGCGGTGAAGTATGCTCTGGGCGTCGGTTCCGATGTACTGGCGCGCAcc                           >  2:312621/1‑251 (MQ=255)
                                                                                                                                                                                                                                   aaaGACTAGATTGGCCTGATCGTTGTGGCAACGACTTCTGCTACGCACTCTTTCTCGAGCGCAGCTTGTCAGATTCAAAGCATGTTGGGCATTAAAGGTTGCCCGGCATTTGACGTTGCAGCAGCCTGCGCAGGTTTCACCTATGCATTAAGCGTAGCCGATCAATACGTGAAATCTGGGGCGGTGAAGTATGCTCTGGTCGTCGGTTCCGATGTACTGGCGCGCACCTGCGATCCAACCGATCGTGGGAc    <  1:296489/251‑1 (MQ=255)
                                                                                                                                                                                                                                     aGACGAGATTGGTCTGATCGTTGTGGCAACGACGTCTGCTACGCACGCTTTCCCGAGCGCAGCTTGTCATGTTCAAAGCATGTTGGGCATTAAAGGTTGCCCGGCATTTGACGTTGCAGCAGCCTGCGCAGGTTTCACCTATGCATTAAGCGTAGCCGATCAATACGTGAAATCTGGGGCGGTGAAGTATGCTCTGGTCGTCGGTTGCGATGTACTGGCGCGCACCTGCGATCCAACTGATCGTGGGACta  <  2:153162/251‑1 (MQ=255)
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                       
CGAAAAGTGACTGAGCGTACATGTATACGAAGATTATTGGTACTGGCAGCTATCTGCCCGAACAAGTGCGGACAAACGCCGATTTGGAAAAAATGGTGGACACCTCTGACGAGTGGATTGTCACTCGTACCGGTATCCGCGAACGCCACATTGCCGCGCCAAACGAAACCGTTTCAACCATGGGCTTTGAAGCGGCGACACGCGCAATTGAGATGGCGGGCATTGAGAAAGACCAGATTGGCCTGATCGTTGTGGCAACGACTTCTGCTACGCACGCTTTCCCGAGCGCAGCTTGTCAGATTCAAAGCATGTTGGGCATTAAAGGTTGCCCGGCATTTGACGTTGCAGCAGCCTGCGCAGGTTTCACCTATGCATTAAGCGTAGCCGATCAATACGTGAAATCTGGGGCGGTGAAGTATGCTCTGGTCGTCGGTTCCGATGTACTGGCGCGCACCTGCGATCCAACCGATCGTGGGACTA  >  NC_000913/1148739‑1149218

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: