Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,984,452 (C)5→4 coding (1141/1182 nt) yqeF ← short chain acyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,984,4480C.96.2% 101.2 / ‑2.5 26coding (1145/1182 nt)yqeFshort chain acyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base . (13/12);  minor base T (0/1);  total (13/13)
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

GGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACA  >  NC_000913/2984319‑2984579
                                                                                                                                 |                                                                                                                                   
gggAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCGCCGAt                                                                                                                           >  1:349191/1‑139 (MQ=255)
  gAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCGCCGATAc                                                                                                                         <  1:220768/139‑1 (MQ=255)
                            ccagccaCTATTATTCAACGGATTATTGACGGTAGATATTTATCGGTAAGGAGGAAATGCAGAGCTTTTCGTCACGTTAAATGGTAAATGCCACACCC‑GGTCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCa                                                                                               <  1:74315/139‑1 (MQ=255)
                             cagctaTAATTATTCAACAGATTAAGGAGGGTGAAAATTTTACGGGAAGGAGAAAATGCACAGCTATTCGTCACGTTAAATGTTCAATGCCAAACCCTAG‑CCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCAt                                                                                              <  1:576277/139‑1 (MQ=255)
                                              cAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCa                                                                             >  1:427966/1‑139 (MQ=255)
                                                 aTTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGa                                                                          >  1:120605/1‑139 (MQ=255)
                                                      tGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCag                                                                     >  1:85376/1‑139 (MQ=255)
                                                               atatTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCa                                                                             >  2:508410/1‑122 (MQ=255)
                                                               atatTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCa                                                                             <  1:508410/122‑1 (MQ=255)
                                                                     aTCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGccccg                                                                                                                                 >  2:211613/1‑64 (MQ=255)
                                                                     aTCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCg                                                      <  2:255087/139‑1 (MQ=255)
                                                                     aTCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTAAATGGTCAATGCCACACCCTGGccccg                                                                                                                                 <  1:211613/65‑2 (MQ=255)
                                                                            aggaggAAATGCAGAGCTATTCGTCACTTTCAATGGTCAATGCCACACCCTGG‑CCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCaga                                                                    <  1:106623/118‑1 (MQ=255)
                                                                            aggaggAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCaga                                                                    >  2:106623/1‑118 (MQ=255)
                                                                                ggaAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAg                                           <  2:85376/139‑1 (MQ=255)
                                                                                ggaAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAg                                           <  2:119599/139‑1 (MQ=255)
                                                                                   aaTGCAGAGCTATTCGTCACGTTCAATGGTCAAAGCCACACCCTGG‑CCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAgcgc                                        >  1:553777/1‑139 (MQ=255)
                                                                                    aTGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGcc                                       >  2:258406/1‑139 (MQ=255)
                                                                                              aTTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGa                             <  1:258406/139‑1 (MQ=255)
                                                                                               ttCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGa                             <  1:72532/138‑1 (MQ=255)
                                                                                               ttCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGa                             >  2:72532/1‑138 (MQ=255)
                                                                                               ttCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGAc                            >  1:59412/1‑139 (MQ=255)
                                                                                               ttCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGAc                            <  2:59412/139‑1 (MQ=255)
                                                                                                          cAATGGTCAATGCCACACCCTGG‑CCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATc                 >  1:422147/1‑139 (MQ=255)
                                                                                                             tgttaaaTGCCCCACCCT‑GACACCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATACGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATcgcg              <  1:226258/134‑1 (MQ=255)
                                                                                                             tGGTCAATGCCACACCCTGG‑CCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATcgcg              <  2:376397/139‑1 (MQ=255)
                                                                                                                         cacCCTGG‑CCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattgaca  >  1:440936/1‑139 (MQ=255)
                                                                                                                          acCCTGG‑CCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAg                                           >  1:595472/1‑97 (MQ=255)
                                                                                                                          acCCTGG‑CCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAg                                           <  2:595472/97‑1 (MQ=255)
                                                                                                                                 |                                                                                                                                   
GGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACA  >  NC_000913/2984319‑2984579

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