Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 3,166,744 A→C 20.1% L272V (TTA→GTA)  U712_16260 ← putative ABC transporter permease protein yurN

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068813,166,7440AC20.1% 55.2 / 7.5 38L272V (TTA→GTA) U712_16260putative ABC transporter permease protein yurN
Reads supporting (aligned to +/- strand):  ref base A (22/8);  new base C (2/6);  total (24/14)
Fisher's exact test for biased strand distribution p-value = 3.37e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCACCTCAATGTGCTGGTCAGGCACCGGCTCCAATCTGAACGAATCTGTCTTCTTCTGAGGCAGCATCAGTTACTCCTCCTTTCCGGTTTTAAAGAACTTCATTTGCATGAGCGATAATGCAAGCGTAATGATTAACACGACAGTGGCAATCGCTGAGGCGTAGCCCATCATATCCTTTGTAAAGGCGCTTTTATATAAAAACGTGCTCAAGACTTCTGACGCTTTTCCCGGCCCGCCGCCTGTCAAAATGTACGGCTCATTAAAGACCGT  >  CP006881/3166609‑3166879
                                                                                                                                       |                                                                                                                                       
tCACCTCAATGTGCTGGTCAGGCACCGGCTCCAATCTGAACGAATCTGTCTTCTTCTGAGGCAGCATCAGTTACTCCTCCTTTCCGGTTTTAAAGAACTTCATTTGCATGAGCGATAATGCAAGCGTAATGATTAacac                                                                                                                                      >  5:16785/1‑139 (MQ=255)
       aaTGTGCTGGTCAGGCACCGGCTCCAATCTGAACGAATCTGTCTTCTTCTGAGGCAGCATCAGTTACTCCTCCTTTCCGGTTTTAAAGAACTTCATTTGCATGAGCGATAATGCAAGCGTAATGATTAACACGACAGTg                                                                                                                               >  5:272476/1‑139 (MQ=255)
       aaTGTGCTGGTCAGGCACCGGCTCCAAACTGAACGAATCTGTCTTCTTCTGAGGCAGCATCAGTTACTCCGCCTTTCCGGTTTTAAAGAACTTCATTTGCATGAGCGATAATGCAAGCCTAATGATTAACACGACAGTg                                                                                                                               >  2:11714/1‑139 (MQ=255)
                gTCAGGCACCGGGACCAATCTGAACTAATCTGTCTTCTTCTGAGGCAGCATCAGTTACTCCTCGTTTCCTGTTTAAAAGAAATTCATTTGCATGAGCTATAATGCAAGCGTAATGATTAACACGACAGTGGCAATCGCt                                                                                                                      <  2:262845/139‑1 (MQ=255)
                gTCAGGCACCGGCTCCAATCTGAACGAATCTGTCTTCTTCTGAGGCAGCATCAGTTACTCCTCCTTTCCGGTTTTAAAGAACTTCATTTGCATGAGCGATAATGCACGCGTAATGATTAACACGACAGTGGCAATCGCt                                                                                                                      >  5:227506/1‑139 (MQ=255)
                gTCAGGCACCGGCTCCAATCTGAACGAATCTGTCTTCTTCTGAGGCAGCATCAGTTACTCCTCCTTTCCGGTTTTAAAGAACTTCATTTGCATGAGCGATAATGCAAGCGTAATGATTAACACGACAGTGGCAATCGCt                                                                                                                      >  1:72711/1‑139 (MQ=255)
                     gCACCGGCTCCAATCTGAACGAATCTGTCTTCTTCTGAGGCAGCATCAGTTACTCCTCCTTTCCGGTTTTAAAGAACTTCATTTGCATGAGCGATAATGCAAGCGTAATGATTAACACGACAGTGGCAATCGCTGAGGc                                                                                                                 >  1:237647/1‑139 (MQ=255)
                          ggCTCCACTCGGACCCAAGCTGGCTTCTTCTGAGGCACCAGCCGTAACTCCTCCTTGCCGGTTTAAAAGAACTTCATTTGCATGAGCTATACTGCAAGCGTAATTATAAACACGACAGTGGCAATCGCTGAGGCGTAGc                                                                                                            <  6:282914/139‑1 (MQ=255)
                          ggCTCCAATCTGAAGTAAGCTGTCTTCTGCTGAGGCCGCATCTGTTACTCCTCCTTGCCGGTTTTAAATAACTTCATTTGCTGGTGCTATAATGCACGCTTAATTATAACCACGACAGTGGCAATCGCTGAGGCGTAGc                                                                                                            <  6:227506/139‑1 (MQ=255)
                                     cacataaTCTGTTTTGTTGTGAGGCAGCATCAGTTACTCCTCCTTTCCGGTTTAAAATAACTTCATTTGATTGAGCTATAATGCATGCTTAATTATTAACACGACAGTGGCAATCGCTGAGGCGTAGCCCATCATATc                                                                                                  <  6:268030/133‑1 (MQ=255)
                                      aaCGAATCTGTCTTCTTCTGAGGCAGCATCAGTTACTCCTCCTTTCCGGTTTTAAAGAACTTCATTTGCATGAGCGATAATGCAAGCGTAATGATTAACACGACAGTGGCAATCTCTGAGGCTTAGCCCATCATATcat                                                                                                >  8:133023/1‑137 (MQ=255)
                                      aaCGAATCTGTCTTCTTCTGAGGCAGCATCAGTTACTCCTCCTTTCCGGTTTTAAAGAACTTCATTTGCATGAGCGATAATGCAAGCGTAATGATTAACACGACAGTGCCACTCGCTGAGGCGTAGCCCATCATATCCt                                                                                                >  2:208168/1‑139 (MQ=255)
                                      aaCGAATCTGTCTTCTTCTGAGGCAGCATCAGTTACTCCTCCTTTCCGGTTTTAAAGAACTTCATTTGCATGAGAGATAATGCAAGCGTAATGATTAACACGACAGTGGCAATCGCTGAGGCGTAGCCCATCATATcat                                                                                                >  8:128098/1‑137 (MQ=255)
                                           aTCTGTCTTCTTCTGAGGCAGCATCAGTTACTCCTCCTTTCCGGTTTTAAAGAACTTCATTTGCATGAGCGATAATGCAAGCGTAATGATTAACACGACAGTGGCAATCGCTGAGGCGTAGCCCATCATTTCCTTGGt                                                                                            >  4:8693/1‑138 (MQ=255)
                                                     ttctGAGGCAGCATCAGTTACTCCTCCTTTCCGGTTTTAAAGAACTTCATTTGCATGAGCGATAATGCAAGCGTAATGATTAACACGACAGTGGCAATCGCTGAGGCGTAGCCCATCATATCCTTTGTAAAGgctgtct                                                                                 >  5:268030/1‑134 (MQ=255)
                                                       ctgcgccagcaCCAGATACACCTCCTTACCGGTTTTAAACAAATTAATTTGCTTGCGCTATAATCCACGCTTATTCATAACCACGACAGTGGCAATCGCTGAGGCGTAGCCCATCATATCCTTTGTAAAGGCGCtt                                                                                  <  2:52102/130‑1 (MQ=255)
                                                              agcaTCAGTTACTCCTCCTTTCCGGTTTTAAAGAACTTCATTTGCATGAGCGATAATGCAAGCGTAATGATTAACACGACAGTGGCAATCGCTGAGGCGTAGCCCATCATATCCTTTGTAAAGGCGCTTTTATATaaa                                                                         >  2:136902/1‑138 (MQ=255)
                                                               gcaTCAGTTACTCCTCCTTTCCGGTTTTAAAGAACTTCATTTGCATGAGCGATAATGCATGAGTAATGATTAACACGACAGTGGCAATCGCTGAGGCGTAGCCCATCATATCCTTAGTAACTGCGCTTTTATATaaagg                                                                       >  4:41521/1‑137 (MQ=255)
                                                                   cAGTTACTCCTCATTTCCGTTTTTAACCAAATTAATTTGCATGAGCTAAAATCCAAGCTTAATTATAACCACGACAGTGGCAATCGCTGAGGCGTAGCCCATCATATCCTTTGTAAAGGCGCTTTTATATAAAAACGTg                                                                   <  4:290654/139‑1 (MQ=255)
                                                                          tcctccTTTCCGGTTTTAAAGAACTTCATTTGCATGAGCGATAATGCAAGCGTAATGATTACCACGACAGTGGCAATCGCTGAGGCGTAGCCCATCAGATCCTTTGTAAAGGCGCTTTTATATAACAACGTGCTCAAGa                                                            >  8:75306/1‑139 (MQ=255)
                                                                          tcctccTTTCCGGTTTTAAAGAACTTCATTTGCATGAGCGATAATGCAAGCGTAATGATTAACACGACAGTGGCAATCGCTGAGGCGTAGCCCATCATATCCTTTGTAAAGGCGCTTTTATATAAAAACGTGCTCAAGa                                                            >  7:353926/1‑139 (MQ=255)
                                                                          tcctccTTTCCGGTTTTAAAGAACTTCATTTGCATGAGCGATAATGCAAGCGTAATGATTAACACGACAGTGGCAATCCCTGAGGCGTAGCCCATCATATCCTTTGTAAAGGCGCTTTTATATAAAAACGTGCTCAAGa                                                            >  8:371593/1‑139 (MQ=255)
                                                                              cctTTCCGGTTTTAAAGAACTTCATTTGCATGAGCGATAATGCAAGCGTAATGATTAACACGACAGTGGCAATCGCTGAGGCGTAGCCCATCATATCCTTTGTAAAGGCGCTTTTATATAAAAACGTGCTCAAGACTTc                                                        >  6:356644/1‑139 (MQ=255)
                                                                              cctTTCCGGTTTTAAAGAACTTCATTTGCATGAGCGATAATGCAAGCGTAATGATTAACACGACAGTGGCAATCGCTGAGGCGTAGCCCATCATATCCTTTGTAAAGGCGCTTTTATATAAAAACGTGCTCAAGACTTc                                                        >  7:306725/1‑139 (MQ=255)
                                                                                 ttgcgctttaAACGAACTTCATTTCAATGAGCGATCAGCAAAGCTTAATTAATAACACGACAGGGGCAATCGCTGAGGCGTAGCCCATCATATCCTTTGTAAAGGCGCTTTTATATAAAAACGTGCTCAAGACTTCt                                                       <  6:100121/131‑1 (MQ=255)
                                                                                   ccGGTTTTAAAGAACTTCATTTGCATGAGCGATAATGCAAGCGTAATGATTAACACGACAGTGGCAATCGCTGAGGCGTAGCCCATCATATCCTTTGTAAAGGCGCTTTTATATAAAAACGTGCTCAAGACTTCTGACg                                                   >  6:288781/1‑139 (MQ=255)
                                                                                                          cATGAGCGATAATGCAAGCGTAATGATTACCACGACAGTGGCAATCGCTGAGGCGTAGCCCATCATATCCTTTGTAAAGGCGCTTTTATATTAAAACGTGCTCAAGACTTCTGACGCTTTTCCCTGCCCGCCCCCTGt                             >  5:49409/1‑138 (MQ=255)
                                                                                                          cATGAGCGATAATGCAAGCGTAATGATTAACACGACAGTGGCAATCGCTGAGGCGTAGCCCATCATATCCTTTGTAAAGGCGCTTTTATATAAAAACGTGCTCAAGACTTCTGACGCTTTTCCCGGCCCGCCGCCTGt                             >  7:67807/1‑138 (MQ=255)
                                                                                                           aTGAGCGATAATGCAAGCGTAATGATAAACACGACAGTGGCAATCGCTGAGGCGTAGCCCATCATATCCTTTGTAAAGGCGCTTTTATATAAAAACGTGCTCAAGACTTCTGACGCTTTTCCCGGCCCGCCGCCTGTCa                           <  3:177271/139‑1 (MQ=255)
                                                                                                              aGCGATAATGCAAGCGTAATGATTAACACGACAGTGGCAATCGCTGAGGCGTAGCCCATCATATCCTTTGTAAAGGCGCTTTTATATAAAAACGTGCTCAAGACTTCTGACGCTTTTCCCGGCCCGCCGCCTGTCaaaa                        >  7:188074/1‑139 (MQ=255)
                                                                                                                  aTAATGAACGCTTAATGATAACCACGACAGTGGCAATCGCTGAGGCGTAGCCCATCATATCCTTTGTAAAGGCGCTTTTATATAAAAACGTGCTCAAGACTTCTGACGCTTTTCCCGGCCCGCCGCCTGTCAAAATGTa                    <  5:345168/139‑1 (MQ=255)
                                                                                                                  aTAATCCAAGCTTAATTATAAACACGACAGTGGCAATCGCTGAGGCGTAGCCCATCATATCCTTTGTAAAGGCGCTTTTATATAAAAACGTGCTCAAGACTTCTGACGCTTTTCCCGGCCCGCCGCCTGTCAAAATGTa                    <  7:291109/139‑1 (MQ=255)
                                                                                                                      tGCAAGCGTAATGATTAACACGACAGTGGCAATCGCTGAGGCGTAGCCCATCATATCCTTTGTAAAGGCGCTTTTATATAAAAACGTGCTCAAGACTTCTGACGCTTTTCCCGGCCCGCCGCCTGTCAAAATGTACGGc                <  3:135786/139‑1 (MQ=255)
                                                                                                                               aaTGATTACCACGACAGTGGCAATCGCTGAGGCGTAGCCCATCATATCCTTTGTAAAGGCGCTTTTATATAAAAACGTGCTCAAGACTTCTGACGCTTTTCCCGGCCCGCCGCCTGTCAAAATGTACGGCTCATTaaa        <  7:295778/138‑1 (MQ=255)
                                                                                                                               aaTGATTAACACGACAGTGGCAATCGCTGAGGCGTAGCCCATCATATCCTTTGTAAAGGCGCTTTTATATAAAAACGTGCTCAAGACTTCTGACGCTTTTCCCGGCCCGCCGCCTGTCAAAATGTACGGCTCATTaaa        <  1:136902/138‑1 (MQ=255)
                                                                                                                                aTGATTAACACGACAGTGGCAATCGCTGAGGCGTAGCCCATCATATCCTTTGTAAAGGCGCTTTTATATAAAAACGTGCTCAAGACTTCTGACGCTTTTCCCGGCCCGCCGCCTGTCAAAATGTACGGCTCATTAAAGa      >  7:369293/1‑139 (MQ=255)
                                                                                                                                aTGATTAACACGACAGTGGCAATCGCTGAGGCGTAGCCCATCATATCCTTTGTAAAGGCGCTTTTATATAAAAACGTGCTCAAGACTTCTGACGCTTTTCCCGGCCCGCCGCCTGTCAAAATGTACGGCTCATTAAAGa      >  6:296147/1‑139 (MQ=255)
                                                                                                                                aTGATTAACACGACAGTGGCAATCGCTGAGGCGTAGCCCATCATATCCTTTGTAAAGGCGCTTTTATATAAAAACGTGCTCAAGACTTCTGACGCTTTTCACGGCCCTCCGCATGTCAAAATGTACGGCTCATTAAAGa      >  6:110793/1‑139 (MQ=255)
                                                                                                                                 tGATAACCACGACAGTGGCAATCGCTGAGGCGTAGCCCATCATATCCTTTGTAAAGGCGCTTTTATATAAAAACGTGCTCAAGACTTCTGACGCTTTTCCCGGCCCGCCGCCTGTCAAAATGTACGGCTCATTAAAGAc     <  6:272476/139‑1 (MQ=255)
                                                                                                                                      aaCACGACAGTGGCAATCGCTGAGGCGTAGCCCATCATATCCTTTGTAAAGGCGCTTTTATATAAAAACGTGCTCAAGACTTCTGACGCTTTTCCCGGCCCGCCGCCTGTCAAAATGTACGGCTCATTAAAGACCGt  >  1:355526/1‑137 (MQ=255)
                                                                                                                                       acacGACAGTGGCAATCGCTGAGGCGTAGCCcatcat                                                                                                     >  7:296033/1‑37 (MQ=255)
                                                                                                                                       acacGACAGTGGCAATCGCTGAGGCGTAGCCcatcat                                                                                                     <  8:296033/37‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                       
TCACCTCAATGTGCTGGTCAGGCACCGGCTCCAATCTGAACGAATCTGTCTTCTTCTGAGGCAGCATCAGTTACTCCTCCTTTCCGGTTTTAAAGAACTTCATTTGCATGAGCGATAATGCAAGCGTAATGATTAACACGACAGTGGCAATCGCTGAGGCGTAGCCCATCATATCCTTTGTAAAGGCGCTTTTATATAAAAACGTGCTCAAGACTTCTGACGCTTTTCCCGGCCCGCCGCCTGTCAAAATGTACGGCTCATTAAAGACCGT  >  CP006881/3166609‑3166879

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