Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 155,546 A→C 13.8% intergenic (‑251/+32) 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,5460AC13.8% 111.7 / 14.6 52intergenic (‑251/+32)yjbB/ygdQNa+/Pi‑cotransporter/transport protein
Reads supporting (aligned to +/- strand):  ref base A (21/23);  new base C (4/3);  total (25/27)
Fisher's exact test for biased strand distribution p-value = 7.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.37e-01

GTCGCCAGCTTCGCATGGCAGTCATTTATTTGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTGCGCATGCGGATGTTGATGGCTTCCACCGCCAGCGAGAAGGCC  >  NC_002947/155400‑155692
                                                                                                                                                  |                                                                                                                                                  
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)
                                                                                                                                                  |                                                                                                                                                  
GTCGCCAGCTTCGCATGGCAGTCATTTATTTGAAACAGATCTTGAAACGAATGACTGCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTGCGCATGCGGATGTTGATGGCTTCCACCGCCAGCGAGAAGGCC  >  NC_002947/155400‑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: ‑