Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,340,191 G→A 7.6% intergenic (‑126/‑84) mlaF ← / → yrbG intermembrane phospholipid transport system, ATP binding subunit MlaF/putative transport protein YrbG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,340,1910GA 7.6% 134.8 / 3.1 53intergenic (‑126/‑84)mlaF/yrbGintermembrane phospholipid transport system, ATP binding subunit MlaF/putative transport protein YrbG
Reads supporting (aligned to +/- strand):  ref base G (26/23);  new base A (2/2);  total (28/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.69e-01

CCACAGACTGCTCCATAATTCACCCTTCGTCTTGCGTTGATTTTCTAAGCATGGCGCTCAATTTAACCTTGAACCCAACATATTTACAGAATATTACCCGCCGTGGTTAGCGAAAGCTGGCATTTGTTTTACTTTTTAGCCGCATAAAGTCAAAATTAAGCATCCGTTACGGCTTTCTGAAAATCTTCAGCGGACCGGCGAGTATACCTGAAGAAAGGACGTTAGATGCTTTTAGCTACGGCACTGTTAATTGTTGGTTTACTTTTGGTCGTTTACAGTGCCGA  >  NC_000913/3340050‑3340333
                                                                                                                                             |                                                                                                                                              
ccACAGACTGCTCCATAATTCACCCTTCGTCTTGCGTTGATTTTCTAAGCATGGCGCTCAATTTAACCTTGAACCCAACATATTTACAGAATATTACCCGCCGTGGTTAGCGAAAGCTGGCATTTGTTTTACTTTTTAGCCGc                                                                                                                                               >  1:609367/1‑143 (MQ=255)
  acaGACTGCTCCATAATTCACCCTTCGTCTTGCGTTGATTTTCTAAGCATGGCGCTCAATTTAACCTTGAACCCAACATATTTACAGAATATTACCCGCCGTGGTTAGCGAAAGCTGGCGTTTGTGTTACTTTGTAGCGGCAt                                                                                                                                             >  1:49026/1‑143 (MQ=255)
        tGCTCCATAATTCACCCTTCGTCTTGCGTTGATTTTCTAAGCATGGCGCTCAATTTAACCTTGAACCCAACCTATTTACAGAAGATTACCCGCCGTGGTGAGCGAAAGCTGGCACTTGTTTTACTTTTTAGGCGCATAAAggc                                                                                                                                       >  2:610627/1‑141 (MQ=255)
            ccATAATTCACCCTTCGTCTTGCGTTGATTTTCTAAGCATGGCGCTCAATTTAACCTTGAACCCAACATATTTACAGAATATTACCCGCCGTGGTTAGCGAAAGCTGGCATTTGTTTTACTTTTTAGCCGCATAAAGTCaaaa                                                                                                                                   <  2:535358/143‑1 (MQ=255)
                aaTTCACCCTTCGTCTTGCGTTGATTTTCTAAGCATGGCGCTCAATTTAACCTTGAACCCAACATATTTACAGAATATTACCCGCCGTGGTTAGCGAAAGCTGGCATTTGTTTTACTTTTTAGCCGCATAAAGTCAAAATTaa                                                                                                                               <  1:581576/143‑1 (MQ=255)
                             tCTTGCGTTGATTTTCTAAGCATGGCGCTCAATTTAACCTTGAACCCAACATATTTACAGAATATTGCCCGCCGTGGTTAGCGAAAGCTGGCATTTGTTTTACTTTTTAGCCGCATAAAGTCAAAATTAAGCATCCGTTACgg                                                                                                                  >  1:471575/1‑143 (MQ=255)
                             gcTTTCGTTGATTTTCTAAGCAGGGCGCTCAAGTTACCCTTGAATCCAACGTATTTACAGCATATTACCTGCCGTGGTTAGCGAAAGCTGGCATGTGTTTTACTTTTTAGCCGCATAAAGTCAAAATTAAGCATCCGTTACgg                                                                                                                  <  2:92100/142‑1 (MQ=255)
                                 gCGTTGATTTTCTAAGCATGGCGCTCAATTTAACCTTGAACCCAACATATTTACAGAATCTTACCCGCAGTGGTTAGCGAAAGCTGGCATTTGTTTTACTTTTTAGCCGCATAAAGTCAAAATTAAGCATCCGTTACGGCttt                                                                                                              >  1:561111/1‑143 (MQ=255)
                                   gTTGATTTTCTAATCCTGGCGCTCAATTTAATCTTGCACCCAACATATTTACAGACTATTACCCGCCGTGGTTAGCGCAAGCTGGCATTTGTTTTACTTTTTAGCCGCATAAAGTCAAAATTAAGCATCCGTTACGGCTTTCt                                                                                                            <  2:135362/143‑1 (MQ=255)
                                   gTTGATTTTCTAAGCATGGCGCTCAATTTAACCTTGAACCCAACATATTTACAGAATATTACCCGCCGTGGTTCGCGAAAGCTGGCCTTTGTTTTACTTTTTAGCCGCATAAAGTCACACTTACGCATCCGTTACGGCTTTCt                                                                                                            >  1:779973/1‑143 (MQ=255)
                                    ttGATTTTCTAAGCATGGCGCTCAATTTAACCTTGAACCCAACATATTTACAGAATATTACCCGCCGTGGTTAGCGAAAGCTGGCATTTGTTTTACTTTTTAGCCGCATAAAGTCAAAATTAAGCATCCGTTACGGCTTTCTg                                                                                                           >  2:18417/1‑143 (MQ=255)
                                      gATTTTCTAGGCATGGCGCTCAATTTAACCTTGAACCCAACATATTTACAGAATATTACCCGCCGTGGTTAGCGAAAGCTGGCATTTGTTTTACTTTTTAGCCGCATAAAGTCAAAATTAAGCATCCGTTACGGCTTTCTGaa                                                                                                         <  1:636511/143‑1 (MQ=255)
                                      gATTTTCTAAGCATGGCGCTCAATTTAACCTTGAACCCAACATATTTACAGAATATTACCCGCCGTGGTTAGCGAAAGCTGGCATTTGTTTTACTTTTTCGCCGCATAAAGTCAAAATTAAGCATCCGTTCCGGCTTTCTGaa                                                                                                         >  1:167349/1‑143 (MQ=255)
                                       aTTTTCTAAGCATGGCGCTCAATTTAACCTTGAACCCAACATATTTACAGAATATTACCCGCCGTGGTTAGCGAAAGCTGGCATTTGTTTTACTTTTTAGCCGCATAAAGTCAAAATTAAGCATCCGTTACGGCTTTCTGaaa                                                                                                        >  1:273110/1‑143 (MQ=255)
                                         tttCTAAGCATGGCGCTCAATTTAACCTTGAACCCAACATATTTACAGAATATTACCCGCCGTGGTTAGCGAAAGCTGGCATTTGTTTTACTTTTTAGCCACATAAAGTCAAAATTAAGCATCCGTTACGGCTTTCTGAAAAt                                                                                                      >  1:598768/1‑143 (MQ=255)
                                         tttCTAAGCATGGCGCTCAATTTAACCTTGAACCCAACATATTTACAGAATATTACCCGCCGTGGTTAGCGAAAGCTGGCATTTGTTTTACTTTTTAGCCACATAAAGTCAAAATTAAGCATCCGTTACGGCTTTCTGAAAAt                                                                                                      >  1:598747/1‑143 (MQ=255)
                                           tCTAAGCATGGCGCTCAATTTAACCTTGAACCCAACATATTTACAGAATATTACCCGCCGTGGTTAGCGAAAGCTGGCATTTGTTTTACTTTTTAGCCGCATAAAGTCAAAATTAAGCATCCGTTACGGCTTTCTGAAAATCt                                                                                                    >  1:126798/1‑143 (MQ=255)
                                              atgCATGGCGTTCCATTTAACCGTGGACCCACCATCTTACCAGAATATTACCCGCCGTGGTTAGCGAAAGCTGGCATTTGTTTTACTTTTTAGCCGCATAAAGTCAAAATTAAGCATCCGTTACGGCTTTCTGAAAATCTTCa                                                                                                 <  1:517661/141‑1 (MQ=255)
                                                      cgcTCAATTTAACCTTGAACCCAACATATTTACAGAATATTACCCGCCGTGGTTAGCGAAAGCTGGCATTTGTTTTCCTTTTTAGCCGCATAAAGTCAAAATTAAGCATCCGTTACGGCTTTCTGAAAATCTTCAGCGGACCg                                                                                         >  2:437765/1‑143 (MQ=255)
                                                        cTCAATTTAACCTTGAACCCAACATATTTACAGAATATTACCCGCCGTGGTTAGCGAAAGCTGGCATTTGTTTTACTTTTTAGCCGCATAAAGTCAAAATTAAGCATCCGTTACGGCTTTCTGAAAATCTTCAGCGGACCGGc                                                                                       >  1:310227/1‑143 (MQ=255)
                                                            aTTTAACCTTGAACCCAACATATTTACAGAATATTACCCGCCGTGGTTAGCGAAAGCTGGCATTTGTTTTACTTTTTAGCCGCATAAAGTCAAAATTAAGCATCCGTTACGGCTTTCTGAAAATCTTCAGCGGACCGGCGAGt                                                                                   <  1:37447/143‑1 (MQ=255)
                                                                aaCCTTGAACCCAACTTATTTACAGAATATTACCCGCCGTGGTTAGCGAAAGCTGGCATTTGTTTTACTTTTTAGCCGCATAAAGTCAAAATTAAGCATCCGTTACGGCTTTCTGAAAATCTTCAGCGGACCGGCGAGTATAc                                                                               <  2:668882/143‑1 (MQ=255)
                                                                    ttGAACCCAACATATTTACAGAATATTACCCGCCGTGGTTAGCGAAAGCTGGCATTTGTTTTACTTTTTAGCCGCATAAAGTCAAAATTAAGCATCCGTTACGGCTTTCTGAAAATCTTTAGCCGACCGGCGAGTATACCTga                                                                           >  2:124958/1‑143 (MQ=255)
                                                                    ttGAACCCAACATATTTACAGAATATTACCCGCCGTGGTTAGCGAAAGCTGGCATTTGTTTTACTTTTTAGCCGCATAAAGTCAAAATTAAGCATCCGTTACGGCTTTCTGAAAATCTTCAGCGGACCGGCGAGTATACCTga                                                                           >  1:100716/1‑143 (MQ=255)
                                                                    ttGAACCCAACATATTTACAGAATATTACCCGCCGTGGTTAGCGAAAGCTGGCATTTGTTTTACTTTTTAGCCGCATAAAGTCAAAATTAAGCATCCGTTACGGCTTTCTGAAAATCTTCAGCGGACCGGCGAGTATACCTga                                                                           >  2:697833/1‑143 (MQ=255)
                                                                         cccAACATGTTTACAGAATAGTACCCGCCGTGGTTAGCGAAAGCTGGCATTTGTTTTACTTTTTAGCCGCATAAAGTCAAAATTAAGCATCCGTTACGGCTTTCTGAAAATCTTCAGCGGACCGGCGAGTATACCTGAAGaaa                                                                      <  1:694014/143‑1 (MQ=255)
                                                                          ccAACATATTTACAGAATATTACCCGCCGTGGTTAGCGAAAGCTGGCATTTGTTTTACTTTTTAGCCGCATAAAGTCAAAATTaa                                                                                                                               <  2:284333/85‑1 (MQ=255)
                                                                          ccAACATATTTACAGAATATTACCCGCCGTGGTTAGCGAAAGCTGGCATTTGTTTTACTTTTTAGCCGCATAAAGTCAAAATTaa                                                                                                                               >  1:284333/1‑85 (MQ=255)
                                                                            aaCATTTTGTCAGAATATTGCCCGCCGTGGGTAGCGAAAGCTGGCATTTGTTTTACTTTTTAGCCGCAGAAAGTCAAAATTAAGCATCCGTTACGGCTTTCGGAAAATCTTCAGCGGACCGGCGAGTATACCTGAAGAAAGGa                                                                   <  2:714077/143‑1 (MQ=255)
                                                                               atatTTACAGAATATTACCCGCCGTGGTTAGCGAAAGCTGGCATTTGTTTTACTTTTTAGCCACATAAAGTCAAAATTAAGCATCCGTTACGGCTTTCTGAAAATCTTCAGCGGACCGGCGAGTATACCTGAAGAAAGGACGt                                                                <  2:598747/143‑1 (MQ=255)
                                                                               atatTTACAGAATATTACCCGCCGTGGTTAGCGAAAGCTGGCATTTGTTTTACTTTTTAGCCACATAAAGTCAAAATTAAGCATCCGTTACGGCTTTCTGAAAATCTTCAGCGGACCGGCGAGTATACCTGAAGAAAGGACGt                                                                <  2:598768/143‑1 (MQ=255)
                                                                                tatCTACAGAATATTACCCGCCGTGGTTAGCGCAAGCTAGCATCTGTTTTACTTTTTAGCCGCATAAAGTCAAAATTAAGCATTCGTTACGGCTTTCTGAAAATCTTCAGCGGACCGGCGAGTATACCTGAAGAAAGGACGtt                                                               <  1:742315/143‑1 (MQ=255)
                                                                                  tttACAGAATATTACCCGCCGTGGTTAGCGAAAGCTGGCATTTGTTTTACTTTTTAGCCGCATAAAGTCAAAATTAAGCATCCGTTACGGCTTTCTGAAAATCTTCAGCGGACCGGCGAGTATACCTGAAGAAAGGACGCTAg                                                             >  2:182831/1‑143 (MQ=255)
                                                                                   ttcctgAATATGACCCTCCCTGCGTAGCCAAAGTCGGCATTTGTTTTACTTTTTAGCCGCATAAAGTCAAAATTAAGCATCCGTTACGGCTTTCTTAAAATCTTCAGCGGACCGGCGAGTATACCTGAAGAAAGGACGTTAGa                                                            <  1:520813/138‑1 (MQ=255)
                                                                                    tACAGAATATTACCCGCCGTGGTTAGCGAAAGCTGGCATTTGTTTTACTTTTTAGCCGCATAAAGTCAAAATTAAGCATCCGTTACGGCTTTCTGAAAATCTTCAGCGGACCGGCGAGTATACCTGAAGAAAGGACGTTAGAt                                                           <  2:273110/143‑1 (MQ=255)
                                                                                      cAGAATATTACCCGCCGTGGTTAGCGAAAGCTGGCATTTGTTTTACTTTTTAGCCGCATAAAGTCAAAATTAAGCATCCGTTACGGCTTTCTGAAAATCTTCAGCGGACCGGCGAGTATACCTGAAGAAAGGACGTTAGATGc                                                         >  1:598120/1‑143 (MQ=255)
                                                                                        gAATATTACCCGCCGTGGTTAGCGAAAGCTGGCATTTGTTTTACTTTTTAGCAGCATCATGTCAAAATTAAGCATCCGTTACGGCTTTCTGCAACTCTTCAGTCGGCCGACGATTTTACCTGAAGCAAGGACGTTGCATGCtt                                                       >  2:310822/1‑143 (MQ=255)
                                                                                          atatTACCCGCCGTGGTTAGCGAAAGCTGGCATTTGTTTTACTTTTTAGCCGCATAAAGTCAAAATTAAGCATCCGTTACGGCTTTCTGAAAATCTTCAGCGGACCGGCGAGTATACCTGAAGAAAGGACGTTAGATGCtttt                                                     <  1:18417/143‑1 (MQ=255)
                                                                                               aCCCGCCGTGGTGAGCGAAAGCTGGCATTTGTTTGACTTGTTAGCCGCATAAAGTCAAAATTAAGCATCCGTTACGGCTTTCTGAAAATCTTCAGCGGACCGGCGAGTATACCTGAAGAAAGGACGTTAGATGCTTTTAGCTa                                                <  2:782542/143‑1 (MQ=255)
                                                                                                     cgTGGTTAGCGAAAGCTGGCATTTGTTTTACTTTTTAGCCGCATAAAGTCAAAATTAAGCATCCGTTACGGCTTTCTGAAAATCTTCAGCGGACCGGCGAGTATACCTGAAGAAAGGACGTTAGATGCTTTGAGCTACGGCAc                                          >  2:783700/1‑143 (MQ=255)
                                                                                                      gTGGTTAGCGTAAGCTGGCATGTGTTTTACTTTTTAGCCGAATATTGTTAGTGTTGGTGGTCTGTCGGGCGGTGTTGAAAATCTTCAGCCGACTGTCGAGTATACCTGAAGAAAGTACGTTAGATGCTTTTAGCTACGGCACt                                         <  2:165752/143‑1 (MQ=255)
                                                                                                           tAGCGAAAGCTGGCATTTGTTTTACTTTTTAGCCGCATAAAGTCAAAATTAAGCAGCCGTTACGGCTTTCTGAAAATCTTCAGCGGACCGGCGAGTATACCTGAAGAAAGGACGTTAGATGCTTTTAGCTACGGCACTGTTaa                                    <  1:663580/143‑1 (MQ=255)
                                                                                                            aGCGAAAGCTGGCATTTGTTTTACTTTTTAGCGGCATAAAGTCAAAATTAAGCATCCGTTATGGCTTGCTGAAAATCTTCAGCGGAGCGGCGAGTATACCTGAAGAAAGGACGTTAGATGCTTTTAGCGACGGCACTGGTAAt                                   >  2:230212/1‑143 (MQ=255)
                                                                                                              cGAAAGCTGGCATTTGTTTTACTTTTTATCCGCCTAAAGTCAAAATTAAGCATCCGTTACGGCTTTCGGAAAATCTTCAGCGGACCGGCGAGTATACCTGAAGAGAGGAAGTACCATGCTTTTAGCTACAGCACTGTTaaatc                                 >  2:627133/1‑140 (MQ=255)
                                                                                                                         aTTTGTTTTACTTTTTAGCCGCATAAAGTCAAAATTAAGCATCCGTTACGACTTTCTGAAAATCTTAAGCGGACGGGAGGGTATAGGTGAAGAAAGGAGGTTTGATGCTGTTATCTACGGCACTGTTAATTGTTGGTTTGCtt                      >  1:782127/1‑143 (MQ=255)
                                                                                                                          tttGTTTTACTTTTTAGCCGCATAAAGTCAAAATTAAGCATCCGTTACGGCTTTCTGAAAATCTTCAGCGGACCGGCGAGTATACCTGAAGAAAGGACGTTAGATGCTTTTAGCTACGGCACTGTTAATTGTTGGTTTACttt                     <  1:640206/143‑1 (MQ=255)
                                                                                                                             tttttACTTTTTAGCCGCATAAAGTCAAAATTAAGCATCCGTTACGGCTTTCTGAAAATCTTCAGCGGACCGGCGAGTATACCTGAAGAAAGGACGTTAGATGCTTTTAGCTACGGCACTGTTAATTGTTGGTTTACTTTTgg                  <  2:95006/142‑1 (MQ=255)
                                                                                                                              ttttACTTTTTAGCCGCATAAAGTCAAAATTAAGCATCCGTTACGGCTTTCTGAAAATCTTCAGCGGACCGGCGAGTATACCTGAAGAAAGGACGTTAGATGCTTTTAGCTACGGCACTGTTAATTGTTGGTTTACTTTTGGt                 <  2:100716/143‑1 (MQ=255)
                                                                                                                              ttttACTTTTTAGCCGCATAAAGTCAAAATTAAGCATCCGTTACGGCTTTCTGAAAATCTTCAGCGGACCGGCGAGTATACCTGAAGAAAGGACGTTAGATGCTTTTAGCTACGGCACTGTTAATTGTTGGTTTACTTTTGGt                 >  1:601276/1‑143 (MQ=255)
                                                                                                                                    tttttAGCCGCATAAAGTCAAAATTAAGCATCCGTTACGGCTTTCTGAAAATCTTCAGCGGACCGGCGAGTATACCTGAAGAAAGGACGTTAGATGCTTTTAGCTACGGCACTGTTAATTGTTGGTTTACTTTTGGTCGTTTa           >  2:42132/1‑143 (MQ=255)
                                                                                                                                      tttAGCCGCATAAAGTCAAAATTAAGCATCCGTTACGGCTTTCTGAAAATCTTCAGCGGACCGGCGAGTATACCTGAAGAAAGGACGTTAGATGCTTTTAGCTACGGCACTGTTAATTGTTGGTTTACTTTTGGTCGTTTACa         <  1:364365/143‑1 (MQ=255)
                                                                                                                                           ccGCATAAAGTCACAATTCAGCATCCGTTACGGCTTTCTGAAAATCTTCAGCGGACCGGCGAGTATACCTGAAGAAAGGACGTTAGATGCTTTTAGCTACGGCACTGTTAATTGTTGGTTTACTTTTGGTCGTTTACAGTGcc    <  2:328584/143‑1 (MQ=255)
                                                                                                                                           ccGCATAAAGTCAAAATTAAGCATCCGTTACGGCTTTCTGAAAATCTTCAGCGGACCGGCGAGTATACCTGAAGAAAGGACGTTAGATGCTTTTAGCTACGGCACTGTTAATTGTTGGTTTACTTTTGGTCGTTTACAGTGcc    >  1:723248/1‑143 (MQ=255)
                                                                                                                                             gCATAAAGTCAAAATTAAGCATCCGTTACGGCTTTCTGAAAATCTTCAGCGGACCGGCGAGTATACCTGAAGAAAGGACGTTAGATGCTTTTAGCTACGGCACTGTTAATTGTTGGTTTACTTTTGGTCGTTTACAGTGCCGa  <  2:206359/143‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                              
CCACAGACTGCTCCATAATTCACCCTTCGTCTTGCGTTGATTTTCTAAGCATGGCGCTCAATTTAACCTTGAACCCAACATATTTACAGAATATTACCCGCCGTGGTTAGCGAAAGCTGGCATTTGTTTTACTTTTTAGCCGCATAAAGTCAAAATTAAGCATCCGTTACGGCTTTCTGAAAATCTTCAGCGGACCGGCGAGTATACCTGAAGAAAGGACGTTAGATGCTTTTAGCTACGGCACTGTTAATTGTTGGTTTACTTTTGGTCGTTTACAGTGCCGA  >  NC_000913/3340050‑3340333

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