Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 184,253 C→A 12.9% intergenic (‑158/+4) yaeH ← / ← yaeI UPF0325 family protein/phosphodiesterase with model substrate bis‑pNPP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913184,2530CA12.9% 84.5 / 3.7 44intergenic (‑158/+4)yaeH/yaeIUPF0325 family protein/phosphodiesterase with model substrate bis‑pNPP
Reads supporting (aligned to +/- strand):  ref base C (23/15);  new base A (2/4);  total (25/19)
Fisher's exact test for biased strand distribution p-value = 3.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.08e-01

CAAAATGAGTCAGTAAATTGTGCTTATTTTAGCATTTGGCCTGGCCCGCCCGGCTCGACTGTGATCGCAGAAAGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATCTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGTTGTGTAGATGTGTCTTTCGCCAAAGGCATTTAATCCGGCGACGTAACGTTTATCTTCGACAGGG  >  NC_000913/184106‑184395
                                                                                                                                                   |                                                                                                                                              
cAAAATGAGTCAGTAAATTGTCCTTATTTTAGCATTTGGCCTGGCCCGCCCGGCTCGACTGTGATCGCAGAAAGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACg                                                                                                                                               <  1:960758/149‑1 (MQ=255)
  aaaTGAGTCAGTAAATTGTGCTTATTTTAGCATTTGGCCTGGCACGCCCGGCTCGAATGTGATCGCAGAAAGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGAt                                                                                                                                             <  1:598912/149‑1 (MQ=255)
   aaTGAGTCAGTAAATTGTGCTTATTTTAGCATTTGGCCTGGCCCGCCCGGCTCGACTGGGATCGCAGAAAGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAACTTACCTGCTTTATTCTGGTGATAAAATTCACGATc                                                                                                                                            <  2:404915/149‑1 (MQ=255)
        gTCATTAAATAGTGCTTATTTTAGCATTTGGCCTGGCCAGCCCGGCTCGACAGTGATCGCAGAAAGACTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATCTacac                                                                                                                                       <  1:84881/149‑1 (MQ=255)
            ttAAATTGTGCTTATTTTAGCATTTGGCCTGGCCCGCCCGGCTCGACTGTGATCGCAGAAAGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATCTACACCAGt                                                                                                                                   >  2:602070/2‑149 (MQ=255)
                     gCTTATTTTAGCATTTGGCCTGGCCCGCCCGGCGCGACTGTGATCGCAGAAAGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTTTTTTTTTTATTGTTTTGTTTAAATTCTAGATCTACACCAGTTCCAGCATc                                                                                                                          <  2:309308/149‑1 (MQ=255)
                        tcttttAGCATTTGGCCTGGCCCGCCCGGCGCGACTGTGATCGCAGAAAGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTTTTTTTTTTATTGTTTTGTTTAAATTCTAGATCTACACCAGTTCCAGCATCGtt                                                                                                                       <  2:163556/147‑1 (MQ=255)
                              aGCATTTGGCCTGGCCCGCCCGGCTCGACTGTGATCGCAGAAAGCCTTTTTCTTCCTGACTTCTATCATGGCGTAAAAATCAAAAAATTACCAGCTTTATTCTGGTGATAAAATTCACGATCTACACCAGTTCCAGCATCGGTACTTcc                                                                                                                 >  2:306414/1‑149 (MQ=255)
                                      gcctggccCGCCCGGCTCGACTGTGATCGCAGAAAGCCTTTTTCTTCCAGACTTCTATCATGGCGTTAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATCTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCa                                                                                                         >  1:471888/1‑149 (MQ=255)
                                           gcccgcccGTCTCGACTGTGATCGCAGAAAGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATCTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCa                                                                                                    <  1:488357/149‑1 (MQ=255)
                                           gcccgcccGGCTCGACTGTGATCGCAGAAAGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATCTACACGAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCa                                                                                                    <  1:274530/149‑1 (MQ=255)
                                              cgcccgGCTCGACTGTGATCGCAGAAAGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATCTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGAc                                                                                                 <  1:659495/149‑1 (MQ=255)
                                                cccgGCTCGACTGTGATCGCAGAAAGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATCTACCCCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGc                                                                                               <  1:594695/149‑1 (MQ=255)
                                                cccgGCTCGACTGTGATCGCAGAAAGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATCTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGc                                                                                               <  1:846489/149‑1 (MQ=255)
                                                 ccgGCTCGACTGTGATCGCAGAAAGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATCACCACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGCCGCa                                                                                              >  2:259570/1‑149 (MQ=255)
                                                        gACTGTGATCGCAGAAAGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATCTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATa                                                                                       >  1:448421/1‑149 (MQ=255)
                                                        gACTGTGATCGCAGAAAGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATCTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATa                                                                                       >  2:267600/1‑149 (MQ=255)
                                                           tgtgATCGCAGAAAGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATCTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACaa                                                                                    >  1:415009/1‑149 (MQ=255)
                                                             tgATCGCAGAAAGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTATAGTGATAAAATTCAAGATCTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACCCAAACCATAAAAAc                                                                                  >  2:600369/1‑149 (MQ=255)
                                                               aTCGCAGAAAGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATCTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTg                                                                                >  2:856436/1‑149 (MQ=255)
                                                                tCGCAGAAAGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATCTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAACCCATACAAACTGc                                                                               >  2:964302/1‑149 (MQ=255)
                                                                tCGCAGAAAGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATCTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGc                                                                               >  1:19142/1‑149 (MQ=255)
                                                                tCGCAGAAAGCCTTTTACTTCCAGACTACTATCAAGGCGTAAAAATCAAAAAATTACCTGCTTTAATCTGGTGATAAAATTCAAGATCTACACCAGTTACAGCATCGTTACTTCCGTGCGGCAGTTCACACGCAAACCATACAAAATGc                                                                               >  1:412712/1‑149 (MQ=255)
                                                                            tttttCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATCTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACggg                                                                   >  2:174627/1‑149 (MQ=255)
                                                                                  ttcCAGACTTCTATCATTGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATCTACACCGGGTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACTCCACGGGTTgtgt                                                             <  1:226178/149‑1 (MQ=255)
                                                                                  ttcCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATCTACACCAGATCCAGCATCGTTACTTCCGGCCGGCAGTTCAGACGCAAACCATACAAACAGCACACGCCACGGGTTgagt                                                             >  2:683219/1‑149 (MQ=255)
                                                                                    ccAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATCTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGTTGTGTAg                                                           >  2:1005517/1‑149 (MQ=255)
                                                                                    ccAGACTTCTATCATGCCGTAAAAATCATAAAATTACCTGCTTTATTCTGGTGAAAAAATACACGAACTACACCAGTTCCAGCCTCGTTACTTCCGGCCGGCAGTTCAGACGCAAACCATACACACTGCCCACGCCACGGGGTGTGTAg                                                           <  1:328539/149‑1 (MQ=255)
                                                                                     cAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATCTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGTTGTGTAGa                                                          <  1:856436/149‑1 (MQ=255)
                                                                                        aCTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATCTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGTTGTGTAGAtgt                                                       <  1:267600/149‑1 (MQ=255)
                                                                                                  tGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATCTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATCCAAACTGCCCACGCCACGGGTTGGGTAGTTGTGTCTTTCGcc                                             >  2:866098/1‑149 (MQ=255)
                                                                                                          aaaTCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATCTACACCAGTTCCAGCATCGTTACTTCCTGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGTTGTGTAGATGTGTCTTTCGCCAAAGGAAt                                     >  1:106998/1‑149 (MQ=255)
                                                                                                               aaaaaaTTACCTGCTTTATTCTGGTGATAAAATTAACGATCTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCTCACCATACAATCTGCCCACGCCACGGGTTGTGTAGATGTGTCTTTCGCCAAAGGCATTTact                                >  1:780632/1‑147 (MQ=255)
                                                                                                                 aaaaTTACCTGCTTTATTCTGGTGATAAAATTCACGATCTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGTTGTGTAGATGTGTCTTTCGCCAAAGGCATTTAATcc                              >  1:873879/1‑149 (MQ=255)
                                                                                                                      tACCTGCTTTATTCTGGTGATAAAATTCACGATCTACAACAGTTCCAGCATCGTTAGTTCCGGGCGGCAGTTCAGCCGCAAACCATACAAACTGCCCACGCCACGGGTTGTGTAGATGTGTCTTTCGCCAAAGGCATTTTATCCGGCGa                         <  1:683219/149‑1 (MQ=255)
                                                                                                                         cTGCTTTATTCTGGTGATAAAATTCACGATCTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCTGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGTTGGGTAGATGTGTCTTTCGCCAAAGGCATTTAATCCGGCGACGt                      >  2:351058/1‑149 (MQ=255)
                                                                                                                           gCTTTATTCTGGTGATAAAATTCACGATCTACACCAGTTCCAGCATCGTTACTTGCGGGCGGCAGTTCAGACGCAAACCATTCAACCTGCCCACGCCACGGGTTGTGTAGATGTGTCTTTCGCCAAAGGCTTTTAATCCGGCGACGTaa                    >  2:618728/1‑149 (MQ=255)
                                                                                                                                aTTCTGGTGATAAAATTCACGATCTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGTTGTGTAGATGTGTCGTTAGCCAAAGGCATTTAATCCCGCGACGTAACGttt               >  2:900755/1‑149 (MQ=255)
                                                                                                                                    tGGTGATAAAATTCACGATCTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACACCACGGGTTGTGTAGATGTGTCTTTCGCCAAAGGCATGTAATCCGGCGACGTATCGTTTATCt           >  1:979445/1‑149 (MQ=255)
                                                                                                                                     ggTGATAAAATTCACGATCTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGTTGTGTAGATGTGTCTTTCGCCAAAGGCATTTAATCCGGCGACGTAACGTTTATCtt          <  1:1026719/149‑1 (MQ=255)
                                                                                                                                       tGATAAAATTCACGATCTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAACCCATACAAACTGCCCCCGCCAAGGGTTGTGTAGATGTGTCTTTCGCCAAAGGCATTTAATACGGCGACGTAACGTTTATCTTCg        >  2:349084/1‑149 (MQ=255)
                                                                                                                                       tGATAAAATTCACGATCTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGTTGTGTAGATGTGTCTTTCGCCAAAGGCATTTAATCCGGCGACGTAACGTTTATCTTCg        >  1:16306/1‑149 (MQ=255)
                                                                                                                                            aaaTTCTAGATCTACACCAGTTCGAGCTTCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGTTGTGTAAATGTGTCTTTCGCCAAAGGCATTTAATCCGGCGACGTAACGTTTATCTTCGACAgg   <  2:412712/149‑1 (MQ=255)
                                                                                                                                             aaTTCTAGATCTACACATGTTCCAGCATCGTTACTTACGGGCGGCAGTTCAGACGCAAACCATACAAACTGACCACGCCACGGGTTGTGTAGATGTGTCTTTCGCCAAAGGAATTTAATCCGGCGACGTAACGTTTATCTTCGACAggg  <  2:902831/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
CAAAATGAGTCAGTAAATTGTGCTTATTTTAGCATTTGGCCTGGCCCGCCCGGCTCGACTGTGATCGCAGAAAGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATCTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGTTGTGTAGATGTGTCTTTCGCCAAAGGCATTTAATCCGGCGACGTAACGTTTATCTTCGACAGGG  >  NC_000913/184106‑184395

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑