Predicted mutation
evidence seq id position mutation annotation gene description
RA CP006881 2,185,850 A→G F417F (TTT→TTC U712_10955 ← 3‑phosphoshikimate 1‑carboxyvinyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068812,185,8500AG94.6% 105.0 / ‑2.3 37F417F (TTT→TTCU712_109553‑phosphoshikimate 1‑carboxyvinyltransferase
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (12/23);  total (13/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

AGCGCATCCCCCCTTTTCGTTTTTTCTTTTCTCTCTATAAGACAAGCTATGTTTTTTTGCCTCGGATTATGATTCTCACAGGCATTCTGCACAAAAAATCCTGAAGGAAAATCTTCAGGATTTGTTTGAAAGCTTATTTAAATGCTCGAAAAAGGTTGGATAAGAAACGTGAATGGCATCTGTCTGCTCAATTTCAATTGGCTCTTCAGTTATACACGAAGCGATACCAAGCATCATTCCGATTCGATGATCTCCATGGCTGGACACTGTTGCGCCGCCTCTC  >  CP006881/2185700‑2185982
                                                                                                                                                      |                                                                                                                                    
aGCACATCCCCCCTTTTCGTTTTTCTATCCGTTCTATAAGACAAGCTATGTTTTTTCGCCTCGAATTATGA‑‑‑‑‑‑‑‑‑‑‑‑‑TCTGCATAAAAAATCCTGAAGTAAAACTTCAGGATTTTTTCGAAAGCTTATTTAAATGCTCGaagaag                                                                                                                                   >  1:135544/1‑139 (MQ=255)
aGCACATCCCCCCTTTTCGTTTTTCTATCCGTTCTATAAGACAAGCTATGTTTTTTCGCCTCGAATTATGA‑‑‑‑‑‑‑‑‑‑‑‑‑TCTGCATAAAAAATCCTGAAGTAAAACTTCAGGATTTTTTCGAAAGCTTATTTAAATGCTCGaagaag                                                                                                                                   >  4:207313/1‑139 (MQ=255)
  cacATCCCCCCTTTTCGTTTTTCTATCCGTTCTATAAGACAAGCTATGTTTTTTCGCCTCGAATTATGA‑‑‑‑‑‑‑‑‑‑‑‑‑TCTGCATAAAAAATCCTGAAGTAAAACTTCAGGATTTTTTCGAAAGCTTATTTAAATGCTCGAAGAAGGt                                                                                                                                 <  5:73831/137‑1 (MQ=255)
  cacATCCCCCCTTTTCGTTTTTCTATCCGTTCTATAAGACAAGCTATGTTTTTTCGCCTCGAATTATGA‑‑‑‑‑‑‑‑‑‑‑‑‑TCTGCATAAAAAATCCTGAAGTAAAACTTCAGGATTTTTTCGAAAGCTTATTTAAATGCTCGAAGAAGGt                                                                                                                                 <  4:408029/137‑1 (MQ=255)
       ccccccTTTTCGTTTTTCTATCCGTTCTATAAGACAAGCTATGTTTTTTCGCCTCGAATTATGA‑‑‑‑‑‑‑‑‑‑‑‑‑TCTGCATAAAAAATCCTGAAGTAAAACTTCAGGATTTTTTCGAAAGCTTATTTAAATGCTCGAAGAAGGTTGGAt                                                                                                                            >  8:80215/1‑139 (MQ=255)
                                             gCTATGTTTTTTTGCCTCGGATTATGATTCTCACAGGCATTCTGCACAAAAAATCCTGAAGGAAAATCTTCAGGATTTGTTTGAAAGCTTATTTAAATGCTCGAAAAAGGt                                                                                                                                 <  6:35321/111‑1 (MQ=255)
                                             gCTATGTTTTTTTGCCTCGGATTATGATTCTCACAGGCATTCTGCACAAAAAATCCTGAAGGAAAATCTTCAGGATTTGTTTGAAAGCTTATTTAAATGCTCGAAAAAGGt                                                                                                                                 >  5:35321/1‑111 (MQ=255)
                                                                                 atgatCTGCATAAAAAATCCTGAAGTAAAACTTCAGGATTTTTTCGAAAGCTTATTTAAATGCTCGAAGAAGGTTGGATAAGAAACGTGAATGGCATCCGTGTGCTCGATTTCAATCGGCTCCTCCGTTATACAGGAAg                                                                >  3:281509/5‑139 (MQ=255)
                                                                                 atgatCTGCATAAAAAATCCTGAAGTAAAACTTCAGGATTTTTTCGAAAGCTTATTTAAATGCTCGAAGAAGGTTGGATAAGAAACGTGAATGGCATCCGTGTGCTCGATTTCAATCGGCTCCTCCGTTATACAGGAAg                                                                >  3:287422/5‑139 (MQ=255)
                                                                                  tgatCTGCATAAAAAATCCTGAAGTAAAACTTCAGGATTTTTTCGAAAGCTTATTTAAATGCTCGAAGAAGGTTGGATAAGAAACGTGAATGGCATCCGTGTGCTCGATTTCAATCGGCTCCTCCGTTATACAGGAAg                                                                <  8:250184/135‑1 (MQ=255)
                                                                                  tgatCTGCATAAAAAATCCTGAAGTAAAACTTCAGGATTTTTTCGAAAGCTTATTTAAATGCTCGAAGAAGGTTGGATAAGAAACGTGAATGGCATCCGTGTGCTCGATTTCAATCGGCTCCTCCGTTATACAGGAAg                                                                <  7:210510/135‑1 (MQ=255)
                                                                                  tgatCTGCATAAAAAATCCTGAAGTAAAACTTCAGGATTTTTTCGAAAGCTTATTTAAATGCTCGAAGAAGGTTGGATAAGAAACGTGAATGGCATCCGTGTGCTCGATTTCAATCGGCTCCTCCGTTATACAGGAAg                                                                <  6:228438/135‑1 (MQ=255)
                                                                                  tgatCTGCATAAAAAATCCTGAAGTAAAACTTCAGGATTTTTTCGAAAGCTTATTTAAATGCTCGAAGAAGGTTGGATAAGAAACGTGAATGGCATCCGTGTGCTCGATTTCAATCGGCTCCTCCGTTATACAGGAAg                                                                <  3:315626/135‑1 (MQ=255)
                                                                                  tgatCTGCATAAAAAATCCTGAAGTAAAACTTCAGGATTTTTTCGAAAGCTTATTTAAATGCTCGAAGAAGGTTGGATAAGAAACGTGAATGGCATCCGTGTGCTCGATTTCAATCGGCTCCTCCGTTATACAGGAAg                                                                <  4:157767/135‑1 (MQ=255)
                                                                                       tGCATAAAAAATCCTGAAGTAAAACTTCAGGATTTTTTCGAAAGCTTATTTAAATGCTCGAAGAAGGTTGGATAAGAAACGTGAATGGCATCCGTGTGCTCGATTTCAATCGGctcct                                                                               >  4:291376/1‑116 (MQ=255)
                                                                                       tGCATAAAAAATCCTGAAGTAAAACTTCAGGATTTTTTCGAAAGCTTATTTAAATGCTCGAAGAAGGTTGGATAAGAAACGTGAATGGCATCCGTGTGCTCGATTTCAATCGGctcct                                                                               <  3:291376/118‑3 (MQ=255)
                                                                                       tGCATAAAAAATCCTGAAGTAAAACTTCAGGATTTTTTCGAAAGCTTATTTAAATGCTCGAAGAAGGTGGGATAAGAAACGTGAATGGCATCCGTGTGCTCGATTTCAATCGGCTCCTCCGTTATACAGGAAGCAATAc                                                          >  3:82397/1‑139 (MQ=255)
                                                                                        gCATAAAAAATCCTGAAGTAAAACTTCAGGATTTTTTCGAAAGCTTATTTAAATGCTCGAAGAAGGTTGGATAAGAAACGTGAATGGCATCCGTGTGCTCGATTTCAATCGGCTCCTCCGtt                                                                          <  4:291729/122‑1 (MQ=255)
                                                                                        gCATAAAAAATCCTGAAGTAAAACTTCAGGATTTTTTCGAAAGCTTATTTAAATGCTCGAAGAAGGTTGGATAAGAAACGTGAATGGCATCCGTGTGCTCGATTTCAATCGGCTCCTCCGtt                                                                          >  3:291729/1‑122 (MQ=255)
                                                                                         cataaaaatTTCTGAATTAAAACTTCAGGATTTTTTCGAAAGCTTATTTAAATGCTCTAAGAAGGTTGGATATGAATCGTGAATGGCATGCGTGTGCTCGATTTCAATCGGCTGCTCCGTTATACAGGAAGCAATAcc                                                         <  1:153485/135‑1 (MQ=255)
                                                                                         cataaaaaaTCCTGAAGTAAAACTTCAGGATTTTTTCGAAAGCTTATTTAAATGCTCGAAGAAGGTTGGATAAGAAACGTGAATGGCATCCGTGTGCTCGATTTCAATCGGCTCCTCCGTTATACAGGAAGCAATAcc                                                         >  5:37440/4‑138 (MQ=255)
                                                                                         cataaaaaaTCCTGAAGTAAAACTTCAGGATTTTTTCGAAAGCTTATTTAAATGCTCGAAGAAGGTTGGATAAGAAACGTGAATGGCATCCGTGTGCTCGATTTCAATCGGCTCCTCCGTTATACAGGAAGCAATAcc                                                         <  7:123030/135‑1 (MQ=255)
                                                                                         cataaaaaaTCCTGAAGTAAAACTTCAGGATTTTTTCGAAAGCTTATTTAAATGCTCGAAGAAGGTTGGATAAGAAACGTGAATGGCATCCGTGTGCTCGATTTCAATCGGCTCCTCCGTTATACAGGAAGCAATAcc                                                         >  1:153181/4‑138 (MQ=255)
                                                                                         cataaaaaaTCCTGAAGTAAAACTTCAGGATTTTTTCGAAAGCTTATTTAAATGCTCGAAGAAGGTTGGATAAGAAACGTGAATGGCATCCGTGTGCTCGATTTCAATCGGCTCCTCCGTTATACAGGAAGCAATAcc                                                         <  3:348155/135‑1 (MQ=255)
                                                                                         cataaaaaaTCCTGAAGTAAAACTTCAGGATTTTTTCGAAAGCTTATTTAAATGCTCGAAGAAGGTTGGATAAGAAACGTGAATGGCATCAGTGTGCTCGATTTCAATCGGCTCCTCCGTTATACAGGAAGCAATACCa                                                        <  3:248194/136‑1 (MQ=255)
                                                                                         cataaaaaaTCCTGAAGTAAATCTTCAGGATTTTTTCGTTAGCTTATTTAAATGCTCGAAGAAGGTTGGATAAGAAACGTGAATGTCATCCGTGTGCTCGATTTCAATCGGCTCCTCCGTTATACAGGAAGCAATAcc                                                         <  2:67735/135‑1 (MQ=255)
                                                                                                   ccTGAAGTAAAACTTCAGGATTTTTTCGAAAGCTTATTTAAATGCTCGAAGAAGGTTGGATAAGAAACGTGAATGGCATCCGTGTGCt                                                                                                 <  2:329713/88‑1 (MQ=255)
                                                                                                   ccTGAAGTAAAACTTCAGGATTTTTTCGAAAGCTTATTTAAATGCTCGAAGAAGGTTGGATAAGAAACGGGAATGGCATCCGTGTGCt                                                                                                 >  1:329713/1‑88 (MQ=255)
                                                                                                        aagtaaaaCTTCAGGATTTTTTCGAAAGCTTATTTAAATGCTCGAAGAAGGTTGGATAAGAAACGTGAATGGCATCCGTGTGCTCGATTTCAATCGGCTCCTCCGTTATACAGGAAGCAATACCAAGCATCATTCCg                                            <  3:207313/133‑1 (MQ=255)
                                                                                                         agtaaaaCTTCAGGATTTTTTCGAAAGCTTATTTAAATGCTCGAAGAAGGTTGGATAAGAAACGTGAATGGCATCCGTGTGCTCGATTTCAATCGGCTCCTCCGTTATACAGGAAGCAATACCAAGCATCATTCCGAtt                                         <  4:220739/136‑1 (MQ=255)
                                                                                                             aaacTTCAGGATTTTTTCGAAAGCTTATTTAAATGCTCGAAGAAGGTTGGATAAGAAACGTGAATGGCATCCGTGTGCTCGATTTCAATCGGCTCCTCCGTTATACAGGAAGCAATACCAAGCATCATTCCGATTCgat                                     >  6:128842/4‑139 (MQ=255)
                                                                                                                  tCAGGATTTTTTCGAAAGCTTATTTAAATGCTCGAAGAAGGTTGGATAAGAAACGTGAATGGCATCCGAGTGCTCGATTTCAATCGGCTCCTCCGTTATACAGGAAGCAATACCAAGCATCATTCCGATTCGATGAtct                                >  7:99958/1‑139 (MQ=255)
                                                                                                                    aGGATTTTTTCGAAAGCTTATTTAAATGCTCGAAGAAGGTTGGATAAGAAACGTGAATGGCATCCGTGTGCTCGATTTCAATCGGCTCCTCCGTTATACAGGAAGCAATACCAAGCATCATTCCGATTCGATGAtctc                               <  6:37440/138‑1 (MQ=255)
                                                                                                                     ggATTTTTTCGAAAGCTTATTTAAATGCTCGAAGAAGGTTGGATAAGAAACGTGAATGGCATCCGTGTGCTCGATTTCAATCGGCTCCTCCGTTATACAGGAAGCAATACCAAGCATCATTCCGATTCGATGATCTccg                             <  2:235710/139‑2 (MQ=255)
                                                                                                                              cgAAAGCTTATTTAAATGCTCGAAGAAGGTTGGATAAGAAACGTGAATGGCATCCGTGTGCTCGATTTCAATCGGCTCCTCCGTTATACAGGAAGCAATACCAAGCATCATTCCGATTCGATGATCTCCGTGGCTGGa                     <  7:199228/137‑1 (MQ=255)
                                                                                                                              cgAAAGCTTATTTAAATGCTCGAAGAAGGTTGGATAAGAAACGTGAATGGCATCCGTGTGCTCGATTTCAATCGGCTCCTCCGTTATACAGGAAGCAATACCAAGCATCATTCCGATTCGATGATCTCCGTGGCTGGa                     <  7:428395/137‑1 (MQ=255)
                                                                                                                              cgAAAGCTTATTTAAATGCTCGAAGAAGGTTGGATAAGAAACGTGAATGGCATCCGTGTGCTCGATTTCAATCGGCTCCTCCGTTATACAGGAAGCAATACCAAGCATCATTCCGATTCGATGATCTCCGTGGCTGGa                     <  8:370249/137‑1 (MQ=255)
                                                                                                                                        tttAAATGCTCGAAGAAGGTTGGATAAGAAAAGTGAATGGCATCCGTGTGCTCGATTTCAATCGGCTCCTCCGTTATACAGGAAGAAATAACAAGCATCATTCCGATTCGATTATCTCCGTGGCTGGACACTgcagctc          <  3:120465/139‑7 (MQ=255)
                                                                                                                                                 tCGAAGAAGGTTGGATAAGAAACGTGAATGGCATCCGTGTGCTCGATTTCAATCGGCTCCTCCGTTATACAGGAAGCAATACCAAGCATCATTCCGATTCGATGATCTCCGTGGCTGGACACTGCAGCGCCGCCtttc  >  2:14393/1‑138 (MQ=255)
                                                                                                                                                      |                                                                                                                                    
AGCGCATCCCCCCTTTTCGTTTTTTCTTTTCTCTCTATAAGACAAGCTATGTTTTTTTGCCTCGGATTATGATTCTCACAGGCATTCTGCACAAAAAATCCTGAAGGAAAATCTTCAGGATTTGTTTGAAAGCTTATTTAAATGCTCGAAAAAGGTTGGATAAGAAACGTGAATGGCATCTGTCTGCTCAATTTCAATTGGCTCTTCAGTTATACACGAAGCGATACCAAGCATCATTCCGATTCGATGATCTCCATGGCTGGACACTGTTGCGCCGCCTCTC  >  CP006881/2185700‑2185982

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