Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,519,036 A→T 7.5% V10D (GTT→GAT)  aroK ← shikimate kinase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,519,0360AT 7.5% 145.1 / 3.7 57V10D (GTT→GAT) aroKshikimate kinase 1
Reads supporting (aligned to +/- strand):  ref base A (26/24);  new base T (2/2);  total (30/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.09e-01

AATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACTATTCGTTAATGATAAACCCGCTTCGCTCAGAGAGCGCCGCAGGACATGAACTGAAACCTCATAAGATATTGCGAGAGTCAGACTGAAAATTAT  >  NC_000913/3518908‑3519178
                                                                                                                                |                                                                                                                                              
aaTCTAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTg                                                                                                                                  <  1:675172/143‑1 (MQ=255)
    gAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGtt                                                                                                                              >  1:879450/1‑143 (MQ=255)
    gAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGtt                                                                                                                              >  2:710079/1‑143 (MQ=255)
     aaaaCCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGatataacgagt                                                                                                                             >  1:554276/1‑136 (MQ=255)
     aaaaCCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGttt                                                                                                                             >  1:610429/1‑143 (MQ=255)
        aCCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTTGCTCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTtctc                                                                                                                          <  1:416820/143‑1 (MQ=255)
          ccAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTg                                                                                                                        >  1:439918/1‑143 (MQ=255)
                   aTCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCAtttttt                                                                                                               <  1:710079/143‑1 (MQ=255)
                       gCTCCGGTTCGTTTCTCAATCTCTTGATCGGAGTGGCGCGCAAACAAATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCGTGTTCCGGCACGCCTAGGACCACCCAGAAAGATATGGCGTTTCTCTGCCAGTTTTTCGGt                                                                                                           >  2:1081209/1‑143 (MQ=255)
                           cGGTTCGTTTCTCATTCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTGTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGtacta                                                                                                       <  2:14669/143‑1 (MQ=255)
                            ggTTCGTTTCTTAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGACGCATAATAGTGCTTTTTCCGGAACCCCTAGGCCCAACATGAAAGATACTGAGTTTCTCTGTGATTTTTTCggacgtaa                                                                                                      >  2:123784/1‑137 (MQ=255)
                             gTTCGTTTCTCAATCTCTGGAGCGGAGTCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAAGAGTGCTTTTTCCGGCACCCATGGGCCCAACCAGAAAGAGATTGCGTTTCTCTGCCATTTTTTCGGTactact                                                                                                     <  2:572942/143‑1 (MQ=255)
                                 gTTTCTTGAGCTCGTGATCGGACTCGTAAAGTTCCATATTGAGTTGTTGAGCTAACTGGCGGCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCGCTGCCATTTTTTCGGTACTACTAAGa                                                                                                 <  2:308833/143‑1 (MQ=255)
                                  tttCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCGCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTCCTACTGAGAc                                                                                                >  2:1107353/1‑143 (MQ=255)
                                     ctcAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACTAt                                                                                             >  1:688/1‑143 (MQ=255)
                                        aaTCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGTAAGATTTTGCGGCTCTCTTCCATTTTTTCGGTACTACTAAGACTATTCg                                                                                          >  2:229576/1‑143 (MQ=255)
                                          tctcTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACTATTCGtt                                                                                        <  1:109203/143‑1 (MQ=255)
                                          tctcTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACTATTCGtt                                                                                        <  1:935934/143‑1 (MQ=255)
                                              ttGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACTATTCGTTAATg                                                                                    <  1:422652/143‑1 (MQ=255)
                                                  tCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCATCCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACTATTCGTTAATGATaa                                                                                >  2:752959/1‑143 (MQ=255)
                                                  tCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCATCCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACTATTCGTTAATGATaa                                                                                >  2:199352/1‑143 (MQ=255)
                                                   cGGAATCGTATGATTCGATATTGAGTTGTTGAGCTGACTGGCGCCCCATAGTGCTTTTTCCGGCGCCCAGAGGCCCAACCAGAAAGATAGTGCGTTTCTCTGCCATTTTTTCGCTACTACTACGACTATTCGTTAATGATaaa                                                                               <  2:81377/143‑1 (MQ=255)
                                                   cGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACTATTCGTTAATGATaaa                                                                               >  2:347295/1‑143 (MQ=255)
                                                     gAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACTATTCGTTAATGATAAAcc                                                                             >  1:1007894/1‑143 (MQ=255)
                                                        tCGTAAATTTCCATATTGCGCTGTTGGGCTAACTGCCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATGGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACTATTCGTTAATGATAAACCCGc                                                                          <  2:935386/143‑1 (MQ=255)
                                                          gTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACTATTCGATAATGATAAACCCGCtt                                                                        >  2:64747/1‑143 (MQ=255)
                                                          gTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATATCATAGGAAGTACTAAGACTATTCGTTAAGGAAAAACCCACtt                                                                        >  2:328291/1‑143 (MQ=255)
                                                           tAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCAATCAACACAGGCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTCAGACTATTAGTTAATGATAAACCCGCTTc                                                                       >  1:415628/1‑143 (MQ=255)
                                                                ttCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACTATTCGTTAATGATAAACCCGCTTCGCTCa                                                                  <  2:879450/143‑1 (MQ=255)
                                                                     tatTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACTATTCGTTAATGACAAACCCGCTTCGCTTagagag                                                             >  2:769096/1‑143 (MQ=255)
                                                                          cgttgttAAACTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATAGTGCGTTTCTCTGCCATTTTCTCGGTACTACTAAGACTATTCGTTAATGATAAACCCGCTTCGCTCAGAGAGcgccg                                                        <  2:237/142‑1 (MQ=255)
                                                                            ttgttgAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGTGTTTCTGTGCCATTTTTTCGGTACTACTAAGACTATTGGTTAATGATAAACGCGCTTCGCGCAGAGAgcgcgggt                                                      >  1:795923/1‑139 (MQ=255)
                                                                            ttgttgAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACTATTCGTTAATGATAAACCCGCTTCGCTCAGAGAGCGCCGCa                                                      >  1:1061694/1‑143 (MQ=255)
                                                                             tgttgAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACTATTCGTTAATGATAAACCCGCTTCGCTCAGAGAGCGCCGCAg                                                     <  1:488095/143‑1 (MQ=255)
                                                                                 gAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACTATTCGTTAATGATAAACCCGCTTCGCTCAGAGAGCGCCGCAGGACa                                                 >  2:142689/1‑143 (MQ=255)
                                                                                  aGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACTATTCGTTAATGATAAACCCGCTTCGCTCCGAGAGCGCCGCAGGACAt                                                >  2:1093051/1‑143 (MQ=255)
                                                                                  aGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACTATTCGTTAATGATAAACCCGCTTCGCTCAGAGAGCGCCGCAGGACAt                                                <  2:421584/143‑1 (MQ=255)
                                                                                   gCTAACTGGCGCCCAATAGTGCTCTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACTATTCGTTAATGATAAACCCGCTTCGCTCAGAGAGCGCCGCAGGACATg                                               <  1:384817/143‑1 (MQ=255)
                                                                                     tAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAACGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACTATTCGTTAATGATAAACCCGCTTCGCTCAGAGAGCGCCGCAGGACATGaa                                             <  1:3271/143‑1 (MQ=255)
                                                                                         tGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACTATTCGTTAATGATAAACCCGCTTCGCTCAGAGAGCGCACCAGGACATGacatga                                         >  2:821692/1‑138 (MQ=255)
                                                                                           gctcccAATAGTGCTTTTTCCGGCACCCATCGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTCCGGTACTACTAAGACTATTCGTTAATGATAAACCCGCTTCGCTCAGAGAGCGCCGCAGGACATGAACTGaaa                                       <  2:54890/140‑1 (MQ=255)
                                                                                            cgcCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACTATTCGTTAATGATAAACCCGCTTCGCTCAGAGAGCGCCGCAGGACATGAACTGAAAc                                      >  1:657525/1‑143 (MQ=255)
                                                                                             gcCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACTATTCGTTAATGATAAACCCGCTTCGCTCAGAGAGCGCCGCAGGACATGAACTGAAAcc                                     >  2:999537/1‑143 (MQ=255)
                                                                                               ccAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACTATTCGTTAATGATAAACCCGCTTCGCTCAGAGCGCGCCGCAGGACAGGAACTGAAACCTc                                   >  2:899261/1‑143 (MQ=255)
                                                                                                 aaTAGTGCTTTTTCCGTCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACGACTAAGCCTATTCGTTAATGATAAACCCGCTTCGCTCAGAGAGCGCCGCAGGACATGAACTGAAACCTCAt                                 <  2:410288/143‑1 (MQ=255)
                                                                                                  aTAGTGCTTTTTCCGGCACCCATAGGCCCAGCCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACTATTCGTTAATGATAAACCCGCTTCGCTCAGAGAGCGCCGCAGGACATGAACTGAAACCTCATa                                <  2:827777/143‑1 (MQ=255)
                                                                                                         tttttCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACTATTCGTTAATGATAAACCCGCTTCGCTCCGAGAGCGCCGCAGGACATGAACTGAAACCTCATAAGATAtt                         >  1:958833/1‑143 (MQ=255)
                                                                                                            ttCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACTATTCGTTAATGATAAACCCGCTTCGCTCAGAGAGCGCCGCAGGACATGAACTGAAACCTCATAAGATATTGCg                      >  2:633809/1‑143 (MQ=255)
                                                                                                            tgccGGCGTCCATCGGCCCAACCAGAAAGATATGGCGTTTCTCTGCCATTTGTTCGGCACTACTAAGACTATGCGTTAATGATAAACCCGCTTCGCTCAGAGAGCGCCGCAGGACATGAACTGAAACCTCATAAGATATTGCg                      <  2:610429/141‑1 (MQ=255)
                                                                                                               cGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACTATTCGTTAATGATAAACCCGCTTCGCTCAGAGAGCGCCGCAGGACATGAACTGAAACCTCATAAGATATTGCgaga                   <  2:97315/143‑1 (MQ=255)
                                                                                                                 gCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCGGCCATTTTTTCGGTACTACTAAGACTATTCGTTAATGATAAACCCGCTTCGCTCAGGGAGCGCCGCAGGACGTGAACGGAAACCGCATCAGATCTCGGGAGAGt                 >  2:863074/1‑143 (MQ=255)
                                                                                                                     ccATATGCCCAACCCGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACTATTCGTTAATGATAAACCCGCTTCGCTCAGAGAGCGCCGCAGGACATGAACTGAAACCTCATAAGATATTGCGAGAGTCAGa             <  1:164065/143‑1 (MQ=255)
                                                                                                                     ccATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACTATTCGTTAATGATAAACCCGCTTCGCTCAGAGAGCGCCGCAGGACATGAACTGAAACCTCATAAGATATTGCGAGAGTCAGa             >  1:870649/1‑143 (MQ=255)
                                                                                                                      cATCGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACTATTCGTTAATGATAAACCCGCTTCGCTCAGAGAGCGCCGCAGGACATGAACTGAAACCTCATAAGATATTGCGAGAGTCAGAc            <  1:188543/143‑1 (MQ=255)
                                                                                                                      cATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACTATTCGTTAATGATAAACCCGCTTCGCTCAGAGAGCGCCGCAGGACATGAACTGAAACCTCATAAGATATTGCGAGAGTCAGAc            <  1:484866/143‑1 (MQ=255)
                                                                                                                         aGGCCCATCCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACTATTCGTTAATGATAAACCCGCTTCGCTCAGAGAGCGCCGCAGGACATGAACTGAAACCTCATAAGATATTGCGAGAGTCAGACTGa         <  1:752959/143‑1 (MQ=255)
                                                                                                                         aGGCCCATCCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACTATTCGTTAATGATAAACCCGCTTCGCTCAGAGAGCGCCGCAGGACATGAACTGAAACCTCATAAGATATTGCGAGAGTCAGACTGa         <  1:199352/143‑1 (MQ=255)
                                                                                                                                aCCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACTATTCGTTAATGATAAACCCGCTTCGCTCAGAGAGCGCCGCAGGACATGAACTGAAACCTCATAAGATATTGCGAGAGTCAGACTGAAAATTAt  <  2:1007894/143‑1 (MQ=255)
                                                                                                                                |                                                                                                                                              
AATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACTATTCGTTAATGATAAACCCGCTTCGCTCAGAGAGCGCCGCAGGACATGAACTGAAACCTCATAAGATATTGCGAGAGTCAGACTGAAAATTAT  >  NC_000913/3518908‑3519178

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