Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 417,391 G→C A84P (GCC→CCC)  phoB → response regulator in two‑component regulatory system with PhoR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913417,3910GC98.3% 164.9 / ‑6.1 58A84P (GCC→CCC) phoBresponse regulator in two‑component regulatory system with PhoR
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base C (31/26);  total (32/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

CTGTGAATCAACTGAATGAACCCTGGCCGGATTTAATTCTCCTCGACTGGATGTTACCTGGCGGCTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGA  >  NC_000913/417254‑417527
                                                                                                                                         |                                                                                                                                        
ttgtgAATCAACAGAAGGAACCCTGGCCGGATTTAATGCACCACGACTGGATGTGACCTGGCGGCTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCCCCaga                                                                                                                                     <  1:378039/142‑1 (MQ=255)
    gAATCGACTGCTTGGACCCTGGCCCGTTTTGATTCTCCTCGTCTGTATGTTACCTGGCGGCGCCGGTATCCAGTGCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCCCCAGAgggg                                                                                                                                 <  2:472558/143‑1 (MQ=255)
     aaTCAACTGAATGAACCCTGGCCGGATTTAATTCTCCTCGACTGGATGTTACCTGGCGGCTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCCCCAGAGGGga                                                                                                                                >  1:862829/1‑143 (MQ=255)
       tctgcTGAATGCACCCTGGTCGGGTTCAATTCCGCTCTGCTGTGTCTGCTCCGCCGGCTCCGGTATCCAGTTCATCAAACACCTCAAGGGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCCCCAGAGGGgaag                                                                                                                              <  2:640711/139‑1 (MQ=255)
       tCAACTGAATGAACCCTGGCCGGATTTAATTCTCCTCGACTGGATGTTACCTGGCGGCTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCCCCAGCGGGgaag                                                                                                                              >  2:589275/1‑143 (MQ=255)
         aaCTGAGGGAACCCTGGCCGGCTTTAATTCTCCTCGACTGGATGTTACCTGGCGGCTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCCCCAGAGGGgaagaa                                                                                                                            <  2:95770/143‑1 (MQ=255)
         aaCTGAATGAACCCTGGCCGGATTTAATTCTCCTCGACTGGATGTTACCTGGCGGCTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCCCCAGAGGGgaagaa                                                                                                                            <  1:448112/143‑1 (MQ=255)
         aaCTGAATGAACCCTGGCCGGATTTAATTCTCCTCGACTGGATGTTACCTGGCGGCTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGAGGACCCGGGATATTCCAGTGGTGATGTTGACCCCCAGAGGGgaagaa                                                                                                                            >  1:76715/1‑143 (MQ=255)
          aCTGAATGAACCCTGGCCGGATTTAATTCTCCTCGACTGGATGTTACTCGGCGGCTCCGGTATCCAGTTCATCAAACACCTCACGCGCGAGACGATGACGCGGGATATTCCAGTGGTGATGTTGACGCCCGGAGGGcaaggag                                                                                                                           >  2:348226/1‑143 (MQ=255)
          aCTGAATGAACCCTGGCCGGATTTAATTCTCCTCGACTGGATGTTACCTGGCGGCTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCCCCAGAGGGgaagaag                                                                                                                           <  2:457309/143‑1 (MQ=255)
           cTGAATGAACCCTGGCCGGATTTAATTCTCCTCGACTGGATGTTACCTGGCGGCTCCGGTATGCAGTTCATCAAAAACCTCAAGCGCGAGTCGATGACCCGGGATATTGCAGTGGTCGTGTTGACCGCCGCAGGGgaaggaga                                                                                                                          >  1:147845/1‑143 (MQ=255)
            tgaatgaaCCCTGGCCGGATTTAATTCTCCTCGACTGGATGTTACCTGGCGGCTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCCCCAGAGGGgaagaagaa                                                                                                                         <  2:697992/143‑1 (MQ=255)
            tgaatgaaCCCTGGCCGGATTTAATTCTCCTCGACTGGATGTTACCTGGCGGCTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCCCCAGAGGGgaagaagaa                                                                                                                         <  1:815520/143‑1 (MQ=255)
                 gaaCCCTGGCCGGATTTAATTCTCCTCGACTGGATGTTACCTGGCGGCTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCCCCAGAGGGGaagaagaagcacg                                                                                                                    >  2:765017/1‑139 (MQ=255)
                      cTGGCCGGATTTAATTCTCCTCGACTGGATGTTACCTGGCGGCTCCGGTATCCAGTTGATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCCCCAGAGGGGAAGAAGAAGATCGCGTgc                                                                                                               >  1:668900/1‑143 (MQ=255)
                      cTGGCCGGATTTAATTCTCCTCGACTGGATGTTACCTGGCGGCTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCCCCAGAGGGGAAGAAGAAGATCGCGTgc                                                                                                               >  1:340709/1‑143 (MQ=255)
                                  aaTTCTCCTCGACTGGATGTTACCTGGCGGCTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGCCGATGACCCGGGATATTCCAGTGGTGATGTTGACCCCCAGAGGGCCAGAAGAAGTGCGCGTGCGCGGCCTTGaaa                                                                                                   >  1:325079/1‑143 (MQ=255)
                                    gtatcttcGACTGGATGTTACCTGGCTGCTCCGGTCTACTAGCCATCAAACACCTCAAGCTCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCCCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAAcc                                                                                                 <  2:181897/137‑1 (MQ=255)
                                          tcGACTGGATGTTACCTGGCGGCTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTTCAGTGGTGATGTTGACCCCCAGACGGGAGGCAGAAGATCGCGAGCGCGTCCTTAAACGCGGagcg                                                                                           >  2:30637/1‑143 (MQ=255)
                                          tcGACTGGATGTTACCTGGCGGCTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCCCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGcgcg                                                                                           >  1:34770/1‑143 (MQ=255)
                                           cGACTGAATGTTACCTGGCGGCTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATCCCAGTGGTGAGGTTGACCCCAAGAGGGGAAGAAGAAGATCGCGTGCGGGGCCTTGAACCAGGCGAgg                                                                                          >  2:819286/1‑143 (MQ=255)
                                             aCTCGAGGTTACCTGGCGGCTCCGGTATCCAGTTCATCAAACACCTCAGGCGCGAGTCGATGGCCCGGGATATTCCAGTGGTGATGTTGACCCCCAGAGGGGAAGAAGAAGAGCGCGTGCGCGGCCTTGAAACCGGCGCGGAt                                                                                        <  2:354796/143‑1 (MQ=255)
                                               gggAGGTTCCCGGGCGGCTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCCCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGa                                                                                      <  2:610259/142‑1 (MQ=255)
                                                              ggcTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGAACCCCAGAGGGGAAGCAGAAGATCGCGTGAGCGGCCTTGAGACCGGCGAGGCTGACTGTATccccgaccc                                                                       >  1:173411/1‑138 (MQ=255)
                                                                 tCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCCCCAGAGGGGCAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCACGGATGACTATATCACCAAGCCGtt                                                                    >  2:367166/1‑143 (MQ=255)
                                                                 tCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCCCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATAACCACGctgcg                                                                    >  1:590711/1‑139 (MQ=255)
                                                                   cGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCCCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGtttt                                                                  >  2:647095/1‑143 (MQ=255)
                                                                         ccAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCCCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGa                                                            >  1:187776/1‑143 (MQ=255)
                                                                         ccAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCCCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGa                                                            <  2:187776/143‑1 (MQ=255)
                                                                           aGTTCAGCAGCCGCCTCAAGCGCGAGTCGATGACCCGGGATATCCCAGTGGTCAGGTTGACCCCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAg                                                          <  2:688865/143‑1 (MQ=255)
                                                                                  cacacacCTCCCGCGCGAGTCGCTGACCCGGGATATTCCAGTGGTGATGTTGACCCCCAGCGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCCAGCCGTTTTCGCCGAAGGAGCtgg                                                   <  1:127282/140‑1 (MQ=255)
                                                                                   aaaCACCTCAAGCGCGAGTCGATGACACGGGATATTCCAGTGGTGATGTTGACCCCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATGTCACCCAGCCGTGTGCGCAGCAGGAGCcggc                                                  >  2:837788/1‑142 (MQ=255)
                                                                                        ccgcAAGCGCGAGTCGATGCCCCGGGATAGTCCAGTGGCGATGTGGACCCCCAGAGGGGAAGAAGAAGACCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGgcgc                                             <  1:307794/140‑1 (MQ=255)
                                                                                            aaGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCCCCAGAGCGGAAGAAAAAGATAGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCTTTTTAGCCGAAAGAACTGGTGGCGCGAAt                                         >  2:278050/1‑143 (MQ=255)
                                                                                              gcgcGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCCCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGGATCa                                       >  1:836486/1‑143 (MQ=255)
                                                                                                gcgAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCCCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCaaa                                     <  1:690888/143‑1 (MQ=255)
                                                                                                gcgAGTCGATGACCCGGGATATCCCCGTGGTGATGTTGACCCCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCaaa                                     <  1:363933/143‑1 (MQ=255)
                                                                                                   aGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCCCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGACGGAGCTGGTGGCGCGAATCAAAGCg                                  >  2:65340/1‑143 (MQ=255)
                                                                                                    gTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCCCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCgg                                 >  1:752444/1‑143 (MQ=255)
                                                                                                    gTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCCCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCgg                                 <  2:475920/143‑1 (MQ=255)
                                                                                                        aTGACCCGGGATATTCCAGTGGTGATGTTGACCCCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAAt                             <  1:716804/143‑1 (MQ=255)
                                                                                                        aTGACCCGGGAGATTCCAGTGGTGATGTGGACCCCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCGGGTGGCGGGAGTGAAAGCGGGaaa                             >  2:138226/1‑142 (MQ=255)
                                                                                                            cgcGGGATGTGCCGGTGGTGAGGTTGACCCCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGCAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATgcgc                         <  2:275079/141‑1 (MQ=255)
                                                                                                              ggggATATGCCGGTGGTGATGTGGACCCCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTGTGCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCg                       <  2:353165/142‑1 (MQ=255)
                                                                                                              cGGGATATTCCAGTGGTGATGTTGACCCCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGACACCGGCGCGGATGACTATATCACCAAGCCGTTTTCTCCGCGGGAGCGGGTGGCGCGAATCAAAGCGCTAATGCCCCg                       >  2:97814/1‑143 (MQ=255)
                                                                                                              cGGGATATTCCAGTGGTGATGTTGACCCCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCg                       <  1:107412/143‑1 (MQ=255)
                                                                                                              cGGGATATTCCAGTGGTGATGTTGACCCCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGAGGACTATATCACCAAGCCGTTTTCGCCGAGGGAGCTGGTGGCGCGACCCAAAGCGGTAATGCCCCg                       >  2:599560/1‑143 (MQ=255)
                                                                                                                    attcgggtggtgATGTTGACCCCCAGGTGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTc                 <  2:256818/137‑1 (MQ=255)
                                                                                                                     ttCCAGTGGTGATGTTGACCCCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCg                <  1:553227/143‑1 (MQ=255)
                                                                                                                     ttCCAGTGGTGATGTTGACCCCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCg                <  1:781960/143‑1 (MQ=255)
                                                                                                                     ttCCAGTGGTGATGTTGACCCCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCGGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTca                >  1:647384/1‑142 (MQ=255)
                                                                                                                      tCCAGTGGTGATGTTGACCCCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATGTCGc               >  1:136869/1‑143 (MQ=255)
                                                                                                                        cAGTGGTGATGTTGACCCCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCa             >  1:414951/1‑143 (MQ=255)
                                                                                                                             gtgATGTTGACCCCCAGAGGGGAAGAAGAAGATCGCGGGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGCCTCACAGCGGTAATGCGCCGAATTTCGCCAACGGc        >  2:572546/1‑143 (MQ=255)
                                                                                                                              tgATGTTGACCCCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTACGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCg       >  2:433688/1‑143 (MQ=255)
                                                                                                                                aTGTTGACCCCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGCAACCGGCGCGGATGACTATATCCCCAAGCCGGTCTCGCCGAAGGCGCGGGTGGCGCGAACAAAAGCGGCAACGCGCCGGAATTCGCCAACGGAggg     >  2:659971/1‑142 (MQ=255)
                                                                                                                                  gTTGACCCCCAGAGGCGAGCAAGAAGATCGCGCGTGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTgg   <  1:576256/143‑1 (MQ=255)
                                                                                                                                   ttGACCCCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATAACCAAGCCAAAATCGCAGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGGGGa  >  1:331580/1‑143 (MQ=255)
                                                                                                                                         |                                                                                                                                        
CTGTGAATCAACTGAATGAACCCTGGCCGGATTTAATTCTCCTCGACTGGATGTTACCTGGCGGCTCCGGTATCCAGTTCATCAAACACCTCAAGCGCGAGTCGATGACCCGGGATATTCCAGTGGTGATGTTGACCGCCAGAGGGGAAGAAGAAGATCGCGTGCGCGGCCTTGAAACCGGCGCGGATGACTATATCACCAAGCCGTTTTCGCCGAAGGAGCTGGTGGCGCGAATCAAAGCGGTAATGCGCCGTATTTCGCCAATGGCGGTGGA  >  NC_000913/417254‑417527

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