Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,234,563 G→A 9.5% D147D (GAC→GAT)  malF ← maltose transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,234,5630GA 9.5% 112.5 / 2.2 42D147D (GAC→GAT) malFmaltose transporter subunit
Reads supporting (aligned to +/- strand):  ref base G (16/22);  new base A (2/2);  total (18/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 = 7.82e-01

TTGTTGCCATCCGGCAGAATGGCGGTAATGTCACTCAGCGCCTGACGATTCTGGGTAATCACGCGCAGATTCGCGCGTTCGCCTTCGGGCTGGGCGGTCGTTTCTTTCAGTTGCAGTTTTTGCTCGCCGCCAAATTTAAAAGCGTCGGAGAGGTAATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCCGCCGGGTAAAGACCAAAGTTATAGGTTTTGCCTGCTTGCCAGGAGCGATCTAACAACACTTCCTGCGCACGTTCAAAAGTCAGCTGG  >  CP009273/4234420‑4234710
                                                                                                                                               |                                                                                                                                                   
ttgttgCCATCCGGCAGAATGGCGGTAATGTCACTCAGCGCCTGACGATTCTGGGTAATCACGCGCGGATTCGCGCGTTCGCCTTCGGGCTGGGCGGTCGTTTCTTTCAGTTGCAGTTTTTGCTCGCCGCCAAATTTAAAAGCGTCGga                                                                                                                                                <  2:357950/149‑1 (MQ=255)
            ggCAGAATGGCGGTAATGTCACTCAGCGCCTGACGATTCTGGGTAATCACGCGCAGATTCGCGCGTTCGCCTTCGGGCTGGGCGGTCGTTTCTTTCAGTTGCAGTTTTTGCTCGCCGCCAAATTTAAAAGCGTCGGAGAGGTAAttttc                                                                                                                                    <  1:153302/149‑2 (MQ=255)
            ggCAGAATGGCGGTAATGTCACTCAGCGCCTGACGATTCTGGGTAATCACGCGCAGATTCGCGCGTTCGCCTTCGGGCTGGGCGGTCGTTTCTTTCAGTTGCAGTTTTTGCTCGCCGCCAAATTTAAAAGCATCGGAGAGGTAAttttt                                                                                                                                    >  2:109978/1‑149 (MQ=255)
            ggCAGAATGGCGGTAATGTCACTCAGCGCCTGACGATTCTGGGTAATCACGCGCAGATTCGCGCGTTCGCCTTCGGGCTGGGCGGTCGTTTCTTTCAGTTGCAGTTTTTGCTCGCCGCCAAATCTAAAAGCATCGGAGAgggacgtttt                                                                                                                                    >  2:463049/1‑141 (MQ=255)
              cAGACTGGCGGTAATGTCGCTCGGCGCCTGCCGCTTCTGGGTAATCACGAGCAGATTCGCGCGTTCGCCTTCGGGCTGGGCGGTCATTTCTTTCAGTTGCAGTTTTTGCTCGCCGCCAAATTTAAAAGCGTCGGAGAGGTAATTTTTGc                                                                                                                                  <  2:181977/149‑1 (MQ=255)
                gAATGGCGGTAATGTCACTCAGCGCCTGGCGATTCTGGGTAATCACGCGCAGATTCGCGCGTTCGCCTTCGGGCTGGGCGGTCGTTTCTTTCAGTTGCAGCGTTTGCTCGCCGACAAATTTAAAAGCGTCGGAGAGGTAATTTTTGCCg                                                                                                                                >  2:763175/1‑149 (MQ=255)
                   tGGCGGTAATGTCACTCAGCGCCTGACGATTCTGGGTAATCACGCGCAGATTCGCGCGTTCGCCTTCGGGCTGGGCGGTCGTTTTTTTCAGTTGCAGTTTTTGCTCGCCGCCAAATTTAAAAGCGTCGGAGGGGTAATTTTTGCCGGtt                                                                                                                             >  1:283353/1‑149 (MQ=255)
                                  tCAGCGCCTGACGATTCTGGGTAATCACGCGCAGATTCGCGCGTTCGCCTTCGGGCTGGGCGGTCGTTTCTTTCAGTTGCAGTTTTTGCTCGCCGCCAAACTTAAAAGCGTCGGAGAGGTAATTTTTGCCGGTTTCGCCGTCGCTGAgc                                                                                                              >  1:822560/1‑149 (MQ=255)
                                    aGCGCCTGACGATTCTGGGTAATCACGCGCAGATTCGCGCGTTCGCCTTCGGGCTGGGCGGTCGTTTCTTTCAGTTGCAGTTTTTGCTCGCCGCCAAATTTAAAAGCGTCGGAGAGGTAATTTTTGCCGGTTTCGCCGTCGCTGAgcgc                                                                                                            >  2:753545/1‑149 (MQ=255)
                                           gACGGTTCTGGGTAATCACGGGCAGATTCGCGCGTTCGCCTTCGGGCTGGGCGGTCGTTTCTTTCAGTTGCAGTTTTTGCTCGCCGCCAAATTTAAAAGCGTCGGAGAGGTAATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGc                                                                                                     <  1:465857/149‑1 (MQ=255)
                                           gACGATTCTGGGTAATCACGCGCAGATTCGCGCGTTCGCCTTCGGGCTGGGCGGTCGTTTCTTTCAGTTGCAGTTTTTGCTCGCCGCCAAATTTAAAAGCGTCGGAGAGGTAATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGc                                                                                                     >  2:14556/1‑149 (MQ=255)
                                             cGATTCTGGGTAATCACGCGCAGATTCGCGCGTTCGCCTTCGGGCTGGGCGGTCGTTTCTTTCAGTTGCAGTTTTTGCTCGCCGCCAAATTTAAAAGCGTCGGAGAGGTAATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCa                                                                                                   <  1:777467/149‑1 (MQ=255)
                                             cGATTCTGGGTAATCACGCGCAGATTCGCGCGTTCGCCTTCGGGCTGGGCGGTCGTTTCTTTCAGTTGCAGTTTTTGCTCGCCGCCAAATTTAAAAGCGTCGGAGAGGTAATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCa                                                                                                   <  1:609500/149‑1 (MQ=255)
                                              gATTCTGGGTAATCACGCGCAGATTCGCGCGTTCGCCTTCGGGCTGGGCGGTCGTTTCTTTCAGTTGCAGTTTTTGCTCGCCGCCAAATTTAAAAGCGTCGGAGAGGTAATTTTTGCCGGTTTCGCCGTCGCTGAGCGACAGTTTCCag                                                                                                  >  2:585799/1‑148 (MQ=255)
                                                 tCTGGGTAATCACGCGCAGATTCGCGCGTTCGCCTTCGGGCTGGGCGGTCGTTTCTTTCAGTTGCAGTTTTTGCTCGCCGCCAAATTTAAAAGCGTCGGAGAGGTAATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCa                                                                                               >  2:788010/1‑149 (MQ=255)
                                                 tCTGGGTAATCACGCGCAGATTCGCGCGTTCGCCTTCGGGCTGGGCGGTCGTTTCTTTCAGTTGCAGGTTTTGCTCGCCGCCAAATTTAAAAGCGTCGGAGAGGTAATTTTTGCCGGCTTCGCCGGCGCTGAGCGCCAGTTGCCCCTCa                                                                                               >  1:121879/1‑149 (MQ=255)
                                                   tGGGTAATCACGCGCAGATTCGCGCGTTCGCCTTCGGGCTGGGCGGTCGTTTCTTTCAGTTGCAGTTTTTGCTCGCCGCCAAATTTAAAAGCGTTGGAGAGGTAATTCTTCAATGCTTGGGCGTCGCTGAACGCCAGTTGCCACTCATc                                                                                             <  2:640124/149‑1 (MQ=255)
                                                          tCACGCGCAGATTCGCGCGTTCGCCTTCGGGCTGGGCGGTCGTTTCTTTCAGTTGCAGTTTTTGCTCGCCGCCAAATTTAAAAGCGTCGGAGAGGTAATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATcgcccgcc                                                                                      >  2:295226/1‑149 (MQ=255)
                                                          tCACGCGCAGATTCGCGCGTTCGCCTTCGGGCTGGGCGGTCGTTTCTTTCAGTTGCAGTTTTTGCTCGCCGCCAAATTTAAAAGCGTCGGAGAGGTAATTTTGGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATcgcccgcc                                                                                      >  2:656222/1‑149 (MQ=255)
                                                             cgcgCAGATTCGCGCGTTCGCCTTCGGGCTGGGCGGTCGTTTCTTTCAGTTGCAGTTTTTGCTCGCCGCCAAATTTAAAAGCGTCGGAGAGGTAATTTTTGCCGGTTTCGCCGTCGCTTAGCGCCAGTTGCCACTCATCGCCCGCCggg                                                                                   <  1:189533/149‑1 (MQ=255)
                                                             cgcgCAGATTCGCGCGTTCGCCTTCGGGCTGGGCGGTCGTTTCTTTCAGTTGCAGTTTTTGCTCGCCGCCAAATTTAAAAGCGTCGGAGAGGTAATTTTTGCCGGTTTCGCCGTCGCCGAGCGCCAGTTGCCACTCATCGCCCGCCggg                                                                                   <  2:852691/149‑1 (MQ=255)
                                                             cgcgCAGATTCGCGCGTTCGCCTTCGGGCTGGACGGTCGTTTCTTTCAGTTGCAGTTTTTGCTCGCCGCCAAATTTAAAAGCGTCGGAGAGGTAATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCCGCCggg                                                                                   <  2:591442/149‑1 (MQ=255)
                                                              gcgcAGATTCGCGCGTTCGCCTTCGGGCTGGGCGGTCGTTTCTTTCAGTTGCAGTTTTTGCTCGCCGCCAAATTTAAAAGCGTCGGAGAGGTAATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCCGCCGGGt                                                                                  >  2:223567/1‑149 (MQ=255)
                                                               cgcAGATTCGCGCGTTCGCCTTCGGGCTGGGCGGTCGTTTCTTTCAGTTGCAGTTTTTGCTCGCCGCCAAATTTAAAAGCGTCGGAGAGGTAATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCCGCCGGGTa                                                                                 <  2:644254/149‑1 (MQ=255)
                                                                        gcgcgTTCGCCTTCGGGCTGGGCGGTCGTTTCTTTCAGTTGCAGTTTTTGCTCGCCGCCAAATTTAAAAGCGTCGGAGAGGTAATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCCGCCGGGTAAAGACCaaa                                                                        >  1:794450/1‑149 (MQ=255)
                                                                         cgcgTTCGCCTTCCGGCTGGGCGGTCGTTTCTTTTAGTTGTAGTTCATACACGCTGTCAAATTTAAAAGCGTCGGAGAGGTAATGTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCCGCCGGGTAAAGACCAAAg                                                                       <  2:24289/149‑1 (MQ=255)
                                                                             ttCGCCTTCGGGCTGGGCGGTCGTTTCTTTCAGTTGCAGTTTTTGCTCGCCGCCAAATTTAAAAGCGTCGGAGAGGTAATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCCGCCGGGTAAAGACCAAAGTtat                                                                   <  2:758476/149‑1 (MQ=255)
                                                                              tCGCCTTCGGGCTGGGCGGTCGTTTCTTTCAGTTGCAGGTTTTGCTCGCCGCCAAATTTAAAAGCGTCGGAGAGGTAATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGCCAGTGCTCGTTGCCCGCCGGGTAAAGACCAAAGTtata                                                                  <  1:404285/149‑1 (MQ=255)
                                                                               cGCCGTCGGGCTGGGCGGTCGTTTCTTTCAGTTGCAGTTTTTGCTCGCCGCCAAATTTAAAAGCGTCGGAGAGGTAATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCCGCCGGGTAAAGACCAAAGTTATAg                                                                 <  1:698305/149‑1 (MQ=255)
                                                                                  cTTCGGGCTGGGCGGTCGCTTCTTTCAGTTGCAGTTTTTGCTCGCCGCCAAATTTAAAAGCGTCGGAGAGGTAAGTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCCGCCGGGTAAAGACCAAAGTTATAGGtt                                                              >  1:203530/1‑149 (MQ=255)
                                                                                      gggCTGGGCGGTCGTTTCTTTCAGTTGCAGTTTTTGCTCGCCGCCAAATTTAAAAGCGTCGGAGAGGTAATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCCGCCGGGTAAAGACCAAAGTTATAGGTTTTGc                                                          <  2:151583/149‑1 (MQ=255)
                                                                                        gCTGGGCGGTCGTTTCTTTCAGTTGCAGTTTTTGCTCGCCGCCAAATTTAAAAGCGTCGGAGAGGTAATTTTTGCCGATTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCCGCCGGGTAAAGACCAAAGTTATAGGTTTTGCCt                                                        <  2:91574/149‑1 (MQ=255)
                                                                                        gCTGGGCGGTCGTTTCTTTCAGCTGCCGTTTTTGCTCGCCGCCAAATTTAAAAGCGGCGGAGAGGTAATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCCGCCGGGTAAAGACCAAAGTTATAGGTTTTGCCt                                                        <  1:281030/149‑1 (MQ=255)
                                                                                         cTGGGCGGTCGTTTCTTTCCGTTCCAGTTTTTGCTCGCCGCCAAATTTAAAAGCGTCGGAGAGGTAATTTTTGCCGGTTTCGCCGTCTCTGAGCGCCAGTTGCCACTCATCGCCCGCCGGGTAAAGACCAAAGTTATAGGTTTTGCCTg                                                       <  1:684796/149‑1 (MQ=255)
                                                                                                        tttcAGTTGCAGTTTCTGCTCGCCGCCAACTTTAAAAGCATCGGAGAGGTAATTTGTGCCGGTTGCGCCGTCGCTGAGCGCCAGTTGCCACTCCTCGCCCGCCGGGTAAAGACCAAAGTTATAGGTTTTGCCTGCTTGCCAGGAGCGAt                                        <  1:463049/149‑1 (MQ=255)
                                                                                                        tttcAGTTGCAGTTGTTGCTCGCCGCCAAATTTAAAAGCATCGGAGAGGTAATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCCGCCGGGTAAAGACCAAAGTTATAGGTTTTGCCTGCTTGCCAGGAGCGAt                                        <  1:109978/149‑1 (MQ=255)
                                                                                                          tcAGTTGCAGTTTTTGCTCGCCGCCAAATTTAAAAGCGTCGGAGAGGTAATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCCGCCGGGTAAAGACCAAAGTTATAGGTTTTGCCTGCTTGCCAGGAGCGATCt                                      <  1:763175/149‑1 (MQ=255)
                                                                                                                   gTTTTTGCTCGCCGCCAAATTTAAAAGCGTCGGAGAGGTAATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCCGCCGGGTAAAGACCAAAGTTATAGGTTTTGCCTGCTTGCCAGGAGCGATCTAACAACACt                             >  2:111111/1‑149 (MQ=255)
                                                                                                                               cgccAAATTTAAAAGCGTCGGAGAGGTACTTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCCGCCGGGTAAAGACCAAAGTTATAGGTTTTGCCTGCTTGCCAGGAGCGATCTAACAACACTTCCTGCGCACGt                 <  1:477104/149‑1 (MQ=255)
                                                                                                                               cgccAAATTTAAAAGCGTCGGAGAGGTAATTTTTGCCGGTTTCGCCGTCGCTGCGCGCCAGTTGCCACTCATCGCCCGCCGGGTAAAGACCAAAGTTATAGGTTTTGCCTGCTTGCCAGGAGCGATCTAACAACACTTCCTGCGCACGt                 <  1:455574/149‑1 (MQ=255)
                                                                                                                                gccAAATTTAAAAGCGTCGGAGAGGTAATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCCGCCGGGTAAAGACCAAAGTTATAGGTTTTGCCTGCTTGCCAGGAGCGATCTAACAACACTTCCTGCGCACGtt                >  2:385225/1‑149 (MQ=255)
                                                                                                                                              tgTCGGGGAGGTAATTTTTGCTGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCCGCCGGGTAAAGACCAAAGTTATAGGTTTTGCCTGCTTGCCAGGAGCGATCTAACAACACTTCCTGCGCACGTTCAAAAGTCAGCTgg  <  1:864008/148‑1 (MQ=255)
                                                                                                                                              cGTCGGAGAGGTAATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCCGCCGGGTAAAGACCAAAGTTATAGGTTTTGCCTGCTTGCCAGGAGCGATCTAACAACACTTCCTGCGCACTTTCAAAAGTCAGCTgg  >  2:211971/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
TTGTTGCCATCCGGCAGAATGGCGGTAATGTCACTCAGCGCCTGACGATTCTGGGTAATCACGCGCAGATTCGCGCGTTCGCCTTCGGGCTGGGCGGTCGTTTCTTTCAGTTGCAGTTTTTGCTCGCCGCCAAATTTAAAAGCGTCGGAGAGGTAATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCCGCCGGGTAAAGACCAAAGTTATAGGTTTTGCCTGCTTGCCAGGAGCGATCTAACAACACTTCCTGCGCACGTTCAAAAGTCAGCTGG  >  CP009273/4234420‑4234710

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: ‑