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
                                                                                                                                 |                                                                                                                                   
gggAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCGCCGAt                                                                                                                           >  1:349191/1‑139 (MQ=255)
  gAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCGCCGATAc                                                                                                                         <  1:220768/139‑1 (MQ=255)
                            ccagccaCTATTATTCAACGGATTATTGACGGTAGATATTTATCGGTAAGGAGGAAATGCAGAGCTTTTCGTCACGTTAAATGGTAAATGCCACACCCGGTCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCa                                                                                               <  1:74315/139‑1 (MQ=255)
                             cagctaTAATTATTCAACAGATTAAGGAGGGTGAAAATTTTACGGGAAGGAGAAAATGCACAGCTATTCGTCACGTTAAATGTTCAATGCCAAACCCTAGCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCAt                                                                                              <  1:576277/139‑1 (MQ=255)
                                              cAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCa                                                                             >  1:427966/1‑139 (MQ=255)
                                                 aTTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGa                                                                          >  1:120605/1‑139 (MQ=255)
                                                      tGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCag                                                                     >  1:85376/1‑139 (MQ=255)
                                                               atatTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCa                                                                             >  2:508410/1‑122 (MQ=255)
                                                               atatTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCa                                                                             <  1:508410/122‑1 (MQ=255)
                                                                     aTCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGccccg                                                                                                                                 >  2:211613/1‑64 (MQ=255)
                                                                     aTCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCg                                                      <  2:255087/139‑1 (MQ=255)
                                                                     aTCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTAAATGGTCAATGCCACACCCTGGccccg                                                                                                                                 <  1:211613/65‑2 (MQ=255)
                                                                            aggaggAAATGCAGAGCTATTCGTCACTTTCAATGGTCAATGCCACACCCTGGCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCaga                                                                    <  1:106623/118‑1 (MQ=255)
                                                                            aggaggAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCaga                                                                    >  2:106623/1‑118 (MQ=255)
                                                                                ggaAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAg                                           <  2:85376/139‑1 (MQ=255)
                                                                                ggaAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAg                                           <  2:119599/139‑1 (MQ=255)
                                                                                   aaTGCAGAGCTATTCGTCACGTTCAATGGTCAAAGCCACACCCTGGCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAgcgc                                        >  1:553777/1‑139 (MQ=255)
                                                                                    aTGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGcc                                       >  2:258406/1‑139 (MQ=255)
                                                                                              aTTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGa                             <  1:258406/139‑1 (MQ=255)
                                                                                               ttCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGa                             <  1:72532/138‑1 (MQ=255)
                                                                                               ttCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGa                             >  2:72532/1‑138 (MQ=255)
                                                                                               ttCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGAc                            >  1:59412/1‑139 (MQ=255)
                                                                                               ttCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGAc                            <  2:59412/139‑1 (MQ=255)
                                                                                                          cAATGGTCAATGCCACACCCTGGCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATc                 >  1:422147/1‑139 (MQ=255)
                                                                                                             tgttaaaTGCCCCACCCTGACACCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATACGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATcgcg              <  1:226258/134‑1 (MQ=255)
                                                                                                             tGGTCAATGCCACACCCTGGCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATcgcg              <  2:376397/139‑1 (MQ=255)
                                                                                                                         cacCCTGGCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattgaca  >  1:440936/1‑139 (MQ=255)
                                                                                                                          acCCTGGCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAg                                           >  1:595472/1‑97 (MQ=255)
                                                                                                                          acCCTGGCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAg                                           <  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: