Predicted mutation
evidence seq id position mutation annotation gene description
RA CP006881 2,185,991 C→T K370K (AAG→AAA U712_10955 ← 3‑phosphoshikimate 1‑carboxyvinyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068812,185,9910CT98.2% 161.3 / ‑5.7 54K370K (AAG→AAAU712_109553‑phosphoshikimate 1‑carboxyvinyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (31/22);  total (31/23)
Fisher's exact test for biased strand distribution p-value = 4.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATAAGAAACGTGAATGGCATCTGTCTGCTCAATTTCAATTGGCTCTTCAGTTATACACGAAGCGATACCAAGCATCATTCCGATTCGATGATCTCCATGGCTGGACACTGTTGCGCCGCCTCTCAGTGTTTGCTTGCCATAAACCTTCATTCCATCTGCTGTCGGTTCAATCTCAGCACCTAGCTTACGAAGCTCAGAAACCACGGTGTCAATGCGGTTTGTTTCTTTTACTTTAAGCTCTGCCGCGTCCTTAATAACGGTGG  >  CP006881/2185859‑2186121
                                                                                                                                    |                                                                                                                                  
aTAAGAAACGTGAATGGCATCCGTGTGCTCGATTTCAATCGGCTCCTCCGTTATACAGGAAGCAATACCAAGCATCATTCCGATTCGATGATCTCCGTGGCTGGACACTGCAGCGCCGCCTTTCAACGTTTGTTTGCCa                                                                                                                              <  2:316471/139‑1 (MQ=255)
   atataCGTGAATGGCATCCGTGTGCTCGATTTCAATCGGCTCCTCCGTTATACAGGAAGCAATACCAAGCATCATTCCGATTCGATGATCTCCGTGGCTGGACACTGCAGCGCCTCCTTTCAACGTTTGTTTGACATa                                                                                                                            <  3:343508/134‑1 (MQ=255)
   aGAAACGTGAATGGCATCCGTGTGCTCGATTTCAATCGGCTCCTCCGTTATACAGGAAGCAATACCAAGCATCATTCCGATTCGATGATCTCCGTGGCTGGACACTGCAGCGCCGCCTTTCAACGTTTGTTTGCCATa                                                                                                                            <  5:128842/138‑1 (MQ=255)
   aGAAACGTGAATGGCATCCGTGTGCTCGATTTCAATCGGCTCCTCCGTTATACAGGAAGCAATACCAAGCATCATTCCGATTCGATGATCTCCGTGGCTGGACACTGCAGCGCCGCCTTTCAACGTTTGTTTGCCATa                                                                                                                            <  7:49515/138‑1 (MQ=255)
   aGAAACGTGAATCGCATCCGTGTGCTCGATTTCAATCGGCTCCTCCGTTATACAGGAAGCAATACCAAGCATCATTCCGATTCGATGATCTCCGTGGCTGGACACTGCAGCGCCGCCTTTCAACGTTTGTTTGCCATa                                                                                                                            <  4:429507/138‑1 (MQ=255)
    gAAACGTGAATGGCATCCGTGTGCTCGATTTCAATCGGCTCCTCCGTTATACAGGAAGCAATACCAAGCATCATTCCGATTCGATGATCTCCGTGGCTGGACACTGCAGCGCCGCCTGTCAACGTTTGTTTGCCATaaa                                                                                                                          >  4:94415/1‑139 (MQ=255)
       aCGTGAATGGCATCCGTGTGCTCGATTTCAATCGGCTCCTCCGTTATACAGGAAGCAATACCAAGCATCATTCCGATTCGATGATCTCCGTGGCTGGACACTGCAGCGCCGCCTTTCAACGTTTGTTTGCCATAAACCt                                                                                                                       <  4:281509/139‑1 (MQ=255)
       aCGTGAATGGCATCCGTGTGCTCGATTTCAATCGGCTCCTCCGTTATACAGGAAGCAATACCAAGCATCATTCCGATTCGATGATCTCCGTGGCTGGACACTGCAGCGCCGCCTTTCAACGTTTGTTTGCCATAAACCt                                                                                                                       <  4:287422/139‑1 (MQ=255)
                       tgtgCTCGATTTCAATCGGCTCCTCCGTTATACAGGAAGCAATACCAAGCATCATTCCGATTCGATGATCTCCGTGGCTGGACACTGCAGCGCCGCCTTTCAACGTTTGTTTGCCATAAACCTTCATTCCATCTGCTGt                                                                                                       <  6:349229/137‑1 (MQ=255)
                           cTCGATTTCAATCGGCTCCTCCGTTATACAGGAAGCAATACCAAGCATCATTCCGATTCGATGATCTCCGTGGCTGGACACTGCAGCGCCGCCTTTCAACGTTTGTTTGCCATAAACCTTCATTCCATCTGCTGTCGGt                                                                                                   >  5:314198/1‑139 (MQ=255)
                                   caatcGGCTCCTCCGTTATACAGGAAGCAATACCAAGCATCATTCCGATTCGATGATCTCCGTGGCTGGACACTGCAGCGCCGCCTTTCAACGTTTGTTTGCCATAAACCTTCATTCCATCTGCTGTCGGTTCAAttt                                                                                            >  7:319080/1‑136 (MQ=255)
                                   caatcGGCTCCTCCGTTATACAGGAAGCAATACCAAGCATCATTCCGATTCGATGATCTCCGTGGCTGGACACTGCAGCGCCGCCTTTCAACGTTTGTTTGCCATAAACCTTCATTCCATCTGCTGTCGGTTCAAttt                                                                                            >  4:34871/1‑136 (MQ=255)
                                        ggCTCCTCCGTTATACAGGAAGCAATACCAAGCATCATTCCGATTCGATGATCTCCGTGGCTGGACACTGCAGCGCCGCCTTTCAACGTTTGTTTGCCATAAACCTTCATTCCATCTGCTGTCGGTTCAATTTCAGCAc                                                                                      <  2:153181/139‑1 (MQ=255)
                                        ggCTCCTCCGTTATAACGGAAGCAATACCAAGCATCATTCCGATTCGATGATCTCCGTGGCTGGACACTGCAGCGCCGCCTTTCAACGTTTGTTTGCCATAAACCTTCATTCCATCTGCTGTCGGTTCAATTTCAGCAc                                                                                      <  8:392068/139‑1 (MQ=255)
                                                  ttATACAGGAAGCAATACCAAGCATCATTCCGATTCGATGATCTCCGTGGCTGGACACTGCAGCGCCGCCTTTCAACGTTTGTTTGCCATAAACCTTCATTCCATCTGCTGTCGGTTCAATTTCAGCAccc                                                                                    <  1:370813/131‑2 (MQ=255)
                                                  ttATACAGGAAGCAATACCAAGCATCATTCCGATTCGATGATCTCCGTGGCTGGACACTGCAGCGCCGCCTTTCAACGTTTGTTTGCCATAAACCTTCATTCCATCTGCTGTCGGTTCAATTTCAGCAccc                                                                                    >  2:370813/1‑130 (MQ=255)
                                                          gAAGCAATACCAAGCATCATTCCGATTCGATGATCTCCGTGGCTGGACACTGCAGCGCCGCCTTTCAACGTTTGTTTGCCATAAACCTTc                                                                                                                     <  2:330341/90‑1 (MQ=255)
                                                          gAAGCAATACCAAGCATCATTCCGATTCGATGATCTCCGTGGCTGGACACTGCAGCGCCGCCTTTCAACGTTTGTTTGCCATAAACCTTc                                                                                                                     >  1:330341/1‑90 (MQ=255)
                                                          gAAGCAATACCAAGCATCATTCCGATTCGATGATCTCCGTGGCTGGACACTGCAGCGCCGCCTTTCAACGTTTGTTTGCCATAAACCTTCATTCCATCTGCTGTCGGTTCAATTTCAGCACCCAGCTTGCGAAGCt                                                                       >  5:1098/1‑136 (MQ=255)
                                                          gAAGCAATACCAAGCATCATTCCGATTCGATGATCTCCGTGGCTGGACACTGCAGCGCCGCCTTTCAACGTTTGTTTGCCATAAACCTTCATTCCATCTGCTGTCGGTTCAATTTCAGCACCCAGCTTGCGAAGCt                                                                       >  6:332528/1‑136 (MQ=255)
                                                          gAAGCAATACCAAGCATCATTCCGATTCGATGATCTCCGTGGCTGGACACTGCAGCGCCGCCTTTCAACGTTTGTTTGCCATAAACCTTCATTCCATCTGCTGTCGGTTCAATTTCAGCACCCAGCTTGCGAAGCt                                                                       >  6:319526/1‑136 (MQ=255)
                                                          gAAGCAATACCAAGCATCATTCCGATTCGATGATCTCCGTGGCTGGACACTGAAGCGCCGCCTTTCAACGTTTGTTTGCCATAAACCTACATTCCATCTGCTGTCGGTTCAATTTCAGCACCCAGCTTGCGAAGca                                                                       >  7:239103/1‑135 (MQ=255)
                                                            aGCAATAGCAAGCATCATTCCGATTCGATGATCTCCGTGGCTGGACACTGCAGCGCCGCCTTTCAACGTTTGTTTGCCATAAACCTTCATTACATCTGCTGTGGGTTCAATTTCAGCACCCAGCATTCGAAGCTCAGa                                                                   >  4:69624/1‑138 (MQ=255)
                                                            aGCAATACCAAGCATCATTCCGATTCGATGATCTCCGTGGCTGGACACTGCAGCGCCGCCTTTCAACGTTTGTTTGCCATAAACCTTCATTCCATCTGCTGTCGGTTCAATTTCAGCACCCAGCTTGCGAAGCTCAGa                                                                   >  3:116379/1‑138 (MQ=255)
                                                            aGCAATACCAAGCATCATTCCGATTCGATGATCTCCGTGGCTGGACACTGCAGCGCCGCCTTTCAACGTTTGTTTGCCATAAACCTTCATTCCATCTGCTGTCGGTTCAATTTCAGCACCCAGCTTGCGAAGCTCAGa                                                                   >  4:131460/1‑138 (MQ=255)
                                                            aGCAATACCAAGCATCATTCCGATTCGATGATCTCCGTGGCTGGACACTGCAGCGCCGCCTTTCAACGTTTGTTTGCCATAAACCTTCATTCCATCTGCTGTCGGTTCAATTTCAGCACCCAGCTTGCGAAGCTCAGa                                                                   >  1:314712/1‑138 (MQ=255)
                                                            aGCAATACCAAGCATCATTCCGATTCGATGATCTCCGTGGCTGGACACTGCAGCGCCGCCTTTCAACGTTTGTTTGCCATAAACCTTCATTCCATCTGCTGTCGGTTCAATTTCAGCACCCAGCTTGCGAAGCTCAGa                                                                   >  6:185780/1‑138 (MQ=255)
                                                               aaTACCAAGCATCATTCCGATTCGATGATCTCCGTGGCTGGACACTGCAGCGCCGCCTTTCAACGTTTGTTTGCCATAAACCTTCATTCCATCTGCTGTCGGTTCAATTTCAGCACCCAGCTTGCGAAGCTCAGAAAc                                                                >  4:407960/2‑138 (MQ=255)
                                                               aaTACCAAGCATCATTCCGATTCGATGATCTCCGTGGCTGGACACTGCAGCGCCGCCTTTCAACGTTTGTTTGCCATAAACCTTCATTCCATCTGCTGTCGGTTCAATTTCAGCACCCAGCTTGCGAAGCTCAGAAAc                                                                >  7:55746/2‑138 (MQ=255)
                                                               aaTACCAAGCATCATTCCGATTCGATGATCTCCGTGGCTGGACACTGCAGCGCCGCCTTTCAACGTTTGTTTGCCATAAACCTTCATTCCATCTGCTGTCGGTTCAATTTCAGCACCCAGCTTGCGAAGCTCAGAAAc                                                                >  8:217664/2‑138 (MQ=255)
                                                                aTACCAAGCATCATTCCGATTCGATGATCTCCGTGGCTGGACACTGCAGCGCCGCCTTTCAACGTTTGTTTGCCATAAACCTTCATTCCATCTGCTGTCGGTTCAATTTCAGCACCCAGCTTGCGAAGCTCAGAAAcaa                                                              <  5:185780/139‑3 (MQ=255)
                                                                aTACCAAGCATCATTCCGATTCGATGATCTCCGTGGCTGGACACTGCAGCGCCGCCTTTCAACGTTTGTTTGCCATAAACCTTCATTCCATCTGCTGTCGGTTCAATTTCAGCACCCAGCTTGCGAAGCTCAGAAAcaa                                                              <  7:217664/139‑3 (MQ=255)
                                                                 tACCAAGCATCATTCCGATTCGATGATCTCCGTGGCTGGACACTGCAGCGCCGCCTTTCAACGTTTGTTTGCCATAAACCTTCATTCCATCTGCTGTCGGTTCAATTTCAGCACCCAGCTTGCGAAGCTCAGAAAcaa                                                              >  8:397866/1‑136 (MQ=255)
                                                                  aCCAAGCATCATTCCGATTCGATGATCTCCGTGGCTGGACACTGCAGCGCCGCCTTTCAACGTTTGTTTGCCATAAACCTTCATTCCATCTGCTGTCGGTTCAATTTCAGCACCCAGCTTGCGAAGCTCAGAAAcaa                                                              >  8:233294/1‑135 (MQ=255)
                                                                      aGCATCATTCCGATTCGATGATCTCCGTGGCTGGACACTGCAGCGCCGCCTTTCAACGTTTGTTTGCCATAAACCTTCATTCCATCTGCTGTCGGTTCAATTTCAGCACCCAGCTTGCGAAGCTCAGAAAcaa                                                              >  1:12034/1‑131 (MQ=255)
                                                                      aGCATCATTCCGATTCGATGATCTCCGTGGCTGGACACTGCAGCGCCGCCTTTCAACGTTTGTTTGCCATAAACCTTCATTCCATCTGCTGTCGGTTCAATTTCAGCACCCAGCTTGCGAAGCTCAGAAAcaa                                                              <  2:12034/133‑3 (MQ=255)
                                                                      aGCATCATTCCGATTCGATGATCTCCGTGGCTGGACACTGCAGCGCCGCCTTTCAACGTTTGTTTGCCATAAACCTTCATTCCATCTGCTGTCGGTTCAATTTCAGCACCCAGCTTGCGAAGCTCAGAAACAACAgtat                                                        >  2:257344/1‑137 (MQ=255)
                                                                           catTCCGATTCGATGATCTCCGTGGCTGGACACTGCAGCGCCGCCTTTCAACGTTTGTTTGCCATAAACCTTCATTCCATCTGCTGTCGGTTCAATTTCAGCACCCAGCTTGCGAAGCTCAGAAACAACAGTATCAAta                                                   >  8:67996/1‑138 (MQ=255)
                                                                             ttCCGATTCGATGATCTCCGTGGCTGGACACTGCAGCGCCGCCTTTCAACGTTTGTTTGCCATAAACCTTCATTCCATCTGCTGTCGGTTCAATTTCAGCACCCAGCTTGCGAAGCTCAGAAACAACAGTATCAAt                                                    <  8:32396/136‑1 (MQ=255)
                                                                             ttCCGATTCGATGATCTCCGTGGCTGGACACTGCAGCGCCGCCTTTCAACGTTTGTTTGCCATAAACCTTCATTCCATCTGCTGTCGGTTCAATTTCAGCACCCAGCTTGCGAAGCTCAGAAACAACAGTATCAAt                                                    >  7:32396/1‑136 (MQ=255)
                                                                             ttCCGATACGATGATCTCCGTGGCTGGACACTGCAGCGCCGCCTTTCAACGTATGTTTGCCATAAACCTTCATTCCATCTGCAGTCGGTTCAATTTAAGCACCCAGCAAGAGAAGCTCATAAACAACAGAATCAAtaac                                                 >  8:195839/1‑136 (MQ=255)
                                                                                 gATTCGATGATCTCCGTGGCTGGACACTGCAGCGCCGCCTTTCAACGTTTGTTTGCCATAAACCTTCATTCCATCTGCTGTCGGTTCAATTTCAGCACCCAGCTTGCGAAGCTCAGAAACAACAGTATCAATACGgttt                                             <  7:233294/139‑1 (MQ=255)
                                                                                        tgatCTCCGTGGCTGGACACTGCAGCGCCGCCTTTCAACGTTTGTTTGCCATAAACCTTCATTCCATCTGCTGTCGGTTCAATTTCAGCACCCAGCTTGCGAAGCTCAGAAACAACAGTATCAATACGGTTTGTTTCtt                                      <  5:332528/139‑1 (MQ=255)
                                                                                           tctcCGTGGCTGGACACTGCAGCGCCGCCGTTCAACGTTTGTTTGCCATAAACCTTCATTCCATCTGCTGTCGGTTCTATTTAAGCCCCCCGCTTGCGAAGCTCAGAAACCACAGTATCAAAACGGTTTGTTTCtttca                                   >  2:196123/1‑137 (MQ=255)
                                                                                            ctcCGTGGCTGGACACTGCAGCGCCGCCTTTCAACGTTTGTTTGCCATAAACCTTCATTCCATCTGCTGTCGGTTCAATTTCAGCACCCAGCTTGCGAAGCTCAGAAACAACAGTATCAATACGGTTTGTTTCTTTCAc                                  >  2:13220/1‑139 (MQ=255)
                                                                                               cgtGGCTGGACACTGCAGCGCCGCCTTTCAACGTTTGTTTGCCATAAACCTTCATTCCATCTGCTGTCGGTTCAATTTCATCACCCAGCTTGCGAATCTCAGAAACAACAGTAACAAGAAGGGTTGTTTATTTCAAttt                               >  8:279771/3‑139 (MQ=255)
                                                                                               cgtGGCTGGACACTGCAGCGCCGCCTTTCAACGTTTGTTTGCCATAAACCTTCATTCCATCTGCTGTCGGTTCAATTTCAGCACCCAGCTTGCGAAGCTCAGAAACAACAGTATCAATACGGTTTGTTTCtttc                                    <  5:319526/132‑2 (MQ=255)
                                                                                                       gACACTGCAGCGCCGCCTTTCAACGTTTGTTTGCCATAAACCTTCATTCCATCTGCTGTCGGTTCAATTTCAGCACCCAGCTTGCGAAGCTCAGAAACAACAGTATCAATACGGTTTGTTTCTTTCACTTTTAGCTCTg                       >  8:57636/1‑139 (MQ=255)
                                                                                                         cacTGCAGCGCCGCCGTTCAACGTTTGCTTGCCATAAACCTTCATTCCATCTGCTGTCGGTTCAATTTCAGCACCCAGCTTGCGAAGCTCAGAAACAACAGTATCAATACGGTTTGTTTCTTTCACTTTTAGCTCTGcc                     <  1:432495/139‑1 (MQ=255)
                                                                                                            tgcagcgcCGCCTTTCAACGTTTGTTTGCCATAAACCTTCATTCCATCTGCTGTCGGTTCAATTTCAGCACCCAGCTTGCGAAGCTCAGAAACAACAGTATCAATACGGTTTGTTTCTTTCACTTTTAGCTCTGCcgcg                  >  7:404003/5‑139 (MQ=255)
                                                                                                            tgcagcgcCGCCTTTCAACGTTTGTTTGCCATAAACCTTCATTCCATCTGCTGTCGGTTCAATTTCAGCACCCAGCTTGCGAAGCTCAGAAACAACAGTATCAATACGGTTTGTTTCTTTCACTTTTAGCTCTGCcgcg                  >  4:144052/5‑139 (MQ=255)
                                                                                                                    cgccTTTCAACGTTTGTTTGCCATAAACCTTCATTCCATCTGCTGTCGGTTCAATTTCAGCACCCAGCTTGCGAAGCTCAGAAACAACAGTATCAATACGGTTTGTTTCTTTCACTTTTAGCTCTGCCGCGTCCTtaat          <  7:397866/139‑1 (MQ=255)
                                                                                                                    cgccTTTCAACGTTTGTTTGCCATAAACCTTCATTCCATCTGCTGTCGGTTCAATTTCAGCACCCAGCTTGCGAAGCTCAGAAACAACAGTATCAATACGGTTTGTTTCTTTCACTTTTAGCTCTGCCGCGTCCTtaat          <  3:34871/139‑1 (MQ=255)
                                                                                                                    cgccTTTCAACGTTTGTTTGCCATAAACCTTCATTCCATCTGCTGTCGGTTCAATTTCAGCACCCAGCTTGCGAAGCTCAGAAACAACAGTATCAATACGGTTTGTTTCTTTCACTTTTAGCTCTGCCGCGTCCTtaat          <  1:13220/139‑1 (MQ=255)
                                                                                                                          tcaacgtttgtttGCCATAAACCTTCATTCCATCTGCTGTCGGTTCAATTTCAGCACCCAGCTTGCGAAGCTCAGAAACAACAGTATCAATACGGTTTGTTTCTTTCACTTTTAGCTCTGCCGCGTCCTTAATAACg      <  4:201102/132‑1 (MQ=255)
                                                                                                                           caacgtttgtttGCCATAAACCTTCATTCCATCTGCTGTCGGTTCAATTTCAGCACCCAGCTTGCGAAGCTCAGAAACAACAGTATCAATACGGTTTGTTTCTTTCACTTTTAGCTCTGcc                     >  7:406116/5‑121 (MQ=255)
                                                                                                                           caacgtttgtttGCCATAAACCTTCATTCCATCTGCTGTCGGTTCAATTTCAGCACCCAGCTTGCGAAGCTCAGAAACAACAGTATCAATACGGTTTGTTTCTTTCACTTTTAGCTCTGcc                     <  8:406116/117‑1 (MQ=255)
                                                                                                                            aacgtttgtttGCCATAAACCTTCATTCCATCTGCTGTCGGTTCAATTTCAGCACCCAGCTTGCGAAGCTCAGAAACAACAGTATCAATACGGTTTGTTTCTTTCACTTTTAGCTCTGCCGCGTCCTTAATAACggtgg  <  7:279771/136‑1 (MQ=255)
                                                                                                                            aacgtttgtttGCCATAAACCTTCATTCCATCTGCTGTCGGGTCAATTTCAGCACCCAGCTTGCGAAGCTCAGAAACAACAGTATCAATACGGTTTGTTTCTTTCACTTTTAGCTCTGCCGCGTCCTTAATAACggtgg  <  3:245276/136‑1 (MQ=255)
                                                                                                                                 ttgtttgcCATAAACCTTCATTCCATCTGCTGTCGGTTCAATTTCAGCACCCAGCTTGCGAAGCTCAGAAACAACAGTATCAATACGGTTTGTTTCTTTCACTTTTAGCTCTGCcgc                   <  2:104993/117‑1 (MQ=255)
                                                                                                                                 ttgtttgcCATAAACCTTCATTCCATCTGCTGTCGGTTCAATTTCAGCACCCAGCTTGCGAAGCTCAGAAACAACAGTATCAATACGGTTTGTTTCTTTCACTTTTAGCTCTGCcgc                   >  1:104993/1‑117 (MQ=255)
                                                                                                                                    |                                                                                                                                  
ATAAGAAACGTGAATGGCATCTGTCTGCTCAATTTCAATTGGCTCTTCAGTTATACACGAAGCGATACCAAGCATCATTCCGATTCGATGATCTCCATGGCTGGACACTGTTGCGCCGCCTCTCAGTGTTTGCTTGCCATAAACCTTCATTCCATCTGCTGTCGGTTCAATCTCAGCACCTAGCTTACGAAGCTCAGAAACCACGGTGTCAATGCGGTTTGTTTCTTTTACTTTAAGCTCTGCCGCGTCCTTAATAACGGTGG  >  CP006881/2185859‑2186121

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: