Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 119,473 A→G 9.8% H380H (CAT→CAC)  yjbB ← possible PnaS family phosphate:sodium (Na+) symporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730119,4730AG 9.8% 127.0 / 4.5 51H380H (CAT→CAC) yjbBpossible PnaS family phosphate:sodium (Na+) symporter
Reads supporting (aligned to +/- strand):  ref base A (29/17);  new base G (2/3);  total (31/20)
Fisher's exact test for biased strand distribution p-value = 3.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

GTTAGCTCTGCTACTTTTAAAATAGAGCGATTGACATCAAACATAACTGCTAAACGCTCAACGTCTGCCTTCGTAATGGCTTTTGTAGAAATTCTAACTAAATAATTTCGAATGCTATCATTGATTGTTTCAACAGCTTGATGCTTTTGTTCAAGCTTTTTGATCAATTTTTTATCGTCTTTTGTAATTTCGCGAATGTCTTCAAACATTGATAAGACAATCTGACCCACATTTTGTAATTCTTTTTGAGTTTCTTGTAATGCAACACCAGGTGCGTGATA  >  CP000730/119333‑119613
                                                                                                                                            |                                                                                                                                            
gTTAGCTCTGCTACTTTTAAAATAGAGCGATTGACATCAAACATAACTGCTAAACGCTCAACGTCTGCCTTCGTAATGGCTTTTGTAGAAATTCTAACTAAATAATTTCGAATGCTATCATTGATTGTTTCAACAGCTTGa                                                                                                                                              <  1:261722/141‑1 (MQ=255)
  tAGCTCTGCTACTTTTAAAATAGAGCGATTGACATCAAACATAACTGCTAAACGCTCAACGTCTGCCTTCGTAATGGCTTTTGTAGAAATTCTAACTAAATAATTTCGAATGCTATCATTGATTGTTTCAACAGCTTGATg                                                                                                                                            <  2:104840/141‑1 (MQ=255)
        tGCTACTTTTAAAATAGAGCGATTGACATCAAACATAACTGCTAAACGCTCAACGTCTGCCTTCGTAATGGCTTTTGTAGAAATTCTAACTAAATAATTTCGAATGCTATCATTGATTGTTTCAACAGCTTGATGCTTTTg                                                                                                                                      >  1:312255/1‑141 (MQ=255)
             cTTTTAAAATAGAGCGATTGACATCAAACATAACTGCTAAACGCTCAACGTCTGCCTTCGTAATGGCTTTTGTAGAAATTCTAACTAAATCATTTCGACTGCTATCATTGATTGTTTCAACAGCTTGATGCTTTTGTTCaa                                                                                                                                 >  1:87019/1‑141 (MQ=255)
             cTTTTAAAATAGAGCGATTGACATCAAACATAACTGCTAAACGCTCAACGTCTGCCTTCGTAATGGCTTTTGTAGAAATTCTAACTAAATAATTTCGAATGCTATCATTGATTGTTTCAACAGCTTGATGCTTTTGTTCaa                                                                                                                                 >  2:451923/1‑141 (MQ=255)
              ttttAAAATAGAGCGATTGACATCAAACATAACTGCTAAACGCTCAACGTCTGCCTTCGTAATGGCTTTTGTAGAAATTCTAACTAAATAATTTCGAATGCTATCATTGATTGTTTCAACAGCTTGATGCTTTTGTTCAAg                                                                                                                                <  2:161654/141‑1 (MQ=255)
                       agagCGATTGACATCAAACATAACTGCTAAACGCTCAACGTCTGCCTTCGTAATGGCTTTTGTAGAAATTCTAACTAAATAATTTCGAATGCTATCATTGATTGTTTCAACAGCTTGATGCTTTTGTTCAAGCTTTTTGAt                                                                                                                       >  1:64180/1‑141 (MQ=255)
                        gagCGATTGACATCAAACATAACTGCTAAACGCTCAACGTCTGCCTTCGTAATGGCTTTTGTAGAAATTCTAACTAAATAATTTCGAATGCTATCATTGATTGTTTCAACAGCTTGATGCTTTTTTTCAAGCTTTTTGATc                                                                                                                      >  2:183059/1‑141 (MQ=255)
                        gagCGATTGACATCAAACATAACTGCTAAACGCTCAACGTCTGCCTTCGTAATGGCTTTTGTAGAAATTCTAACTAAATAATTTCGAATGCTATCATTGATTGTTTCAACAGCTTGATGCTTTTGTTCAAGCTTTTTGATc                                                                                                                      >  1:361704/1‑141 (MQ=255)
                        gagCGATTGACATCAAACATAACTGCTAAACGCTCAACGTCTGCCTTCGTAATGGCTTTTGTAGAAATTCTAACTAAATAATTTCGAATGCTATCATTGATTGTTTCAACAGCTTGATGCTTTTGTTCAAGCTTTTTGATc                                                                                                                      <  1:32039/141‑1 (MQ=255)
                                  cATCAAACATAACTGCTAAACGCTCAACGTCTGCCTTCGTAATGGCTTTTGTAGAAATTCTAACTAAATAATTTCGAATGCTATCATTGATTGTTTCAACAGCTTGATGCTTTTGTTCAAGCTTTttttttatttttttat                                                                                                            >  2:274397/1‑141 (MQ=255)
                                    tCAAACATAACTGCTAAACGCTCAACGTCTGCCTTCGTAATGGCTTTTGTAGAAATTCTAACTAAATAATTTCGAATGCTATCATTGATTGTTTCAACAGCTTGATGCTTTTGTTCAAGCTTTTTGATCAATTTTTTATCg                                                                                                          >  2:339503/1‑141 (MQ=255)
                                      aaaCATAACTGCTAAACGCTCAACGTCTGCCTTCGTAATGGCTTTTGTAGAAATTCTAACTAAATAATTTCGAATGCTATCATTGATTGTTTCAACAGCTTGATGCTTTTGTTCAAGCTTTTTGATCAATTTTTTATCGTc                                                                                                        >  1:226533/1‑141 (MQ=255)
                                       aaCATAACTGCTAAACGCTCAACGTCTGCCTTCGTAATGGCTCTTGTAGAAATTCTAACTAAATAATTTCGAATGCTATCATTGATTGTTTCAACAGCTTGATGCTTTTGTTCAAGCTTTTTGATCAATTTTTTATCGTCt                                                                                                       >  2:205167/1‑141 (MQ=255)
                                          aTAACTGCTAAACGCTCAACGTCTGCCTTCGTAATGGCTTTTGTAGAAATTCTAACTAAATAATTTCGAATGCTATCATTGATTGTTTCAACAGCTTGATGCTTTTGTTCAAGCTTTTTGATCAATTTTTTATCGTCtttt                                                                                                    >  1:408119/1‑141 (MQ=255)
                                                      cGCTCAACGTCTGCCTTCGTAATGGCTTTTGTAGAAATTCTAACTAAATAATTTCGAATGCTATCATTGATTGTTTCAACAGCTTGGTGCTTTTGTTCAAGCTTCTTGATCAATTTTTTATCGTCTTTTGTAATTTCGCGa                                                                                        <  2:416203/141‑1 (MQ=255)
                                                      cGCTCAACGTCTGCCTTCGTAATGGCTTTTGTAGAAATTCTAACTAAATAATTTCGAATGCTATCATTGATTGTTTCAACAGCTTGATGCTTTTGTTCAAGCTTTTTGATCAATTTTTTATCGTCTTTTGTAATTTCGCGa                                                                                        <  2:354688/141‑1 (MQ=255)
                                                        cTCAACGTCTGCCTTCGTAATGGCTTTTGTAGAAATTCTAACTAAATAATTTCGAATGCTATCATTGATTGTTTCAACAGCTTGATGCTTTTGTTCAAGCTTTTTGATCAATTTTTTATCGTCTTTTGTAATTTCGCGAAt                                                                                      >  1:389328/1‑141 (MQ=255)
                                                        cTCAACGTCTGCCTTCGTAATGGCTTTTGTAGAAATTCTAACTAAATAATTTCGAATGCTATCATTGATTGTTTCAACAGCTTGATGCTTTTGTTCAAGCTTTTTGATCAATTTTTTATCGTCTTTTGTAATTTCGCGAAt                                                                                      >  2:369840/1‑141 (MQ=255)
                                                         tCAACGTCTGCCTTCGTAATGGCTTTTGTAGAAATTCTAACTAAATAATTTCGAATGCTATCATTGATTGTTTCAACAGCTTGATGCTTTTGTTCAAGCTTTTTGATCAATTTTTTATCGTCTTTTGTAATTTCGCGAATg                                                                                     <  2:8168/141‑1 (MQ=255)
                                                             cggcTGCCTTCGTAATGGCCATTGTATAAATTCTAACTAAATAGTTTCGAATGGTATCAGTCATCGTTTCAACACCTGGATGCTTTTGTTCAAGCTCTTCGAGCAATTTTTTATCGTCTTTTGTAATGTAGCTCATGTCtt                                                                                 >  2:358374/4‑141 (MQ=255)
                                                             cGTCTGCCTTCGTAATGGCTTTTGTAGAAATTCTAACTAAATAATTTCGAATGCTATCATTGATTGTTTCAACAGCTTGATGCTTTTGTTCAAGCTTTTTGATCAATTTTTTATCGTCTTTTGTAATTTCGCGAATGTCtt                                                                                 >  1:35209/1‑141 (MQ=255)
                                                             cGTCTGCCTTCGTAATGGCTTTTGTAGAAATTCTAACTAAATAATTTCGAATGCTATCATTGATTGTTTCAACAGCTTGATGCTTTTGTTCAAGCTTTTTGATCAATTTTTTATCGTCTTTTGTAATTTCGCGAATGTCtt                                                                                 >  1:351703/1‑141 (MQ=255)
                                                             cGTCTGCCTTCGTAATGGCTTTTGTAGAAATTCTAACTAAATAATTTCGAATGCTATCATTGATTGTTTCAACAGCTTGATGCTTTTGTTCAAGCTTTTTGATCAATTTTTTATCGTCTTTTGTAATTTCGCGAATGTCtt                                                                                 >  2:187276/1‑141 (MQ=255)
                                                              gTCTGCCTTCGTAATGGCTTTTGTAGAAATTCTAACTAAATAATTTCGAATGCTATCATTGATTGTTTCAACAGCTTGATGCTTTTGTTCAAGCTTTTTGATCAATTTTTTATCGTCTTTTGTAATTTCGCGAATGTCTTc                                                                                <  1:350546/141‑1 (MQ=255)
                                                               tCTGCCTTCGTAATGGCTTTTGTAGAAATTCTAACTAAATAATTTCGAATGCTATCATTGATTGTTTCAACAGCTTGATGCTTTTGTTCAAGCTTTTTGATCAATTTTTTATCGTCTTTTGTAATTTCGCGAATGTCTTCa                                                                               <  2:409117/141‑1 (MQ=255)
                                                               tCTGCCTTCGTAATGGCTTTTGTAGAAATTCTAACTAAATAATTTCGAATGCTATCATTGATTGTTTCAACAGCTTGATGCTTTTGTTCAAGCTTTTTGATCAATTTTTTATCGTCTTTTGTAATTTCGCGAATGTCTTCa                                                                               <  2:423009/141‑1 (MQ=255)
                                                               tCTGCATTCGTAATGGCTTTTGTAGAAATTCTAACTAAATAATTTCGAATGCTATCATTGATTGTTTCAACAGCTTGATGCTTTTGTTCAAGCTTTTTGATCAATTTTTTATCGTCTTTTGTAATTTCGCGAATGTCTTCa                                                                               <  2:87019/141‑1 (MQ=255)
                                                                 tGCCTTCGTAATGGCTTTTGTAGAAATTCTAACTAAATAATTTCGAATGCTATCATTGATTGTTTCAACAGCTTGATGCTTTTGTTCAAGCTTTTTGATCAATTTTTTATCGTCTTTTGTAATTTCGCGAATGTCTTCaaa                                                                             <  1:77243/141‑1 (MQ=255)
                                                                   ccTTCGTAATGGCTTTTGTAGAAATTCTAACTAAATAATTTCGAATGCTATCATTGATTGTTTCAACAGCTTGGTGCTTTTGTTCAAGCTTCTTGATCAATTTTTTATCGTCTTTTGTAATTTCGCGAATGTCTTCAAACa                                                                           <  2:339587/141‑1 (MQ=255)
                                                                            tGGCTTTTGTAGAAATTCTAACTAAATAATTTCGAATGCTATCATTGATTGTTTCAACAGCTTGATGCTTTTGTTCAAGCTTTTTGATCAATTTTTTATCGTCTTTTGTAATTTCGCGAATGTCTTCAAACATAGATAAGa                                                                  >  1:412557/1‑141 (MQ=255)
                                                                              gCTTTTGTAGAAATTCTAACTAAATAATTTCGAATGCTATCATTGATTGTTTCAACAGCTTGATGCTTTTGTTCAAGCTTTTTGATCAATTTTTTATCGTCTTTTGTAATTTCGCGAATGTCTTCAAACATTGATAAGACa                                                                >  2:62309/1‑141 (MQ=255)
                                                                                  ttGTAGAAATTCTAACTAAATAATTTCGAATGCTATCATTGATTGTTTCAACAGCTTGATGCTTTTGTTCAAGCTTTTTGATCAATTTTTTATCGTCTTTTGTAATTTCGCGAATGTCTTCAAACATTGATAAGACAATCt                                                            >  2:356269/1‑141 (MQ=255)
                                                                                   tGTAGAAATTCTAACTAAATAATTTCGAATGCTATCATTGATTGTTTCAACAGCTTGGTGCTTTTGTTCAAGCTTCTTGATCAATTTTTTATCGTCTTTTGTAATTTCGCGAATGTCTTCAAACATTGATAAGACAATCTg                                                           >  2:381295/1‑141 (MQ=255)
                                                                                   tGTAGAAATTCTAACTAAATAATTTCGAATGCTATCATTGATTGTTTCAACACCTTGGTGCTTTTGTTCAAGCTTCTTGATCAATTTTTTATCGTCTTTTGTAATTTCGCGAATGTCTTCAAACATTGATAAGACAATCTg                                                           >  2:64239/1‑141 (MQ=255)
                                                                                       gAAATTCTAACTAAATAATTTCGAATGCTATCATTGATTGTTTCAACAGCTTGATGCTTTTGTTCAAGCTTTTTGATCAATTTTTTATCGTCTTTTGTAATTTCGCGAATGTCTTCAAACATTGATAAGACAATCTGAccc                                                       >  1:290966/1‑141 (MQ=255)
                                                                                         aaTTCTAACTAAATAATTACGAATGCTATCATTGATTGTTTCAACACCTTGATACTTTTGTTCAAACTTTTTGAACAATTTTTTATCGTCTTTTGTTATTGCACGATTGTCTTAAAACATTTCTATGACACTTTGAccatc                                                     >  2:402205/1‑138 (MQ=255)
                                                                                             ctaactaaATAATTTCGAATGCTATCATTGATTGTTTCAACAGCTTGATGCTTTTGTTCAAGCTTTTTGATCAATTTTTTATCGTCTTTTGTAATTTCGCGAATGTCTTCAAACATTGATAAGACAATCTGACCCACAttt                                                 >  2:180201/1‑141 (MQ=255)
                                                                                                actaaATAATTTCGAATGCTATCATTGATTGTTTCAACAGCTTGATGCTTTTGTTCAAGCTTTTTGATCAATTTTTTATCGTCTTTTGTAATTTCGCGAATGTCTTCAAACATTGATAAGACAATCTGACCCACATTTTGt                                              >  2:199962/1‑141 (MQ=255)
                                                                                                actaaATAATTTCGAATGCTATCATTGATTGTTTCAACAGCTTGATGCTTTTGTTCAAGCTTTTTGATCAATTTTTTATCGTCTTTTGTAATTTCGCGAATGTCTTCAAACATTGATAAGACAATCGGACCCACATTTTGt                                              <  2:422417/141‑1 (MQ=255)
                                                                                                   aaaTAATTTCGAATGCTATCATTGATTGTTTCAACAGCTTGATGCTTTTGTTCAAGCTTTTTGATCAATTTTTTATCGTCTTTTGTAATTTCGCGAATGTCTTCAAACATTGATAAGACAATCTGACCCACATTTTGTAAt                                           <  2:312255/141‑1 (MQ=255)
                                                                                                   aaaTAATTTCGAATGCTATCATTGATTGTTTCAACAGCTTGATGCTTTTGTTCAAGCTTTTTGATCAATTTTTTATCGTCTTTTGTAATTTCGCGAATGTCTTCAAACATTGATAAGACAATCTGACCCACATTTTGTAAt                                           <  2:80688/141‑1 (MQ=255)
                                                                                                     ataatTTCGAATGCTATCTTTCATTGTTTCAACAGCTTGATGCTTTTGTTCAAGCTTTTTGATCAATTTTTTATCGTCTGTTGTAATTTCGCGAATGTCTTCAAACATTGATAAGACAATCTGACCCACATTTTGTAATTc                                         <  1:115625/141‑1 (MQ=255)
                                                                                                        atTTCGAATGCTATCATTGATTGTTTCAACAGCTTGGTGCTTTTGTTCAAGCTTCTTGATCAATTTTTTATCGTCTTTTGTAATTTCGCGAATGTCTTCAAACATTGATAAGACAATCTGACCCACATTTTGTAATTCttt                                      <  2:402376/141‑1 (MQ=255)
                                                                                                           tCGAATGCTATCATTGATTGTTTCAACAGCTTGATGCTTTTGTTCAAGCTTTTTGATCAATTTTTTATCGTCTTTTGTAATTTCGCGAATGTCTTCAAACATTGATAAGACAATCTGACCCACATTTTGTAATTCTTTTTg                                   >  1:141432/1‑141 (MQ=255)
                                                                                                              aaTGCTATCATTGATTGTTTCAACAGCTTGATGCTTTTGTTCAAGCTTTTTGATCCAATTTTTATCATCTTTTGTAATTAGGCGAATGTCTTCAAACATTGATAAGACAATCTGACCCACATTTTGTAATTCTTTTTGAGt                                <  1:200559/141‑1 (MQ=255)
                                                                                                               aTGCTATCATTGATTGTTTCAACAGCTTGATGCTTTTGTTCAAGCTTTTTGATCAATTTTTTATCGTCTTTTGTAATTTCGCGAATGTCTTCAAACATTGATAAGACAATCTGACCCACATTTTGTAATTCTTTTTGAGtt                               <  1:420390/141‑1 (MQ=255)
                                                                                                                   tATCATTGATTGTTTCAACAGCTTGATGCTTTTGTTCAAGCTTTTTGATCAATTTTTTATCGTCTTTTGTAATTTCGCGAATGTCTTCAAACATTGATAAGACAATCTGACCCACATTTTGTAATTCTTTTTGAGTTTCtt                           <  1:50692/141‑1 (MQ=255)
                                                                                                                       attgattgTTTCAACAGCTTGATGCTTTTGTTCAAGCTTTTTGATCAATTTTTTATCGTCTTTTGTAATTTCGCGAATGTCTTCAAACATTGATAAGACAATCTGACCCACATTTTGTAATTCTTTTTGAGTTTCTTGTaa                       >  1:21604/1‑141 (MQ=255)
                                                                                                                            ttgTTTCAACAGCTTGATGCTTTTGTTCAAGCTTTTTGATCAATTTTTTATCGTCTCTTGTAATTTCGCGAATGTCTTCAAACATTGATAAGACAATCTGACCCACATTTTGTAATTCTTTTTGAGTTTCTTGTAATGCaa                  >  2:431509/1‑141 (MQ=255)
                                                                                                                                 tCAACAGCTTGATGCTTTTGTTCAAGCTTTTTGATCAATTTTTTATCGTCTTTTGTAATTTCGCGAATGTCTTCAAACATTGATAAGACAATCTGACCCACATTTTGTAATTCTTTTTGAGTTTCTTGTAATGCAACACCa             >  1:295121/1‑141 (MQ=255)
                                                                                                                                  cAACAGCTTGATGCTTTTGTTCAAGCTTTTTGATCAATTTTTTATCGTCTTTTGTAATTTCGCGAATGTCTTCAAACATTGATAAGACAATCTGACCCACATTTTGTAATTCTTTTTGAGTTTCTTGTAATGCAACACCAg            >  1:241867/1‑141 (MQ=255)
                                                                                                                                            aTGCTTTTGTTCAAGCTTTTTGATCAATTTTTTATCGTCTTTTGTAATTTCGCGAATGTCTTCAAACATTGATAAGACAATCTGACCCACATTTTGTAATTCTTTTTGAGTTTCTTGTAATGCAACACCAGGTGCGTGATa  >  2:387361/1‑141 (MQ=255)
                                                                                                                                            |                                                                                                                                            
GTTAGCTCTGCTACTTTTAAAATAGAGCGATTGACATCAAACATAACTGCTAAACGCTCAACGTCTGCCTTCGTAATGGCTTTTGTAGAAATTCTAACTAAATAATTTCGAATGCTATCATTGATTGTTTCAACAGCTTGATGCTTTTGTTCAAGCTTTTTGATCAATTTTTTATCGTCTTTTGTAATTTCGCGAATGTCTTCAAACATTGATAAGACAATCTGACCCACATTTTGTAATTCTTTTTGAGTTTCTTGTAATGCAACACCAGGTGCGTGATA  >  CP000730/119333‑119613

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑