Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,959,515 T→A 21.6% intergenic (+70/+17) ilvC → / ← ppiC ketol‑acid reductoisomerase (NADP(+))/peptidyl‑prolyl cis‑trans isomerase C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,959,5150TA21.6% 56.8 / 12.9 32intergenic (+70/+17)ilvC/ppiCketol‑acid reductoisomerase (NADP(+))/peptidyl‑prolyl cis‑trans isomerase C
Reads supporting (aligned to +/- strand):  ref base T (15/10);  new base A (3/4);  total (18/14)
Fisher's exact test for biased strand distribution p-value = 6.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.62e-01

AGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCGGAA  >  NC_000913/3959384‑3959657
                                                                                                                                   |                                                                                                                                              
aggcaggtAAGCGGCTGCGCGGCTATCTGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCt                                                                                                                                     <  2:551138/143‑1 (MQ=38)
 ggtaggtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCGCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCtt                                                                                                                                    <  1:594004/143‑1 (MQ=255)
  gtaggtaAGAAACTGCGGGGGTATAGGACAGATAGGAAACGTATTGCTGTTGCGGGTTAAGAGCGCGCTGATGCCCTCACCCCGACCCTCTCAAAAAAAAAAAAAAAAAAAACACTCAAGGCCTTCTCCAGGAGAAGGCCTTg                                                                                                                                   <  2:753637/143‑1 (MQ=255)
    aggtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGGAAACGCTCAAGGCCTTCTCCTGGAGAAGGCGTTGcg                                                                                                                                 >  1:189879/1‑142 (MQ=255)
    aggtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCAGGAGAAGGCCTTgag                                                                                                                                 <  1:165853/143‑3 (MQ=255)
        aaGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACGCAAGGCCTTCTCCAGGAGAAGGCCTTgagtgtt                                                                                                                             <  1:468387/143‑7 (MQ=255)
             aCTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGc                                                                                                                        <  2:570089/143‑1 (MQ=255)
              gtGCGCGGCGATAGGACAGAGATGAGACGTATTGCGGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCGCCCACAGGGAGAGGGAGAAAACACTCGAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCg                                                                                                                       <  2:641382/142‑1 (MQ=37)
                                            aTTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGCAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGaa                                                                                         >  1:678309/1‑143 (MQ=255)
                                                 tGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGCAGGCCTCGCTATTCGTTGCGGTACAGCCCCTCAACGATGGGATATCCGACCTgga                                                                                    >  1:693694/1‑142 (MQ=37)
                                                               tGCGCGGTGATGCCCTCACCCCGACAATCTTCGGCCGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCg                                                                      <  1:211496/143‑1 (MQ=255)
                                                                             ctctCCCCGGTCCTCTCCCACAGGGAGCGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATCAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTAc                                                        <  2:673905/143‑1 (MQ=255)
                                                                                  cccGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGAc                                                   >  1:786017/1‑143 (MQ=255)
                                                                                        ccTCTCCCTGTGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAgaga                                             >  1:160155/1‑143 (MQ=255)
                                                                                              ccACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGAACCTTAATGATGTGATATCCGAACTGGGTCTTCTGCCGGCCGGTCGGCTCCAGTCCCGGACAGACGAAAACCa                                       >  1:225829/1‑143 (MQ=255)
                                                                                              ccACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCAATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGCGCGCGAAGCGGGCTGCTCGGCTCCAGCACCGGACAAGCCAACACCa                                       >  2:163503/1‑143 (MQ=255)
                                                                                              ccACAGGGAGAGGGAGAAAACACTCAAGGCCATCTACAAGAAGATGCCTTGCTATTAGTTGCGGTACAGCCCCATAAGCCTGTGATATCCGAACTTGGTGTGCAGCGGGCCGGTCTGCGCCAGTACCTGACAAGAGAAGACCa                                       >  2:597841/1‑143 (MQ=255)
                                                                                                    gcagaggcacaAAACACTCACCGCCGTCTCCTCGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAAGGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTAt                                 <  2:574268/141‑1 (MQ=255)
                                                                                                      agagggagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCg                               >  1:45889/1‑143 (MQ=255)
                                                                                                        agggagaAAACGCTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTTATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGaa                             <  1:10138/143‑1 (MQ=255)
                                                                                                           gagaAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGc                          >  1:348024/1‑143 (MQ=255)
                                                                                                               aaaCACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCAGGCCGGTCGGCTCCAGTAGCGCACAAGAGAAAACCACTTCATCGAAAGCCGGa                      >  2:363735/1‑143 (MQ=255)
                                                                                                                taatagcAAGGCCTTCTCCAGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGcc                         >  1:453498/7‑139 (MQ=255)
                                                                                                                taatagcAAGGCCTTCTCCAGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGcc                         <  2:453498/133‑1 (MQ=255)
                                                                                                                   acTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGCACCa                  >  1:24020/1‑143 (MQ=255)
                                                                                                                        cggCCTCCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGTACCATCTGa             <  1:234353/142‑1 (MQ=255)
                                                                                                                          gCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACTCCGGAAACaccaggcg           >  1:606923/1‑136 (MQ=255)
                                                                                                                          gCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGAcc           >  1:285788/1‑143 (MQ=255)
                                                                                                                          gCCTTCTCCAGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTGGGCGCCAGTACCGGCCAAGAGAAAACGACGTTATCGAGCGCCGGCACGAGCTGAcc           >  2:185201/1‑143 (MQ=255)
                                                                                                                            cTTATCCTGGAGGAGGCCTTGCTATTAGTTAAAAAAAAGCACCTTAATGATGTGATATACGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTAACGGACAAGAGACAACAACTTCATGGAACGCAGGAACCATCGAACCCt         >  2:185696/1‑143 (MQ=255)
                                                                                                                              tctcCTCGAGAAGGCCTTGCTGTTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTgg       <  2:360292/143‑1 (MQ=255)
                                                                                                                                  cTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCGGa   >  1:642183/1‑143 (MQ=255)
                                                                                                                                   tGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCGGaa  <  1:508124/143‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                              
AGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACTTTATCGAACGCCGGAACCATCTGACCCTGGCGGAA  >  NC_000913/3959384‑3959657

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