Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,778,470 A→G Y278H (TAC→CAC)  ydiP ← putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,778,4700AG98.4% 227.8 / NA 64Y278H (TAC→CAC) ydiPputative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (34/29);  total (34/30)
Fisher's exact test for biased strand distribution p-value = 4.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTAAATTGTCTACGGTAATCGCTGGGCGAACAACCGACATGGCGCAAAAACAGTTTGGCAAAGTGATCGACATTTTCATAACCCACACGCCAGGAGATCTCTGCCTGTGATAA  >  NC_000913/1778371‑1778566
                                                                                                   |                                                                                                
aCAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTAAATTGTCTACGGTa                                                                                                  >  2:278433/1‑100 (MQ=255)
   aaCTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTAAATTGTCTACGGTGATc                                                                                               <  1:124892/100‑1 (MQ=255)
     cTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTAAATTGTCTACGGTGATCGc                                                                                             <  1:1207887/100‑1 (MQ=255)
     cTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTAAATTGTCTACGGTGATCGc                                                                                             >  1:853171/1‑100 (MQ=255)
        cgccgcCAGATTTATCCAACAAGACTTTCCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTAAATTGTCTACGGTGATCGCTgg                                                                                          <  1:450606/100‑1 (MQ=255)
         gccgccAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTAAATTGTCTACGGTGATCGCTggg                                                                                         <  2:1216562/100‑1 (MQ=255)
          ccgccAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTAAATTGTCTACGGTGATCGCTGGGc                                                                                        >  1:845828/1‑100 (MQ=255)
           cgccAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTAAATTGTCTACGGTGATCGCTGGGCg                                                                                       >  2:823876/1‑100 (MQ=255)
           cgccAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTAAATTGTCTACGGTGATCGCTGGGCg                                                                                       >  2:355510/1‑100 (MQ=255)
           cgccAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTAAATTGTCTACGGTGATCGCTGGGCg                                                                                       >  1:13762/1‑100 (MQ=255)
            gccAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTAAATTGTCTACGGTGATCGCTGGGCg                                                                                       <  2:999179/99‑1 (MQ=255)
            gccAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTAAATTGTCTACGGTGATCGCTGGGCGa                                                                                      >  1:221882/1‑100 (MQ=255)
             gcAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTAAATTGTCTACGGTGATCGCTGGGCGaa                                                                                     >  2:237102/2‑100 (MQ=255)
              cAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTAAATTGTCTACGGTGATCGCTGGGCGaa                                                                                     >  2:881845/1‑99 (MQ=255)
                  tttATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTAAATTGTCTACGGTGATCGCTGGGCGaacaacn                                                                                <  1:278644/100‑2 (MQ=255)
                  tttATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTAAATTGTCTACGGTGATCGCTGGGCGAACAAcc                                                                                <  1:321740/100‑1 (MQ=255)
                  tttATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTAAATTGTCTACGGTGATCGCTGGGCGAACAAcc                                                                                >  1:423942/1‑100 (MQ=255)
                     aTCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTAAATTGTCTACGGTAATCGCTGGGCGAACAAcc                                                                                <  1:278436/97‑1 (MQ=255)
                               aCTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTAAATTGTCTACGGTGATCGCTGGGCGAACAACCGACATGGCGCaaa                                                                   >  1:615074/1‑100 (MQ=255)
                                  tACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTAAATTGTCTACGGTGATCGCTGGGCGAACAACCGACATGGCGCAAAAAc                                                                >  2:741568/1‑100 (MQ=255)
                                  tACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTAAATTGTCTACGGTGATCGCTGGGCGAACAACCGACATGGCGCAAAAAc                                                                >  2:159875/1‑100 (MQ=255)
                                      ggTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTAAATTGTCTACGGTGATCGCTGGGCGAACAACCGACATGGCGCAAAAACAGtt                                                            <  2:1009671/100‑1 (MQ=255)
                                        ttGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTAAATTGTCTACGGTGATCGCTGGGCGAACAACCGACATGGCGCAAAAACAGTTTg                                                          <  2:1302781/100‑1 (MQ=255)
                                         tGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTAAATTGTCTACGGTGATCGCTGGGCGAACAACCGACATGGCGCAAAAACAGTTTgg                                                         <  2:1035891/100‑1 (MQ=255)
                                         tGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTAAATTGTCTACGGTGATCGCTGGGCGAACAACCGACATGGCGCAAAAACAGTTTgg                                                         >  1:280818/1‑100 (MQ=255)
                                          gAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTAAATTGTCTACGGTGATCGCTGGGCGAACAACCGACATGGCGCAAAAACAGTTTgg                                                         <  1:221213/99‑1 (MQ=255)
                                            ggAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTAAATTGTCTACGGTGATCGCTGGGCGAACAACCGACATGGCGCAAAAACAGTTTGGCaa                                                      <  1:1105785/100‑1 (MQ=255)
                                              aaaTTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTAAATTGTCTACGGTGATCGCTGGGCGAACAACCGACATGGCGCAAAAACAGTTTGGCAAAg                                                    <  1:182751/100‑1 (MQ=255)
                                              aaaTTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTAAATTGTCTACGCTGATCGCTGGGCGAACAACCGACATGGCGCAAAAACAGTTTGGCAAAg                                                    >  2:1295348/1‑100 (MQ=255)
                                                     gATAGGATTTCTTGTTCCGCAAAACAGTTTTTAAATTGTCTACGGTGATCGCTGGGCGAACAACCGACATGGCGCAAAACCAGTTTGACAAAGTGATCGa                                             >  1:1253425/1‑100 (MQ=255)
                                                      aTAGGATTTCTTGTTCCGCAAAACAGTTTTTAAATTGTCTACGGTGATCGCTGGGCGAACAACCGACATGGCGCAAAAACAGTTTGGCAAAGTGATCGAc                                            >  1:230980/1‑100 (MQ=255)
                                                         ggATTTCTTGTTCCGCAAAACAGTTTTTAAATTGTCTACGGTGATCGCTGGGCGAACAACCGACATGGCGCAAAAACAGTTTGGCAAAGTGATCGACAt                                          >  1:1080917/1‑99 (MQ=255)
                                                             ttCTTGTTCCGCAAAACAGTTTTTAAATTGTCTACGGTGATCGCTGGGCGAACAACCGACATGGCGCAAAAACAGTTTGGCAAAGTGATCGACATTTTCa                                     <  2:515275/100‑1 (MQ=255)
                                                               cTTGTTCCGCAAAACAGTTTTTAAATTGTCTACGGTGATCGCTGGGCGAACAACCGACATGGCGCAAAAACAGTTTGGCAAAGTGATCGACATTTTCATa                                   >  2:504092/1‑100 (MQ=255)
                                                               cTTGTTCCGCAAAACAGTTTTTAAATTGTCTACGGTGATCGCTGGGCGAACAACCGACATGGCGCAAAAACAGTTTGGCAAAGTGATCGACATTTTCATa                                   >  1:861682/1‑100 (MQ=255)
                                                                ttGTTCCGCAAAACAGTTTTTAAATTGTCTACGGTGATCGCTGGGCGAACAACCGACATGGCGCAAAAACAGTTTGGCAAAGTGATCGACATTTTCATaa                                  >  2:1179041/1‑100 (MQ=255)
                                                                  gTTCCGCAAAACAGTTTTTAAATTGTCTACGGTGATCGCTGGGCGAACAACCGACATGGCGCAAAAACAGTTTGGCAAAGTGATCGACATTTTCATAAcc                                >  1:413777/1‑100 (MQ=255)
                                                                    tCCGCAAAACAGTTTTTAAATTGTCTACGGTGATCGCTGGGCGAACAACCGACATGGCGCAAAAACAGTTTGGCAAAGTGATCGACATTTTCATAACCca                              >  2:236274/1‑100 (MQ=255)
                                                                     ccGCAAAACAGTTTTTAAATTGTCTACGGTGATCGCTGGGCGAACAACCGACATGGCGCAAAAACAGTTTGGCAAAGTGATCGACATTTTCATAACCcac                             >  2:587935/1‑100 (MQ=255)
                                                                     ccGCAAAACAGTTTTTAAATTGTCTACGGTGATCGCTGGGCGAACAACCGACATGGCGCAAAAACAGTTTGGCAAAGTGATCGACATTTTCATAACCcac                             <  2:413773/100‑1 (MQ=255)
                                                                      cGCAAAACAGTTTTTAAATTGTCTACGGTGATCGCTGGGCGAACAACCGACATGGCGCAAAAACAGTTTGGCAAAGTGATCGACATTTTCATAACCcaca                            <  2:221880/100‑1 (MQ=255)
                                                                          aaaCAGTTTTTAAATTGTCTACGGTGATCGCTGGGCGAACAACCGACATGGCGCAAAAACAGTTTGGCAAAGTGATCGACATTTTCATAACCCACACGcc                        >  2:1278297/1‑100 (MQ=255)
                                                                          aaaCAGTTTTTAAATTGTCTACGGTGATCGCTGGGCGAACAACCGACATGGCGCAAAAACAGTTTGGCAAAGTGATCGACATTTTCATAACCCACACGcc                        <  1:618400/100‑1 (MQ=255)
                                                                            aCAGTTTTTAAATTGTCTACGGTGATCGCTGGGCGAACAACCGACATGGCGCAAAAACAGTTTGGCAAAGTGATCGACATTTTCATAACCCACACGCCAg                      <  1:1054603/100‑1 (MQ=255)
                                                                             cAGTTTTTAAATTGTCTACGGTGATCGCTGGGCGAACAACCGACATGGCGCAAAAACAGTTTGGCAAAGTGATCGACATTTTCATAACCCACACGCCAgg                     >  2:66905/1‑100 (MQ=255)
                                                                               gTTTTTAAATTGTCTACGGTGATCGCTGGGCGAACAACCGACATGGCGCAAAAACAGTTTGGCAAAGTGATCGACATTTTCATAACCCACACGCTAGgag                   <  1:454172/100‑1 (MQ=255)
                                                                               gTTTTTAAATTGTCTACGGTGATCGCTGGGCGAACAACCGACATGGCGCAAAAACAGTTTGGCAAAGTGATCGACATTTTCATAACCCACACGCCAGgag                   <  2:230978/100‑1 (MQ=255)
                                                                               gTTTTTAAATTGTCTACGGTGATCGCTGGGCGAACAACCGACATGGCGCAAAAACAGTTTGGCAAAGTGATCGACATTTTCATAACCCACACGCCAGgag                   >  2:888828/1‑100 (MQ=255)
                                                                               gTTTTTAAATTGTCTACGGTGATCGCTGGGCGAACAACCGACATGGCGCAAAAACAGTTTGGCAAAGTGATCGACATTTTCATAACCCACACGCCAGgag                   <  1:261455/100‑1 (MQ=255)
                                                                               gTTTTTAAATTGTCTACGGTGATCGCTGGGCGAACAACCGACATGGCGCAAAAACAGTTTGGCAAAGTGATCGACATTTTCATAACCCACACGCCAGgag                   <  2:571221/100‑1 (MQ=255)
                                                                               gTTTTTAAATTGTCTACGGTGATCGCTGGGCGAACAACCGACATGGCGCAAAAACAGTTTGGCAAAGTGATCGACATTTTCATAACCCACACGCCAGgag                   >  2:601823/1‑100 (MQ=255)
                                                                                 ttttAAATTGTCTACGGTGATCGCTGGGCGAACAACCGACATGGCGCAAAAACAGTTTGGCAAAGTGATCGACATTTTCATAACCCACACGCCAGGAGAt                 <  1:428003/100‑1 (MQ=255)
                                                                                  tttAAATTGTCTACGGTGATCGCTGGGCGAACAACCGACATGGCGCAAAAACAGTTTGGCAAAGTGATCGACATTTTCATAACCCACACGCCAGGAGAtc                >  1:803246/1‑100 (MQ=255)
                                                                                     aaaTTGTCTACGGTGATCGCTGGGCGAACAACCGACATGGCGCAAAAACAGTTTGGCAAAGTGATCGACATTTTCATAACCCACACGCCAGGAGATctct             >  1:484524/1‑100 (MQ=255)
                                                                                     aaaTTGTCTACGGTGATCGCTGGGCGAACAACCGACATGGCGCAAAAACAGTTTGGCAAAGTGATCGACATTTTCATAACCCACACGCCAGGAGATctct             >  1:944860/1‑100 (MQ=255)
                                                                                        ttGTCTACGGTGATCGCTGGGCGAACAACCGACATGGCGCAAAAACAGTTTGGCAAAGTGATCGACATTTTCATAACCCACACGCCAGGAGATCTCTGcc          <  1:2648/100‑1 (MQ=255)
                                                                                        ttGTCTACGGTGATCGCTGGGCGAACAACCGACATGGCGCAAAAACAGTTTGGCAAAGTGATCGACATTTTCATAACCCACACGCCAGGAGATCTCTGcc          <  2:405991/100‑1 (MQ=255)
                                                                                        ttGTCTACGGTGATCGCTGGGCGAACAACCGACATGGCGCAAAAACAGTTTGGCAAAGTGATCGACATTTTCATAACCCACACGCCAGGAGATCTCTGcc          <  2:13762/100‑1 (MQ=255)
                                                                                        ttGTCTACGGTGATCGCTGGGCGAACAACCGACATGGCGCAAAAACAGTTTGGCAAAGTGATCGACATTTTCATAACCCACACGCCAGGAGATCTCTGcc          <  1:741573/100‑1 (MQ=255)
                                                                                            cTACGGTGATCGCTGGGCGAACAACCGACATGGCGCAAAAACAGTTTGGCAAAGTGATCGACATTTTCATAACCCACACGCCAGGAGATCTCTGCCtgtg      <  1:1278308/100‑1 (MQ=255)
                                                                                              aCGGTGATCGCTGGGCGAACAACCGACATGGCGCAAAAACAGTTTGGCAAAGTGATCGACATTTTCATAACCCACACGCCAGGAGATCTCTGCCTGTGAt    >  1:419477/1‑100 (MQ=255)
                                                                                              aCGGTGATCGCTGGGCGAACAACCGACATGGCGCAAAAACAGTTTGGCAAAGTGATCGACATTTTCATAACCCACACGCCAGGAGATCTCTGCCTGTGAt    >  2:338169/1‑100 (MQ=255)
                                                                                              aCGGTGATCGCTGGGCGAACAACCGACATGGCGCAAAAACAGTTTGGCAAAGTGATCGACATTTTCATAACCCACACGCCAGGAGATCTCTGCCTGTGAt    >  1:1205680/1‑100 (MQ=255)
                                                                                               cGGTGATCGCTGGGCGAACAACCGACATGGCGCAAAAACAGTTTGGCAAAGTGATCGACATTTTCATAACCCACACGCCAGGAGATCTCTGCCTGTGATa   <  1:823882/100‑1 (MQ=255)
                                                                                                ggTGATCGCTGGGCGAACAACCGACATGGCGCAAAAACAGTTTGGCAAAGTGATCGACATTTTCATAACCCACACGCCAGGAGATCTCTGCCTGTGATaa  >  1:447144/1‑100 (MQ=255)
                                                                                                   |                                                                                                
ACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTAAATTGTCTACGGTAATCGCTGGGCGAACAACCGACATGGCGCAAAAACAGTTTGGCAAAGTGATCGACATTTTCATAACCCACACGCCAGGAGATCTCTGCCTGTGATAA  >  NC_000913/1778371‑1778566

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