Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 9,146 T→G 46.5% A303A (GCT→GCG talB → transaldolase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092739,1460TG46.5% ‑0.9 / 40.4 28A303A (GCT→GCGtalBtransaldolase B
Reads supporting (aligned to +/- strand):  ref base T (7/8);  new base G (5/8);  total (12/16)
Fisher's exact test for biased strand distribution p-value = 7.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.81e-01

AGCGAAGGGGCTATCGAACGTAAACTGTCTTACACCGGCGAAGTGAAAGCGCGTCCGGCGCGTATCACTGAGTCCGAGTTCCTGTGGCAGCACAACCAGGATCCAATGGCAGTAGATAAACTGGCGGAAGGTATCCGTAAGTTTGCTATTGACCAGGAAAAACTGGAAAAAATGATCGGCGATCTGCTGTAATCATTCTTAGCGTGACCGGGAAGTCGGTCACGCTACCTCTTCTGAAGCCTGTCTGTCACTCCCTTCGCAGTGTATCATTCTGTTTAACGAGACTGTTTAAAC  >  CP009273/9000‑9293
                                                                                                                                                  |                                                                                                                                                   
agCGAAGGCGCCGTTGAGCGTAAGCTGGTATATACCGGTGAAGTGAAAGCGCGTCCGGAGCGTATTACCGAATCTGAGTTCCTGTGGCAGCATAACCAGGATCCGATGGCGGTAGACAAACTGGCGGAAGGTATCCGTAAGTTTGCGAtc                                                                                                                                                  <  1:134147/150‑2 (MQ=255)
   gAAGGGGCTATCGAACGTAAACTGTCTTACACCGGCGAAGTGAAAGCGCGTCCGGCGCGTATCACTGAGTCCGAGTTCCTGTGGCAGCACAACCAGGATCCAATGGCAGTAGATAAACTGGCGGAAGGTATCCGTAAGTTTGCTATTGAc                                                                                                                                               >  2:123044/1‑150 (MQ=255)
        ggCTATCGAACGTAAACTGTCTTACACCGGCGAAGTGAAAGCGCGTCCGGCGCGTATCACTGAGTCCGAGTTCCTGTGGCAGCACAACCAGGATCCAATGGCAGTAGATAAACTGGCGGAAGGTATCCGTAAGTTTGCTATTGACCAGGa                                                                                                                                          <  1:72400/150‑1 (MQ=255)
         gCTATCGAACGTAAACTGTCTTACACCGGCGAAGTGAAAGCGCGTCCGGCGCGTATCACTGAGTCCGAGTTCCTGTGGCAGCACAACCAGGATCCAATGGCAGTAGATAAACTGGCGGAAGGTATCCGTAAGTTTGCTATTGACCAGGaa                                                                                                                                         <  1:125336/150‑1 (MQ=255)
             tCGAACGTAAACTGTCTTACACCGGCGAAGTGAAAGCGCGTCCGGCGCGTATCACTGAGTCCGAGTTCCTGTGGCAGCACAACCAGGATCCAATGGCAGTAGATAAACTGGCGGAAGGTATCCGTAAGTTTGCTATTGACCAGGaaaaa                                                                                                                                      <  1:145340/149‑1 (MQ=255)
                           gtatatacCGGTGAAGTGAAAGCGCGTCCGGAGCGTATTACCGAATCTGAGTTCCTGTGGCAGCATAACCAGGATCCGATGGCGGTAGACAAACTGGCGGAAGGTATCCGTAAGTTTGCGATCGACCAGGAAAAACTGGAAAAAATGATc                                                                                                                       >  1:10501/7‑150 (MQ=255)
                            tatatacCGGTGAAGTGAAAGCGCGTCCGGAGCGTATTACCGAATCTGAGTTCCTGTGGCAGCATAACCAGGATCCGATGGCGGTAGACAAACTGGCGGAAGGTATCCGTAAGTTTGCGATCGACCAGGAAAAACTGGAAAAAATGATCg                                                                                                                      <  1:138557/145‑1 (MQ=255)
                              tatacCGGTGAAGTGAAAGCGCGTCCGGAGCGTATTACCGAATCTGAGTTCCTGTGGCAGCATAACCAGGATCCGATGGCGGTAGACAAACTGGCGGAAGGTATCCGTAAGTTTGCGATCGACCAGGAAAAACTGGAAAAAATGATCgac                                                                                                                    <  1:190215/147‑3 (MQ=255)
                                tacCGGTGAAGTGAAAGCGCGTCCGGAGCGTATTACCGAATCTGAGTTCCTGTGGCAGCATAACCAGGATCCGATGGCGGTAGACAAACTGGCGGAAGGTATCCGTAAGTTTGCGATCGACCAGGAAAAACTGGAAAAAATGATCGACg                                                                                                                   >  2:152619/2‑149 (MQ=255)
                                   cGGCGAAGTGAAAGCGCGTCCGGCGCGTATCACTGAGTCCGAGTTCCTGTGGCAGCACAACCAGGATCCAATGGCAGTAGATAAACTGGCGGAAGGTATCCGTAAGTTTGCTATTGACCAGGAAAAACTGGAAAAAATGATCGGCGATct                                                                                                               <  1:159742/150‑1 (MQ=255)
                                    ggtgAAGTGAAAGCGCGTCCGGAGCGTATTACCGAATCTGAGTTCCTGTGGCAGCATAACCAGGATCCGATGGCGGTAGACAAACTGGCGGAAGGTATCCGTAAGTTTGCGATCGACCAGGAAAAACTGGAAAAAATGATCGACGAActg                                                                                                              >  1:71862/4‑150 (MQ=255)
                                    ggtgAAGTGAAAGCGCGTCCGGAGCGTATTACCGAATCTGAGTTCCTGTGGCAGCATAACCAGGATCCGATGGCGGTAGACAAACTGGCGGAAGGTATCCGTAAGTTTGCGATCGACCAGGAAAAACTGGAAAAAATGATCGACGAActg                                                                                                              <  2:183261/147‑1 (MQ=255)
                                      tgAAGTGAAAGCGCGTCCGGAGCGTATTACCGAATCTGAGTTCCTGTGGCAGCATAACCAGGATCCGATGGCGGTAGACAAACTGGCGGAAGGTATCCGTAAGTTTGCGATCGACCAGGAAAAACTGGAAAAAATGATCGACGAActgct                                                                                                            <  2:92404/149‑1 (MQ=255)
                                      cGAAGTGAAAGCGCGTCCGGCGCGTATCACTGAGTCCGAGTTCCTGTGGCAGCACAACCAGGATCCAATGGCAGTAGATAAACTGGCGGAAGGTATCCGTAAGTTTGCTATTGACCAGGAAAAACTGGAAAAAATGATCGGCGATctgct                                                                                                            >  2:160600/1‑150 (MQ=255)
                                       gAAGTGAAAGCGCGTCCGGAGCGTATTACCGAATCTGAGTTCCTGTGGCAGCATAACCAGGATCCGATGGCGGTAGACAAACTGGCGGAAGGTATCCGTAAGTTTGCGATCGACCAGGAAAAACTGGAAAAAATGATCGACGAActgctg                                                                                                           >  1:123210/1‑150 (MQ=255)
                                        aaGTGAAAGCGCGTCCGGAGCGTATTACCGAATCTGAGTTCCTGTGGCAGCATAACCAGGATCCGATGGCGGTAGACAAACTGGCGGAAGGTATCCGTAAGTTTGCGATCGACCAGGAAAAACTGGAAAAAATGATCGACGAACTGCTGt                                                                                                          <  2:82091/150‑1 (MQ=255)
                                                gcgcGTCCGGAGCGTATTACCGAATCTGAGTTCCTGTGGCAGCATAACCAGGATCCGATGGCGGTAGACAAACTGGCGGAAGGTATCCGTAAGTTTGCGATCGACCAGGAAAAACTGGAAAAAATGATCGACGAACTGCTGTaagagaat                                                                                                  <  1:63609/150‑7 (MQ=255)
                                                    gTCCGGAGCGTATTACCGAATCTGAGTTCCTGTGGCAGCATAACCAGGATCCGATGGCGGTAGACAAACTGGCGGAAGGTATCCGTAAGTTTGCGATCGACCAGGAAAAACTGGAAAAAATGATCGACGAActgct                                                                                                            >  2:64725/1‑136 (MQ=255)
                                                     tCCGGAGCGTATTACCGAATCTGAGTTCCTGTGGCAGCATAACCAGGATCCGATGGCGGTAGACAAACTGGCGGAAGGTATCCGTAAGTTTGCGATCGACCAGGAAAAACTGGAAAAAATGATCGACGAActgctg                                                                                                           <  1:64725/136‑1 (MQ=255)
                                                                  aCTGAGTCCGAGTTCCTGTGGCAGCACAACCAGGATCCAATGGCAGTAGATAAACTGGCGGAAGGTATCCGTAAGTTTGCTATTGACCAGGAAAAACTGGAAAAAATGATCGGCGATCTGCTGTAATCATTCTTAGCGTGACCGGGAAGt                                                                                >  1:191113/1‑150 (MQ=255)
                                                                          cGAGTTCCTGTGGCAGCACAACCAGGATCCAATGGCAGTAGATAAACTGGCGGAAGGTATCCGTAAGTTTGCTATTGACCAGGAAAAACTGGAAAAAATGATCGGCGATCTGCTGTAATCATTCTTAGCGTGACCGGGAAGTCGGTCACg                                                                        <  1:190167/150‑1 (MQ=255)
                                                                           gAGTTCCTGTGGCAGCACAACCAGGATCCAATGGCAGTAGATAAACTGGCGGAAGGTATCCGTGAGTTTGCTATTGACCAGGAAAAACTGGAAAAAATGATCGGCGATCTGCTGTAATCATTCTTAGCGTGACCGGGAAGTCGGTCACGc                                                                       <  2:4973/150‑1 (MQ=255)
                                                                                               ccAGGATCCAATGGCAGTAGATAAACTGGCGGAAGGTATCCGTAAGATTGCTATTGACCAGGAAAAACTGGAAAAAATGATCGGCGATCTGCTGTAATCATTCTTAGCGTGACCGGGAAGTCGGTCACGCTACCTCTTCTGAAGCctgtc                                                   >  2:61712/1‑150 (MQ=255)
                                                                                                                  gATAAACTGGCGGAAGGTATCCGTAAGTTTGCTATTGACCAGGAAAAACTGGAAAAAATGATCGGCGATCTGCTGTAATCATTCTTAGCGTGACCGGGAAGTCGGTCACGCTACCTCTTCTGAAGCCTGTCTGTCACTCCCTTCGCAgtg                                >  1:182/1‑150 (MQ=255)
                                                                                                                   aTAAACTGGCGGAAGGTATCCGTAAGTTTGCTATTGACCAGGAAAAACTGGAAAAAATGATCGGCGATCTGCTGTAATCATTCTTAGCGTGACCGGGAAGTCGGTCACGCTACCTCTTCTGAAGCCTGTCTGTCACTCCCTTCGCAgtg                                >  2:1489/1‑149 (MQ=255)
                                                                                                                        cTGGCGGAAGGTATCCGTAAGTTTGCTATTGACCAGGAAAAACTGGAAAAAATGATCGGCGATCTGCTGTAATCATTCTTAGCGTGACCGGGAAGTCGGTCACGCTACCTCTTCTGAAGCCTGTCTGTCACTCCCTTCGCAGTGTATCAt                          <  2:182/150‑1 (MQ=255)
                                                                                                                                       cGTAAGTTTGCTATTGACCAGGAAAAACTGGAAAAAATGATCGGCGATCTGCTGTAATCATTCTTAGCGTGACCGGGAAGTCGGTCACGCTACCTCTTCTGAAGCCTGTCTGTCACTCCCTTCGCAGTGTATCATTCTGTTTAACGAGAc           <  1:123044/150‑1 (MQ=255)
                                                                                                                                                gCTATTGACCAGGAAAAACTGGAAAAAATGATCGGCGATCTGCTGTAATCATTCTTAGCGTGACCGGGAAGTCGGTCACGCTACCTCTTCTGAAGCCTGTCTGTCACTCCCTTTGCAGTGTATCATTCTGTTTAACGAGACTGTTTAAAc  >  1:39163/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
AGCGAAGGGGCTATCGAACGTAAACTGTCTTACACCGGCGAAGTGAAAGCGCGTCCGGCGCGTATCACTGAGTCCGAGTTCCTGTGGCAGCACAACCAGGATCCAATGGCAGTAGATAAACTGGCGGAAGGTATCCGTAAGTTTGCTATTGACCAGGAAAAACTGGAAAAAATGATCGGCGATCTGCTGTAATCATTCTTAGCGTGACCGGGAAGTCGGTCACGCTACCTCTTCTGAAGCCTGTCTGTCACTCCCTTCGCAGTGTATCATTCTGTTTAACGAGACTGTTTAAAC  >  CP009273/9000‑9293

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: