Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 155,543 G→T 13.9% intergenic (‑248/+35) yjbB ← / ← ygdQ Na+/Pi‑cotransporter/transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947155,5430GT13.9% 131.4 / 9.9 51intergenic (‑248/+35)yjbB/ygdQNa+/Pi‑cotransporter/transport protein
Reads supporting (aligned to +/- strand):  ref base G (20/23);  new base T (4/3);  total (25/26)
Fisher's exact test for biased strand distribution p-value = 6.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.69e-01

CCAGCTGTCGCCAGCTTCGCATGGCAGTCATTTATTTGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTGCGCATGCGGATGTTGATGGCTTCCACCGCCAGCGAGAAGGCC  >  NC_002947/155394‑155692
                                                                                                                                                     |                                                                                                                                                     
ccAGCTGTCGCCAGCTTCGCATGGCAGTCATTTATTTGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTg                                                                                                                                                       <  1:92086/150‑1 (MQ=255)
      gTCGCCAGCTTCGCATGGCAGTCATTTATTTGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAG‑AAAAAGGCATTAAAAAGGGCCCTGCGAAggg                                                                                                                                                <  2:732680/150‑1 (MQ=255)
       tCGCCAGCTTCGCATGGCAGTCATTTATTTGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAggg                                                                                                                                                <  2:436348/150‑1 (MQ=255)
       tCGCCAGCTTCGCATGGCAGTCATTTATTTGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCACTGCGAAggg                                                                                                                                                >  1:598893/1‑150 (MQ=255)
          ccAGCTTCGCATGGCAGTCATTTATTTGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGccc                                                                                                                                             >  2:429984/1‑150 (MQ=255)
             gCTTCGCATGGCAGTCATTTATTTGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTTCGCAGGGCCCttt                                                                                                                                          <  2:892187/150‑2 (MQ=255)
              cTTCGCATGGCAGTCATTTATTTGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCttgt                                                                                                                                         <  1:606372/150‑1 (MQ=255)
               ttCGCATGGCAGTCATTTATTTGAAACAGAACTTGAAACGAATGACTGCCCGTAAGGGCAGCTTCGCGGGCAAGCAATTACAGTGAAAAAAAAAAAAAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCttgtt                                                                                                                                        <  2:271005/150‑1 (MQ=255)
                tCGCATGGCAGTCATTTATTTGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTTCGCAGGGCCCttttta                                                                                                                                       <  1:777677/150‑2 (MQ=255)
                  gCATGGCAGTCATTTATTTGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTtgtg                                                                                                                                     <  2:475167/150‑1 (MQ=255)
                     tGGCAGTCATTTATTTGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGtgctg                                                                                                                                  >  1:1197166/1‑150 (MQ=255)
                           tCATTTATTTGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTc                                                                                                                            <  1:960564/150‑1 (MQ=255)
                               ttatttGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTg                                                                                                                        <  1:479092/150‑1 (MQ=255)
                                     gAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCGTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTa                                                                                                                  >  1:1151909/1‑150 (MQ=255)
                                      aaaCAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTAt                                                                                                                 >  1:902765/1‑150 (MQ=255)
                                       aaCAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTAtt                                                                                                                >  2:145789/1‑150 (MQ=255)
                                         cAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCAGCACAGTGACAACGCGCAGGAAAAAACGCACAGAAGAAAGGCATTAAAAAGCGCCCTCCGAACGGCCCTTGTTGTGCTGCACTCCCTTGGCCTTATTGa                                                                                                              >  2:157351/1‑150 (MQ=255)
                                            aTCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGAccc                                                                                                           >  1:810177/1‑150 (MQ=255)
                                               ttGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGa                                                                                                        >  2:483034/1‑150 (MQ=255)
                                                 gAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAAt                                                                                                      >  2:80299/1‑150 (MQ=255)
                                                          gACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTg                                                                                                                        >  1:398730/1‑123 (MQ=255)
                                                          gACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTg                                                                                                                        <  2:398730/123‑1 (MQ=255)
                                                            cTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCa                                                                                           <  2:1151909/150‑1 (MQ=255)
                                                             tGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTAtt                                                                                                                <  1:260406/128‑1 (MQ=255)
                                                             tGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTAtt                                                                                                                >  2:260406/1‑128 (MQ=255)
                                                                cagTGGGGACTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGttt                                                                                       <  1:807200/148‑1 (MQ=255)
                                                                  gTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCa                                                                                     <  1:483034/150‑1 (MQ=255)
                                                                   tGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCAc                                                                                    <  1:485441/150‑1 (MQ=255)
                                                                        gcgtctTCGGGGGCAAGCCCCCACAGTGACAACGCGCAGCAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCCCTTCCTTGGCCTTATTGACCCGGAATGTCCTTGATCAGTTTCACTGGtt                                                                               <  2:598893/146‑1 (MQ=255)
                                                                                  gggCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAATGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCtt                                                                     <  2:810177/150‑1 (MQ=255)
                                                                                      aaGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCa                                                                 >  2:732713/1‑150 (MQ=255)
                                                                                        gCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTAt                                                                                                                 >  1:954241/1‑100 (MQ=255)
                                                                                        gCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTAt                                                                                                                 <  2:954241/100‑1 (MQ=255)
                                                                                             cacaGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGc                                                          <  1:732713/150‑1 (MQ=255)
                                                                                              acaGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCTTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCg                                                         <  1:429984/150‑1 (MQ=255)
                                                                                                          cgcAGGAAAAAACGCACAGAAAAAAGGCATTCAACAGGGCCCAGCGGGTTGTCCTTGTTGTGCGGCACTTCCCTAGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACAGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGt                                             <  1:80299/150‑1 (MQ=255)
                                                                                                            cAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTgc                                           <  2:1197166/150‑1 (MQ=255)
                                                                                                             aGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCTCAGTTTCCCCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTgcg                                          >  1:670902/1‑150 (MQ=255)
                                                                                                                         acaGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTGCGCATGCGGATGtt                              >  1:597577/1‑150 (MQ=255)
                                                                                                                           aGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTGCGCATGCGGATGTTGa                            >  2:421228/1‑150 (MQ=255)
                                                                                                                             aaaaaaGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTGCGCATGCGGATGTAGATg                          >  2:637075/1‑150 (MQ=255)
                                                                                                                                      aTTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGa                                                                          >  2:964223/1‑93 (MQ=255)
                                                                                                                                      aTTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGa                                                                          <  1:964223/93‑1 (MQ=255)
                                                                                                                                       acaacaaGGGCCCTTCGCAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTGCGCATGCGGATGTTGATGGCTTCCACCg                >  2:747098/6‑150 (MQ=255)
                                                                                                                                          acaaGGGCCCTTCGCAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTGCGCATGCGGATGTTGATGGCTTCCACCGCCa             >  2:525739/3‑150 (MQ=255)
                                                                                                                                          acaaGGGCCCTTCGCAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTGCGCATGCGGATGTTGATGGCTTCCACCGCCa             >  2:525250/3‑150 (MQ=255)
                                                                                                                                          aaaaGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTGCGCATGCGGATGTTGATGGCTTCCACCGcct             <  2:902765/150‑2 (MQ=255)
                                                                                                                                             aGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTGCGCATGCGGATGTTGATGGCTTCCACCGCCAGCg          >  1:82444/1‑150 (MQ=255)
                                                                                                                                              gggCCCTTCGCAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTGCGCATGCGGATGTTGATGGCTTCCACCGCCAGCga         <  1:747098/150‑1 (MQ=255)
                                                                                                                                              gggCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCTGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTGCGCATGCGGATGTTGATGGCTTCCACCGCCAGCga         <  1:145789/150‑1 (MQ=255)
                                                                                                                                                gCCCTTCGCAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTGCGCATGCGGATGTTGATGGCTTCCACCGCCAGCgaga       >  2:904522/1‑150 (MQ=255)
                                                                                                                                                    tGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTGCGCATGCGGATGTTGATGGCTTCCACCGCCAGCGAGAAGGc   >  2:589967/1‑150 (MQ=255)
                                                                                                                                                     gCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTGCGCATGCGGATGTTGATGGCTTCCACCGCCAGCGAGAAGGcc  <  1:1185593/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCAGCTGTCGCCAGCTTCGCATGGCAGTCATTTATTTGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTGCGCATGCGGATGTTGATGGCTTCCACCGCCAGCGAGAAGGCC  >  NC_002947/155394‑155692

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