Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,290,448 A→G 9.2% intergenic (‑268/‑52) yejO ← / → narP pseudogene, autotransporter outer membrane homology;putative transport; Not classified; putative ATP‑binding component of a transport system/response regulator in two‑component regulatory system with NarQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,290,4480AG 9.2% 161.3 / 6.0 54intergenic (‑268/‑52)yejO/narPpseudogene, autotransporter outer membrane homology;putative transport; Not classified; putative ATP‑binding component of a transport system/response regulator in two‑component regulatory system with NarQ
Reads supporting (aligned to +/- strand):  ref base A (23/26);  new base G (3/2);  total (26/28)
Fisher's exact test for biased strand distribution p-value = 6.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

TTGATTTTTTAGTTGTTTTTCTTGATGAGAAGCTGATGCAAAATTCCGTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCAACGCTACGCTCAAACATAATGATTCTAATAAAACCTCAGGAGACTACTATGCCTGAAGCAACACCTTTTCAGGTGATGATTGTGGATGATCATCCACTTATGCGACGCGGTGTTCGTCAGTTACTGGAGCTTGATCCTGGCT  >  NC_000913/2290300‑2290593
                                                                                                                                                    |                                                                                                                                                 
ttgattTTTTAGTTGTTTTTCTTGATGAGAAGCTGATGCAAAATTCCGTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTa                                                                                                                                                   >  4:118120/1‑149 (MQ=255)
     tttttAGTTGTTTTTCTTGATGAGAAGCTGATGCAAAATTCCGTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCAac                                                                                                                                              <  2:641321/149‑1 (MQ=255)
            ttgttTTTCTTGATGAGAAGCTGATGCAAAATTCCGTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTGGTAacgctacca                                                                                                                                       >  4:99710/1‑147 (MQ=255)
             tgttTTTCTTGATGAGAAGCTGATGCAAAATTCCGTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTGGTAacgctaccag                                                                                                                                      <  3:99710/149‑4 (MQ=255)
              gttTTTCTTGATGAGAAGCTGATGCAAAATTCCGTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCaa                                                                                                                                               >  1:188346/1‑139 (MQ=255)
              gttTTTCTTGATGAGAAGCTGATGCAAAATTCCGTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCaa                                                                                                                                               <  2:188346/139‑1 (MQ=255)
              gttTTTCTTGATGAGAAGCTGATGCAAAATTCCGTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCAACGCTACGCTc                                                                                                                                     <  3:338413/149‑1 (MQ=255)
                    cTTGATGAGAAGCTGATGCAAAATTCCGTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCAACGCTACGCTCAAACAt                                                                                                                               >  1:266402/1‑149 (MQ=255)
                      tgatgaGAAGCTGATGCAAAATTCCGTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCAACGCTACGCTCAAACATaa                                                                                                                             <  1:35332/149‑1 (MQ=255)
                      tgatgaGAAGCTGATGCAAAATTCCGTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCAACGCTACGCTCAAACATaa                                                                                                                             <  3:118120/149‑1 (MQ=255)
                               gCTGATGCAAAATTCCGTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTTCACTACATTTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCAACGCTACGCTCAAACATAATGATTCtaa                                                                                                                    >  1:182652/1‑149 (MQ=255)
                               gCTGATGCAAAATTCCGTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCAACGCTACGCTCAAACATAATGATTCtaa                                                                                                                    <  1:319398/149‑1 (MQ=255)
                                cTGATGCAAAATTCCGTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCAACGCTACGCTCAAACATAATGATTCtaat                                                                                                                   <  1:591642/149‑1 (MQ=255)
                                         aaTTCCGTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCAACGCTACGCTCAAACATAATGATTCTAATAAAACCTCa                                                                                                          <  1:322999/149‑1 (MQ=255)
                                               gTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCAACGCTACGCTCAAACATAATGATTCTAATAAAACCTCAGGAGac                                                                                                    >  3:624049/1‑149 (MQ=255)
                                                   ttATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTGGTAacgctacc                                                                                                                                        >  2:352129/1‑108 (MQ=255)
                                                   ttATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTGGTAacgctacc                                                                                                                                        <  1:352129/109‑2 (MQ=255)
                                                         tGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCAACGCTACGCTCAAACATAATGATTCTAATAAAACCTCAGGAGACTACTATGCCt                                                                                          >  1:109132/1‑149 (MQ=255)
                                                          gAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCAACGCTACGCTCAAACATAATGATTCTAATAAAACCTCAGGAGACTACTATGCCTg                                                                                         >  4:329253/1‑149 (MQ=255)
                                                            aaaTGATGCCAAAGCGAACGACAAGGTTGTAGTTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCAACGCTACGCTCAAACATAATGATTCTAATAAAACCTCAGGAGACTACTATGCCTGaa                                                                                       <  4:624049/149‑1 (MQ=255)
                                                            aaaTGATGCCAAAGCGAACGACAAGGTTGTAGTTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCAACGCTACGCTCAAACATAATGATTCTAATAAAACCTCAGGAGACTACTATGCCTGaa                                                                                       <  2:266402/149‑1 (MQ=255)
                                                                           gAACGACAAGGTTGTAGTTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCAACGCTACGCTCAAACATAATGATTCTAATAAAACCTCAGGAGACTACTATGCCTGAAGCAACACCTTTTCAg                                                                        <  4:171437/149‑1 (MQ=255)
                                                                            aaCGACAAGGTTGTAGTTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCAACGCTACGCTCAAACATAATGATTCTAATAAAACCTCAGGAGACTACTATGCCTGAAGCAACACCTTTTCAgg                                                                       >  4:24404/1‑149 (MQ=255)
                                                                               gACAAGGTTGTAGTTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCAACGCTACGCTCAAACATAATGATTCTAATAAAACCTCAGGAGACTACTATGCCTGAAGCAACACCTTTTCAGGtga                                                                    >  1:177436/1‑149 (MQ=255)
                                                                                      ttGTAGTTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCAACGCTACGCTCAAACATAATGATTCTAATAAAACCTCAGGAGACTACTATGCCTGAAGCAACACCTTTTCAGGTGATGATtgt                                                             <  3:24404/149‑1 (MQ=255)
                                                                                            ttttCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCAACGCTACGCTCAAACATAATGATTCTAATAAAACCTCAGGAGACTACTATGCCTGAAGCAACACCTTTTCAGGTGATGATTGTGgatga                                                       <  1:377710/149‑1 (MQ=255)
                                                                                                cACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCAACGCTACGCTCAAACATAATGATTCTAATAAAACCTCAGGAGACTACTATGCCTGAAGCAACACCTTTTCAGGTGATGATTGTGGATGatcat                                                   >  1:146216/1‑149 (MQ=255)
                                                                                                 aCTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCAACGCTACGCTCAAACATAATGATTCTAATAAAACCTCAGGAGACTACTATGCCTGAAGCAACACCTTTTCAGGTGATGATTGTGGATGatcatc                                                  >  3:313192/1‑149 (MQ=255)
                                                                                                    aCATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCAACGCTACGCTCAAACATAATGATTCTAATAAAACCTCAGGAGACTACTATGCCTGAAGCAACACCTTTTCAGGTGATGATTGTGGATGATCATCCAc                                               <  2:146216/149‑1 (MQ=255)
                                                                                                      aTGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCAACGCTACGCTCAAACATAATGATTCTAATAAAACCTCAGGAGACTACTATGCCTGAAGCAACACCTTTTCAGGTGATGATTGTGGATGATCATCCAcaa                                             <  3:607083/149‑3 (MQ=255)
                                                                                                      aTGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCAACGCTACGCTCAAACATAATGATTCTAATAAAACCTCAGGAGACTACTATGCCTGAAGCAACACCTTTTCAGGTGATGATTGTGGATGATCATCCACtt                                             >  2:671869/1‑149 (MQ=255)
                                                                                                         tCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCAACGCTACGCTCAAACATAATGATTCTAATAAAACCTCAGGAGACTACTATGCCTGAAGCAACACCTTTTCAgg                                                                       <  2:535681/120‑1 (MQ=255)
                                                                                                         tCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCAACGCTACGCTCAAACATAATGATTCTAATAAAACCTCAGGAGACTACTATGCCTGAAGCAACACCTTTTCAgg                                                                       >  1:535681/1‑120 (MQ=255)
                                                                                                           catacataAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCAACGCTACGCTCAAACATAATGATTCTAATAAAACCTCAGGAGACTACTATGCCTGAAGCAACACCTTTTCAGGTGATGATTGTGGATGATCATCCACTTATGCg                                        >  3:64194/1‑149 (MQ=255)
                                                                                                              acataAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCAACGCTACGCGCAAACATAATGATTCTAATAAAACCTCAGGAGACTACTATGCCTGAAGCAACACCTTTTCAGGTGATGATTGTGGATGATCATCCACTTATGCGAcg                                     >  4:239285/1‑149 (MQ=255)
                                                                                                               cataAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCAACGCTACGCTCAAACATAATGCTTCTAATAAAACCTCAGGAGTCTACTATGCCTGAAGCAACACCTTTTCAGGTGATGATTGTGGATGATCATCCACTTATGCGAcgc                                    >  1:615386/1‑149 (MQ=255)
                                                                                                                ataAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCAACGCTACGCTCAAACATAATGATTCTAATAAAACCTCAGGAGACTACTATGCCTGAAGCAACACCTTTTCAGGTGATGATTGTGGATGATCATCCACTTATGCGAcgcg                                   >  1:544574/1‑149 (MQ=255)
                                                                                                                    aaTGGGGTAACATTCACGCGCCTGGTAGCGTTACCAACGCTACGCTCAAACATAATGATTCTAATAAAACCTCAGGAGACTACTATGCCTGAAGCAACACCTTTTCAGGTGATGATTGTGGATGATCATCCACTTATGCGACGCGgtgt                               >  2:490498/1‑149 (MQ=255)
                                                                                                                     aTGGGGTAACATTCACGCGCCTGGTAGCGTTACCAACGCTACGCTCAAACATAATGATTCTAATAAAACCTCAGGAGACTACTATGCCTGAAGCAACACCTTTTCAGGTGATGATTGTGGATGATCATCCACTTATGCGACGCGGTGtt                              <  2:177436/149‑1 (MQ=255)
                                                                                                                      tGGGGTAACATTCACGCGCCTGGTAGCGTTACCAACGCTACGCTCAAACATAATGATTCTAATAAAACCTCAGGAGACTACTATGCCTGAAGCAACACCTTTTCAGGTGATGATTGTGGATGATCATCCACTTATGCGACGCGGTGTTc                             <  4:430412/149‑1 (MQ=255)
                                                                                                                      tGGGGTAACATTCACGCGCCTGGTAGCGTTACCAACGCTACGCTCAAACATAATGATTCTAATAAAACCTCAGGAGACTACTATGCCTGAAGCAACACCTTTTCAGGTGATGATTGTGGATGATCATCCACTTATGCGACGCGGTGTTc                             <  3:49968/149‑1 (MQ=255)
                                                                                                                      tGGGGTAACATTCACGCGCCTGGTAGCGTTACCAACGCTACGCTCAAACATAACGATTCTAATAAAACCTCAGGAGACTACTATGCCTGAAGCAACACCTTTTCAGGTGATGATTGTGGATGATCATCCACTTATGCGACGCGGTGTTc                             <  1:490498/149‑1 (MQ=255)
                                                                                                                         ggTAACATTCACGCGCCTGGTAGCGTTACCAACGCTACGCTCAAACATAATGATTCTAATAAAACCTCAGGAGACTACTATGCCTGAAGCAACACCTTTTCAGGTGATGATTGTGGATGATCATCCACTTATGCGACGCGGTGTTCGTc                          <  2:544574/149‑1 (MQ=255)
                                                                                                                          gTAACATTCACGCGCCTGGTAGCGTTACCAACGCTACGCTCAAACATAATGATTCTAATAAAACCTCAGGAGACTACTATGCCTGAAGCAACACCTTTTCAGGTGATGATTGTGGATGATCATCCACTTATGCGACGCGGTGTTCGTc                          <  2:458925/148‑1 (MQ=255)
                                                                                                                          gTAACATTCACGCGCCTGGTAGCGTTACCAACGCTACGCTCAAACATAATGATTCTAATAAAACCTCAGGAGACTACTATGCCTGAAGCAACACCTTTTCAGGTGATGATTGTGGATGATCATCCACTTATGCGACGCGGTGTTCGTc                          >  1:458925/1‑148 (MQ=255)
                                                                                                                             aCATTCACGCGCCTGGTAGCGTTACCAACGCTACGCTCAAACATAATGATTCTAATAAAACCTCAGGAGACTACTATGCCTGAAGCAACACCTTTTCAGGTGATGATTGTGGATGATCATCCACTTATGCGACGCGGTGTTCGTCAGtt                      <  1:671869/149‑1 (MQ=255)
                                                                                                                               aTTCACGCGCCTGGTAGCGTTACCAACTCTACGCTCAAACATAATGCTTCTAATAAAACCTCAGTAGACTACTATGCCTGAAGCAACACCTTTTCAGGTGATGATTGTGgatga                                                       >  1:517415/1‑114 (MQ=255)
                                                                                                                               aTTCACGCGCCTGGTAGCGTTACCAACGCTACGCTCAAACATAATGATTCTAATAAAACCTCAGGAGACTACTATGCCTGAAGCAACACCTTTTCAGGTGATGATTGTGgatga                                                       <  2:517415/114‑1 (MQ=255)
                                                                                                                                     gcgcCTGGTAGGGTTACCAACGCTACGCTCAAACATAATGATTCTAATAAAACCTCAGGAGACTACTATGCCTGAAGCAACACCTTTTCAGGTGATGATTGTGGATGATCATCCACTTATGCGATGCGGTGTTCGTCAGTTACTGGAGc              >  3:456040/1‑149 (MQ=255)
                                                                                                                                         cTGGTAGCGTTACCAACGCTACGCTCAAACATAATGATTCTAATAAAACCTCAGGAGACTACTATGCCTGAAGCAACACCTTTTCAGGTGATGATTGTGGATGATCATCCACTTATGCGACGCGGTGTTCGTCAGTTACTGGAGCTTGa          >  1:626062/1‑149 (MQ=255)
                                                                                                                                           ggTAGCGTTACCAACGCTACGCTCAAACATAATGATTCTAATAAAACCTCAGGAGACTACTATGCCTGAAGCAACACCTTTTCAGGTGATGATTGTGGATGATCATCCACTTATGCGACGCGGTGTTCGTCAGTTACTGGAGCTTGATc        >  4:622092/1‑149 (MQ=255)
                                                                                                                                           cgTAGCGTTGGTAACGCTACGCTCAAACATAATGATTCTAATAAAACCTCAGGAGACTACTATGCCTGAAGCAACACCTTTTCAGGTGATGATTGTGGATGATCATCCACTTATGCGACGCGGTGTTCGTCAGTTACTGGAGCTTGATc        >  2:96810/2‑149 (MQ=255)
                                                                                                                                             tAGCGTTACCAACGCTACGCTCAAACATAATGATTCTAATAAAACCTCAGGAGACTACTATGCCTGAAGCAACACCTTTTCAGGTGATGATTGTGGATGATCATCCACTTATGCGACGCGGTGTTCGTCAGTTACTGGAGCTTGATCCt      <  4:191600/149‑1 (MQ=255)
                                                                                                                                              aGCGTTACCAACGCTACGCTCAAACATAATGATTCTAATAAAACCTCAGGAGACTACTATGCCTGAAGCAACACCTTTTCAGGTGATGATTGTGGATGATCATCCACTTATGCGACGCGGTGTTCGTCAGTTACTGGAGCTTGATCCTg     <  2:109132/149‑1 (MQ=255)
                                                                                                                                                 gTTACCAACGCTACGCTCAAACATAATGATTCTAATAAAACCTCAGGAGACTACTATGCCTGAAGCAACACCTTTTCAGGTGATGATTGTGGATGATCATCCACTTATGCGACGCGGTGTTCGTCAGTTACTGGAGCTTGATCCTGGct  >  4:201423/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
TTGATTTTTTAGTTGTTTTTCTTGATGAGAAGCTGATGCAAAATTCCGTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCAACGCTACGCTCAAACATAATGATTCTAATAAAACCTCAGGAGACTACTATGCCTGAAGCAACACCTTTTCAGGTGATGATTGTGGATGATCATCCACTTATGCGACGCGGTGTTCGTCAGTTACTGGAGCTTGATCCTGGCT  >  NC_000913/2290300‑2290593

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