Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,151,873 T→A 14.3% intergenic (+22/‑66) acpP → / → fabF acyl carrier protein (ACP)/3‑oxoacyl‑[acyl‑carrier‑protein] synthase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,151,8730TA14.3% 76.2 / 6.1 28intergenic (+22/‑66)acpP/fabFacyl carrier protein (ACP)/3‑oxoacyl‑[acyl‑carrier‑protein] synthase II
Reads supporting (aligned to +/- strand):  ref base T (13/11);  new base A (2/2);  total (15/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 = 7.43e-01

TTGACACCGTTGAGCTGGTAATGGCTCTGGAAGAAGAGTTTGATACTGAGATTCCGGACGAAGAAGCTGAGAAAATCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGCGTCGTGTAGTTGTGACCGGACTGGGCATGTTGTCTCCTGTCGGCAATACCGTAGAGTCTACCTGGAAAGCTC  >  NC_000913/1151727‑1152020
                                                                                                                                                  |                                                                                                                                                   
ttGACACCGTTGAGCTGGTAATGGCTCTGGAAGAAGAGTTTGATACTGAGATTCCGGACGAAGAAGCTGAGAAAATCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTc                                                                                                                                                   >  2:486083/1‑149 (MQ=255)
               tGGTAATGGCTCTGGAAGAAGAGTTTGATACTGAGATTCCGGACGAAGAAGCTGAGAAAATCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTc                                                                                                                                                   >  1:253513/1‑134 (MQ=255)
               tGGTAATGGCTCTGGAAGAAGAGTTTGATACTGAGATTCCGGACGAAGAAGCTGAGAAAATCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTc                                                                                                                                                   <  2:253513/134‑1 (MQ=255)
                 gTAATGGCTCTGGAAGAAGAGTTTGATACTGAGATTCCGGACGAAGAAGCTGAGAAAATCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTAt                                                                                                                                  >  1:367822/1‑149 (MQ=255)
                            ggAAGAAGAGTTTGATACTGAGATTCCGGACGAAGAAGCTGAGAAAATCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTccc                                                                                                                       <  1:133907/149‑1 (MQ=255)
                              aagaagAGTTTGATACTGAGATTCCGGACGAAGAAGCTGAGAAAATCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCAc                                                                                                                     >  4:129556/1‑149 (MQ=255)
                                              tGAGATTCCGGACGAAGAAGCTGAGAAAATCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCAtt                                                                                                           >  4:416172/1‑143 (MQ=255)
                                              tGAGATTCCGGACGAAGAAGCTGAGAAAATCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCAtt                                                                                                           <  3:416172/143‑1 (MQ=255)
                                                    tCCGGACGAAGAAGCTGAGAAAATCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTccctccct                                                                                               <  4:620712/149‑1 (MQ=255)
                                                       ggACGAAGAAGCTGAGAAAATCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTgga                                                                                            <  1:486083/149‑1 (MQ=255)
                                                                  cTGAGAAAATCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACgtg                                                                                 <  3:129556/149‑1 (MQ=255)
                                                                                             ccATTGATTACATCAACGGCCACCAGGCGTAAGTGTACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCttt                                                                                                                              <  1:276117/77‑1 (MQ=255)
                                                                                             ccATTGATTACATCAACGGCCACCAGGCGTAAGTGTACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCttt                                                                                                                              >  2:276117/1‑77 (MQ=255)
                                                                                              cATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTGCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGCGTCGTGTAGTTGTGACCGGAc                                                     >  3:34730/1‑149 (MQ=255)
                                                                                              cATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGCGTCGTGTAGTTGTGACCGGAc                                                     >  1:297604/1‑149 (MQ=255)
                                                                                                    ttACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGCGTCGTGTAGTTGTGACCGGACTGGGCa                                               >  3:453910/1‑149 (MQ=255)
                                                                                                       cATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGCGTCGTGTAGTTGTGACCGGACTGGGCAtgt                                            >  4:460829/1‑149 (MQ=255)
                                                                                                        aTCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGCGTCGTGTAGTtgt                                                             <  1:20257/131‑1 (MQ=255)
                                                                                                        aTCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGCGTCGTGTAGTtgt                                                             >  2:20257/1‑131 (MQ=255)
                                                                                                                       gCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGCGTCGTGTAGTTGTGACCGGACTGGGCATGTTGTCTCCCGTCGGCaa                            <  2:297604/149‑1 (MQ=255)
                                                                                                                           aaGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGCGTCGTGTAGTTGTGACCGGACTGGGCATGTTGTCTCCTGTCGGCAATAcc                        <  2:367822/149‑1 (MQ=255)
                                                                                                                                           gCGGTCGAACGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGc                                                                          <  3:180774/83‑1 (MQ=255)
                                                                                                                                           gCGGTCGAACGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGc                                                                          >  4:180774/1‑83 (MQ=255)
                                                                                                                                             ggTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGCGTCGTGTAGTTGTGACCGGACTGGGCATGTTGTCTCCTGTCGGCAATACCGTAGAGTCTACCTGGaaa      <  3:460829/149‑1 (MQ=255)
                                                                                                                                             ggTCGAACGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGCGTCGTGTAGTTGTGACCGGACt                                                    >  3:101654/1‑103 (MQ=255)
                                                                                                                                             ggTCGAACGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGCGTCGTGTAGTTGTGACCGGACt                                                    <  4:101654/103‑1 (MQ=255)
                                                                                                                                                 gTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGCGTCGTGTAGTTGTGACCGGACTGGGCATGTTGTCTCCTGTCGGCAATACCGTAGAGTCTACCTGGAAAGctc  >  3:258372/1‑149 (MQ=255)
                                                                                                                                                 gTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGCGTCGTGTAGTTGTGACCGGACTGGGCATGTTGTCTCCTGTCGGCAATACCGTAGAGTCTACCTGGAAAGctc  >  2:117771/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TTGACACCGTTGAGCTGGTAATGGCTCTGGAAGAAGAGTTTGATACTGAGATTCCGGACGAAGAAGCTGAGAAAATCACCACCGTTCAGGCTGCCATTGATTACATCAACGGCCACCAGGCGTAAGTGAACATCTCCAGGCGGTCGTTCGACCGCCTGAGTTTTATCTTTTTGTCCCACTAGAATCATTTTTTCCCTCCCTGGAGGACAAACGTGTCTAAGCGTCGTGTAGTTGTGACCGGACTGGGCATGTTGTCTCCTGTCGGCAATACCGTAGAGTCTACCTGGAAAGCTC  >  NC_000913/1151727‑1152020

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: