Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 119,491 T→C 9.6% K374K (AAA→AAG)  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,4910TC 9.6% 131.1 / 6.1 52K374K (AAA→AAG) yjbBpossible PnaS family phosphate:sodium (Na+) symporter
Reads supporting (aligned to +/- strand):  ref base T (30/17);  new base C (2/3);  total (32/20)
Fisher's exact test for biased strand distribution p-value = 3.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01

AGAGCGATTGACATCAAACATAACTGCTAAACGCTCAACGTCTGCCTTCGTAATGGCTTTTGTAGAAATTCTAACTAAATAATTTCGAATGCTATCATTGATTGTTTCAACAGCTTGATGCTTTTGTTCAAGCTTTTTGATCAATTTTTTATCGTCTTTTGTAATTTCGCGAATGTCTTCAAACATTGATAAGACAATCTGACCCACATTTTGTAATTCTTTTTGAGTTTCTTGTAATGCAACACCAGGTGCGTGATAAAC‑AAGATCTTTGTTTA  >  CP000730/119356‑119630
                                                                                                                                       |                                                                                                                                            
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                                                                                                        >  2:369840/1‑141 (MQ=255)
                                 cTCAACGTCTGCCTTCGTAATGGCTTTTGTAGAAATTCTAACTAAATAATTTCGAATGCTATCATTGATTGTTTCAACAGCTTGATGCTTTTGTTCAAGCTTTTTGATCAATTTTTTATCGTCTTTTGTAATTTCGCGAAt                                                                                                        >  1:389328/1‑141 (MQ=255)
                                  tCAACGTCTGCCTTCGTAATGGCTTTTGTAGAAATTCTAACTAAATAATTTCGAATGCTATCATTGATTGTTTCAACAGCTTGATGCTTTTGTTCAAGCTTTTTGATCAATTTTTTATCGTCTTTTGTAATTTCGCGAATg                                                                                                       <  2:8168/141‑1 (MQ=255)
                                      cggcTGCCTTCGTAATGGCCATTGTATAAATTCTAACTAAATAGTTTCGAATGGTATCAGTCATCGTTTCAACACCTGGATGCTTTTGTTCAAGCTCTTCGAGCAATTTTTTATCGTCTTTTGTAATGTAGCTCATGTCtt                                                                                                   >  2:358374/4‑141 (MQ=255)
                                      cGTCTGCCTTCGTAATGGCTTTTGTAGAAATTCTAACTAAATAATTTCGAATGCTATCATTGATTGTTTCAACAGCTTGATGCTTTTGTTCAAGCTTTTTGATCAATTTTTTATCGTCTTTTGTAATTTCGCGAATGTCtt                                                                                                   >  1:351703/1‑141 (MQ=255)
                                      cGTCTGCCTTCGTAATGGCTTTTGTAGAAATTCTAACTAAATAATTTCGAATGCTATCATTGATTGTTTCAACAGCTTGATGCTTTTGTTCAAGCTTTTTGATCAATTTTTTATCGTCTTTTGTAATTTCGCGAATGTCtt                                                                                                   >  1:35209/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:80688/141‑1 (MQ=255)
                                                                            aaaTAATTTCGAATGCTATCATTGATTGTTTCAACAGCTTGATGCTTTTGTTCAAGCTTTTTGATCAATTTTTTATCGTCTTTTGTAATTTCGCGAATGTCTTCAAACATTGATAAGACAATCTGACCCACATTTTGTAAt                                                             <  2:312255/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)
                                                                                                                   tgaacaaaagcatCAAGCTTTTTGATCAATTTTTTATCGTCTTTTGTAATTTCGCGAATGTCTTCAAACATTGATAAGACAATCTGACCCACATTTTGTAATTCTTTTTGAGTTTCTTGTAATGCAACACCAGGTGCGTGa                      >  1:358833/13‑141 (MQ=255)
                                                                                                                     aTGCTTTTGTTCAAGCTTTTTGATCAATTTTTTATCGTCTTTTGTAATTTCGCGAATGTCTTCAAACATTGATAAGACAATCTGACCCACATTTTGTAATTCTTTTTGAGTTTCTTGTAATGCAACACCAGGTGCGTGATa                    >  2:387361/1‑141 (MQ=255)
                                                                                                                        cTTTTGTTCAAGCTTTTTGATCAATTTTTTATCGTCTTTTGTAATTTCGCGAATGTCTTCAAACATTGATAAGACAATCTGACCCACATTTTGTAATTCTTTTTGAGTTTCTTGTAATGCAACACCAGGTGCGTGATAAAc                 <  1:358374/141‑1 (MQ=255)
                                                                                                                          tttGTTCAAGCTTTTTGATCAATTTTTTATCGTCTTTTGTAATTTCGCGAATGTCTTCAAACATTGATAAGACAATCTGACCCACATTTTGTAATTCTTTTTGAGTTTCTTGTAATGCAACACCAGGTGCGTGATAAAC‑aa              >  2:328307/1‑141 (MQ=255)
                                                                                                                              ttCAAGCTTTTTTATTAATTTTTTATCGTCTTTTGTAATTTCGCGAATGTCTTCAAACATTGATAAGACAATCTGACCCACATTTTGTAATTCTTTTTGAGTTTCTTGTAATGCAACACCAGGTGCGTGATAAAC‑AAGATc          >  2:276151/1‑141 (MQ=255)
                                                                                                                                 aaGCTTTTTGATCAATTTTTTATCGTCTTTTGTAATTTCGCGAATGTCTTCAAACATTGATAAGACACTCTGACCCACATTTTGTAATTCTTTTTGAGTTTCTTGTAATGCAACACCAGGTGCGTGATAAACAAAGATCtt        <  2:99813/141‑1 (MQ=255)
                                                                                                                                     tttttGATCAATTTTTTATCGTCTTTTGTACTTTCGCGAATGTCTTCAAACATTGATAAGCCAATCTGACCCACATTTTGTAATTCTTTTTGAGTTTCTTGTAATGCAACACCAGGTGCGTGATAAAC‑AAGATCTTTGttt   <  2:53316/141‑1 (MQ=255)
                                                                                                                                      ttttGATCAATTTTTTATCGTCTTTTGTAATTTCGCGAATGTCTTCAAACATTGATAAGACAATCTGACCCACATTTTGTAATTCTTTTTGAGTTTCTTGTAATGCAACACCAGGTGCGTGATAAAC‑AAGATCTTTGTTTa  <  1:181613/141‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                            
AGAGCGATTGACATCAAACATAACTGCTAAACGCTCAACGTCTGCCTTCGTAATGGCTTTTGTAGAAATTCTAACTAAATAATTTCGAATGCTATCATTGATTGTTTCAACAGCTTGATGCTTTTGTTCAAGCTTTTTGATCAATTTTTTATCGTCTTTTGTAATTTCGCGAATGTCTTCAAACATTGATAAGACAATCTGACCCACATTTTGTAATTCTTTTTGAGTTTCTTGTAATGCAACACCAGGTGCGTGATAAAC‑AAGATCTTTGTTTA  >  CP000730/119356‑119630

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