Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,598,439 T→C 12.2% intergenic (+73/+115) panM → / ← livJ PanD autocleavage accelerator, panothenate synthesis/branched‑chain amino acid ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,598,4390TC12.2% 105.1 / 7.1 49intergenic (+73/+115)panM/livJPanD autocleavage accelerator, panothenate synthesis/branched‑chain amino acid ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base T (26/17);  new base C (3/3);  total (29/20)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.82e-01

TCGCGGTGTGATGACGGCGTTTATGCAGGCGCTGGGGTTTACGGCACAACAGGGCGGCTGGGAGAAGTGTTAATCGTCAAGTTTGATTTCAAAAGTGATATTGCCTGATGCGCTACGCTTATCAGGCCTACAATGTGTGTTGCAATTTACTGATTTCTTTGGATCTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAAGCAACGTACTCGATATTAGCAATTTGGCGGCAACCCAAAGTTGCCGATTAATGATTACTTCGCATCGGTCGCCGTGCCGTTGGCGTGC  >  NC_000913/3598294‑3598587
                                                                                                                                                 |                                                                                                                                                    
tCGCGGTGTGATGACGGCGTTTATGCAGGCGCTGGGGTTTACGGCACAACAGGGCGGCTGGGAGAAGTGTTAATCGTCAAGTTTGATTTCAAAAGTGATATTGCCTGATGCGCTACGCTTATCAGGCCTACAATGTGTGTTGCAATTTa                                                                                                                                                   >  2:486288/1‑149 (MQ=255)
tCGCGGTGTGATGACGGCGTTTATGCAGGCGCTGGGGTTTACGGCACAACAGGGCGGCTGGGAGAAGTGTTAATCGTCAAGTTTGATTTCAAAAGTGATATTGCCTGATGCGCTACGATTATCAGGCCTACAATGTGTGTTGCGATTTa                                                                                                                                                   >  2:219027/1‑149 (MQ=255)
tCGCGGTGTGATGACGGCGGTTATGCAGGCGCTGGGGTTTACGGCACAACAGGGCGGCTGGGAGAAGTGTTAATCGTCAAGTTTGATTTCAAAAGTGATATTGCCTGATGCGCTACGCTTATCAGGCCTACAATGTGTGTTGCAATTTa                                                                                                                                                   <  2:602945/149‑1 (MQ=255)
          atgaCGGCGTTTATGCAGGCGCTGGGGTTTACGGCACAACAGGGCGGCTGGGAGAAGTGTTAATCGTCAAGTTTGATTTCAAAAGTGATATTGCCTGATGCGCTACGCTTATCAGGCCTACAATGTGTGTTGCAATTTACTGATTTCtt                                                                                                                                         <  1:354431/149‑1 (MQ=255)
                     tATGCAGGCGCTGGGGTTTACGGCACAACAGGGCGGCTGGGAGAAGTGTTAATCGTCAAGTTTGATTTCAAAAGTGATATTGCCTGATGCGCTACGCTTATCAGGCCTACAATGTGTGTTGCAATTTACTGATTTCTTTGGATCTTGTa                                                                                                                              <  2:578619/149‑1 (MQ=255)
                         cAGGGGCTGGGGTTTACGGCACAACAGGGCGGCTGGGAGAAGTGTTAATCGTCAAGTTTGATTTCAAAAGTGATATTGCCTGATGCGCTACGCTTATCAGGCCTACAATGTGTGTTGCAATTTACTGATTTCTTTGGATCTTGTAGGcc                                                                                                                          <  1:207007/149‑1 (MQ=255)
                          aGGCGCTGGGGTTTACGGCACAACAGGGCGGCTGGGAGAAGTGTTAATCGTCAAGTTTGATTTCAAAAGTGATATTGCCTGATGCGCTACGCTTATCAGGCCTACAATGTGTGTTGCAATTTACTGATTTCTTTGGTTCTTGTAGGCCg                                                                                                                         >  2:40459/1‑149 (MQ=255)
                             cgcTGGGGTTTACGGCACAACAGGGCGGCTGGGAGAAGTGTTAATCGTCAAGTTTGATTTCAAAAGTGATATTGCCTGATGCGCTACGCTTATCAGGCCTACAATGTGTGTTGCAATTTACTGATTTCTTTGGATCTTGTAGGCCGGAt                                                                                                                      >  2:705634/1‑149 (MQ=255)
                                tGGGGTTTACGGCACAACAGGGCGGCTGGGAGAAGTGTTAATGGTCAAGTTTGATTTCAAAAGTGATATTGCCTGATGCGCTACGCTTATCAGGCCTACAATGTGTGTTGCAATTTACTGATTTCTTTGGATCTTGTAGGCCGGATAAg                                                                                                                   >  2:472546/1‑149 (MQ=255)
                                tGGGGTTTACGGCACAACAGGGCGGCTGGGAGAAGTGTTAATCGTCAAGTTTGATTTCAAAAGTGATATTGCCTGATGCGCTACGCTTATCAGGCCTACGATGTGTGTTGCAATTTACTGATTTCTTTGGATCTTGTAGGCCGGCTAAg                                                                                                                   >  1:88406/1‑149 (MQ=255)
                                tGGGGTTTACGGCACAACAGGGCGGCTGGGAGAAGTGTTAATCGTCAAGTTTGATTTCAAAAGTGATATTGCCTGATGCGCTACGCTTATCAGGCCTACGATGTGTGTTGCAATTTACTGATTTCTTTGGATCTTGTAGGCCGGATAAg                                                                                                                   >  1:99795/1‑149 (MQ=255)
                                  gggTTTGCGGCACGACGGGGCGGCTGGGAGAAGTGTTAATCGTCAAGTTTGATTTCAAAAGTGATATGGCCTGATGCGCTACGCTTATCAGGCCTACAATGTGTGTTGCAATTTACTGATTTCTTTGGATCTTGTAGGCCGGATAAGGc                                                                                                                 <  2:500058/149‑1 (MQ=255)
                                    gTTTACGGCACAACAGGGCGGCTGGGAAAAGTGTTAATCGTCAAGTTTGATTTCAAAAGTGATATTGCCTGATGCGCTAAGATTATCAGGCCTGCAATGTGTGTTGCAATGTACTGGTTTAGATCGATCTTGTAAGCTGGATAAGGCgg                                                                                                               >  2:5623/1‑148 (MQ=255)
                                        aCGGCACAACAGGGCGGGTGGGAGAAGTGTTAATCGTCAAGTTTGATTTCAAAAGTGATATTGCCTGATGCGCTACGCTTATCAGGCCTACAATGTGTGTTGCAATTTACTGATTTCTTTGGATCTTGTAGGCCGGATAAGGCGTTTAc                                                                                                           <  1:540188/149‑1 (MQ=255)
                                          ggCACAACAGGGGGGCTGGGAGAAGTGTGAATCGTCAAGTTTGATTTCAAAAGTGATATTGCCTGATGCGCTACGCTTATCAGGCCTACAATGTGTGTTGCAATTTACTGATTTCTTTGGATCTTGTAGGCCGGATAAGGCGTTTAcgc                                                                                                         <  1:548828/149‑1 (MQ=255)
                                           gCACAACAGGGCGGCTGGGAGAAGGTGTAGTCTTCAAGTGTGATTTCAAAAGTGATATTGCCTGATGCGCTACGCTTATCAGGTCTACAATGTGTGTTGCAATTTACTGATTTCTTTGGATCTTGTAGGCCGGATAAGGCGTTTAcgcc                                                                                                        <  1:390446/149‑1 (MQ=38)
                                                                   tgtTAATCGTCAAGTTTGATTTCAAAAGTGATATTGCCTGATGCGCTACGCTTATCAGGCCTACAATGTGTGTTGCAGTTTACTGATTTCTTTGGATCTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCCTGAAGACACGTACt                                                                                >  2:454731/1‑149 (MQ=255)
                                                                        aTCGTCAAGTTTGATTTCAAAAGTGATATTGCCTGATGCGCTACGCTTATCAGGCCTACAATGTGTGTTGCAACTTACTGATTTCTTTGGATCTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAAGCAACGTACTCGata                                                                           >  2:584380/1‑149 (MQ=255)
                                                                        aTCGTCAAGTTTGATTTCAAAAGTGATATTGCCTGATGCGCTACGCTTATCAGGCCTACAATGTGTGTTGCAACTTACTGATTTCTTTGGATCTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAAGCAACGTACTCGata                                                                           >  2:701151/1‑149 (MQ=255)
                                                                        aTCGTCAAGTTTGATTTCAAAAGTGATATTGCCTGATGCGCTACGCTTATCAGGCCTACAATGTGTGTTGCAACTTACTGATTTCTTTGGATCTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAAGCAACGTACTCGata                                                                           >  2:124822/1‑149 (MQ=255)
                                                                                 tttGATTTCAAAAGTGATATTGCCTGATGCGCTACGCTTATCAGGCCTACAATGTGTGTTGCAATTTACTGATTTCTTTGGATCTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAAGCAACGTACTCGATATTAGCAAtt                                                                  >  2:221757/1‑149 (MQ=255)
                                                                                 tttGATTTCAAAAGTGATATTGCCTGATGCGCTACGCTTATCAGGCCGACAATGTGTGTTGCAATTTACTGATTTCTTTGGATCTTGTAGGCCGGATAAGGCGTTTACGACGCATCCGGCATGAAGCCACGTAAGCGATATTAGCAAtt                                                                  >  2:45308/1‑149 (MQ=38)
                                                                                      tttCAAAAGTGATATTGCCTGATGCGCTACGCTGATCAGGCCTACAATGTGTTTTGCAATTTATTGATTTCTTTGGATCTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCTGCGTGAAGCAACGTACTCGATATTAGCAATTTggcg                                                             <  2:608188/149‑1 (MQ=37)
                                                                                        tCAAAAGTGATATTGCCTTGTGCGCTACGCTTATCAGCCCTACCGCGTTGGGTGCCCTTTAATGATTTCTTTGGATCTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAAGCAACGTACCCGAGATTAGCAATTTggcggc                                                           >  1:574239/1‑149 (MQ=25)
                                                                                                     tGCCTGATGCGCTCCGCTTATCAGGCCTACACTGTGTGTTGCAATTTACTGACCTCCTCGGACCTCCTACGCCGGCTAAGGCGTCTCCGCCGCCCCCGCCATGAAGCAACGTACCCGCCACTAGCAATCTGCCCGCCACCCCAAGtccc                                              >  2:509839/1‑146 (MQ=255)
                                                                                                        cTGATGGGCCACGCTTTTCAGGCCTACAATGTGTGTTGCAATCTACTGATTTCTTTGGATCTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAAGCAACGTACTCGATATTACCAATTTGGCGGCAACCCAAAGTTGCCGa                                           <  1:419395/149‑1 (MQ=38)
                                                                                                               gcTACGCTTATCTGGCCTACAATGTGGGTTGGAATTTACTGATTTCTTTGGAACATGTAAACCGTATAAGGCGTTTACGCCGCATCCGGCATGAAGCAACGTACTCGATATTAGCAATTTGGCGGCAACCCAAAGTTGCCGATTAATGa                                    <  1:182486/149‑1 (MQ=255)
                                                                                                               gcTACGCTTATCAGGCCTACAATGTGTGTTGCAATTTACTGATTTCTTTGGATCTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAAGCAACGTACTCGATATTAGCAATTTGGCGGCACCCCAAAGTTGCCGATCAATgc                                    >  1:447079/1‑148 (MQ=255)
                                                                                                                 tACGCTTATCAGGCCTACAATGTGTGTTGCAATTTACTGATTTCTTTGGATCTTGTAGGCCGGCTAAGGCGTTTACGCCGCATCCGGCCTGGCGCAACGTACTCGATAATAGAAATCTGGCGGCAACCCCACGCTGCCTATTAATGCtt                                  >  2:154436/1‑149 (MQ=37)
                                                                                                                  aCGCTTATCAGGCCTACAAGGTGTGTTGCAATTTACTGATTTCTTTGGATCTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAAGCAACGTACTCGATATTAGCAATTTGGCGGCAACCCAAAGTTGCCGATTAATGATTa                                 <  1:275557/149‑1 (MQ=255)
                                                                                                                    gCTTATCAGGCCTACAATGTGTGTTGCAATTTACTGATTTCTTTGGATCTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAAGCAACGTACTCGATATTAGCAATTTGGCGGCAACCCAAAGTTGCCGATTAATGATTACt                               >  2:21917/1‑149 (MQ=255)
                                                                                                                       tATCAGGCCTACAATGTGTGTTGCAATTTACTGATTTCTTTGGATCTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAAGCAACGTACTCGATATTAGCAATTTGGCGGCAACCCAAAGTTGCCGATTAATGATTACTTCg                            <  1:705634/149‑1 (MQ=255)
                                                                                                                        aTCAGGCCTACAATGTGTGTTGCAATTTACTGATTTCTTTGGATCTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAAGCAACGTACTCGATATTAGCAATTTGGCGGCAACCCAAAGTTGCCGATCAATGATATCTTCGc                           >  2:423515/1‑149 (MQ=255)
                                                                                                                         tCAGGCCTACAATGTGTGTTGCAATTTACTGATTTCTTTGGATCTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAAGCAACGTACTCGATATTAGCAATTTGGCGGCAACCCAAAGTTGCCGATTAATGATTACTTCGCa                          <  1:194232/149‑1 (MQ=255)
                                                                                                                            ggCCTACAATGTGTGTTGCAATTTACTGATTTCTTTGGATCTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAAGCAACGTACTCGATATTAGCAATTTGGCGGCAACCCAAAGTTGCCGATTAATGATTACTTCGCATCg                       >  1:395584/1‑149 (MQ=255)
                                                                                                                             gCCTACAATGTGTGTTGCAATTTACTGATTTCTTTGGATCTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAAGCAACGTACTCGATATTAGCAATTTGGCGGCAACCCAAAGTTGCCGATTAATGATTACTTCGCATCgg                      <  1:301228/149‑1 (MQ=255)
                                                                                                                                tACAATGTGTGTTGCAATTTACTGATTTCTTTGGATCTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAAGCAACGTACTCGATATTAGCAATTTGGCGGCAACCCAAAGTTGCCGGTTAATGATATCTTCGCATCGGTCg                   >  1:596748/1‑149 (MQ=255)
                                                                                                                                tACAATGTGTGTTGCAATTTACTGATTTCTTTGGATCTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAAGCAACGTACTCGATATTAGCAATTTGGCGGCAACCCAAAGTTGCCGATTAATGATTACTTCGCATCGGTCg                   >  2:597912/1‑149 (MQ=255)
                                                                                                                                tACAATGTGTGTTGCAATTTACTGATTTCTTTGGATCTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAAGCAACGTACTCGATATTAGCAATTTGGCGGCAACCCAAAATTGCCGATTAATAATTACTTCCAATCGGTcc                   >  1:668082/1‑148 (MQ=255)
                                                                                                                                   aaTGTGTGTTGCATTTTCCTGCTTTCTTTGGACCTTGTAGGCCGGAAAAAGGGTTTGCGCCGCATCCGGCATGAAGCAACGTACTCGATATTAGCAATTTGGCGGCAACCCAAAGTTGCCGATTAATGATTACTTCGCATCGGTCgccg                <  1:219027/149‑1 (MQ=255)
                                                                                                                                   aaTGTGTGTTGCAATTTGCTGATTTCTTTGGATCGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAAGCAACGTACTCGATATTAGCAATTTGGCGGCAACCCAAAGTTGCCGATTAATGATTACTTCGCATCGGTCgccg                <  1:486288/149‑1 (MQ=255)
                                                                                                                                         tgtTGCAATTTACTGATTTCTTTGGATCTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAAGCAACGTACTCGATATTAGCAATTTGGCGGCAACCCAAAGTTGCCGATTAATGATTACTTCGCATCGTTCgccgtgccgt          >  2:316374/1‑149 (MQ=255)
                                                                                                                                         tgtTGCAATTTACTGATTTCTTTGGATCTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAAGCAACGTACTCGATATTAGCAATTTGGCGGCAACCCAAAGTTGCCGATTAATGATTACTTCGCATCGGTCgccgtgccgt          >  2:780192/1‑149 (MQ=255)
                                                                                                                                          gtTGCAACTTACTGATTTCTTTGGATCTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAAGCAACGTACTCGATATTAGCAATTTGGCGGCAACCCAAAGTTGCCGATTAATGATTACTTCGCATCGGTCGCCGTGCCGtt         <  1:584380/149‑1 (MQ=255)
                                                                                                                                          gtTGCAACTTACTGATTTCTTTGGATCTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAAGCAACGTACTCGATATTAGCAATTTGGCGGCAACCCAAAGTTGCCGATTAATGATTACTTCGCATCGGTCGCCGTGCCGtt         <  1:124822/149‑1 (MQ=255)
                                                                                                                                          gtTGCAACTTACTGATTTCTTTGGATCTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAAGCAACGTACTCGATATTAGCAATTTGGCGGCAACCCAAAGTTGCCGATTAATGATTACTTCGCATCGGTCGCCGTGCCGtt         <  1:701151/149‑1 (MQ=255)
                                                                                                                                             gCAATTTACTGATTTCTTTGGGTCTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAAGCAACGTACTCGATATTAGCAATTTGGCGGCAACCCAAAGTTGCCGATTAATGATTACTTCGCATCGGTCGCCGTGCCGTTGGc      >  2:255984/1‑149 (MQ=255)
                                                                                                                                              cAATTTACTGATTTCTTTGGCTCTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAAGCAACGTACTCGATATTAGCAATTTGGCGGCAACCCAAAGTTGCCGATTAATGATTACTTCGCATCGGTCCCCGTGCCGTTGGCg     >  1:384373/1‑149 (MQ=255)
                                                                                                                                              cAATTTACTGATTTCTTTGGATCTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAAGCAACGTACTCGATATTAGCAATTTGGCGGCAACCCAACGTTGCCGCTTAATGGTTGCTTCGCATGGGTGGCCGTGCCGTTGGCg     >  1:719140/1‑149 (MQ=38)
                                                                                                                                                 tttACTGATTTCTTTGGATCTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAAGCAACGTACTCGATATTAGCAATTTGGCGGCAACCCAAAGTTGCCGATTAATGATTACTTCGCATCGGTCGCCGTGCCGTTGGCGTGc  <  1:634888/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
TCGCGGTGTGATGACGGCGTTTATGCAGGCGCTGGGGTTTACGGCACAACAGGGCGGCTGGGAGAAGTGTTAATCGTCAAGTTTGATTTCAAAAGTGATATTGCCTGATGCGCTACGCTTATCAGGCCTACAATGTGTGTTGCAATTTACTGATTTCTTTGGATCTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGAAGCAACGTACTCGATATTAGCAATTTGGCGGCAACCCAAAGTTGCCGATTAATGATTACTTCGCATCGGTCGCCGTGCCGTTGGCGTGC  >  NC_000913/3598294‑3598587

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: