Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 635,361 G→A 14.8% Q6* (CAA→TAA) 
C403C (TGC→TGT
ybdM ←
ybdN ←
ParB‑like nuclease domain‑containing protein YbdM
putative PAPS reductase/DUF3440 domain‑containing protein YbdN

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913635,3610GA14.8% 54.9 / 3.9 27Q6* (CAA→TAA) 
C403C (TGC→TGT
ybdM
ybdN
ParB‑like nuclease domain‑containing protein YbdM
putative PAPS reductase/DUF3440 domain‑containing protein YbdN
Reads supporting (aligned to +/- strand):  ref base G (11/12);  new base A (2/2);  total (13/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.08e-01

CACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACGTTTTCCTTTTATCCGTTCGTTATAACGCTGATAGTTCTTCGATTTTGTGGGGCTAAATGATAATGCCCGACACCAGTAATCGTTATTGAGTAAAACTTTGCATATCCGCCGCCA  >  NC_000913/635222‑635498
                                                                                                                                           |                                                                                                                                         
cacaGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCAt                                                                                                                                        >  2:262669/1‑143 (MQ=255)
 acaGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATa                                                                                                                                       >  1:204863/1‑143 (MQ=255)
       cGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGTCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATcc                                                                                                                                 <  1:445971/143‑1 (MQ=255)
       cGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAGATTGCGTTACATCCTGCGTTACTCGTTGTTGCATAGTacgc                                                                                                                                 >  2:106409/1‑140 (MQ=255)
         cAATCCACGGTCTGCTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATcccc                                                                                                                               <  1:208272/143‑1 (MQ=255)
              cACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTc                                                                                                                          <  1:390048/143‑1 (MQ=255)
              cACGGGCTCTTCACGAAAGGGGCGAACCTGATGGATCGCCAGGCGGATTTCATTAAGCGCTTTGATGCGATCATCATCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTc                                                                                                                          <  2:8562/143‑1 (MQ=255)
                       ttCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTAACGttt                                                                                                                 >  1:381797/1‑143 (MQ=255)
                       ttCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCAGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATTCCCCATTCTTGACGttt                                                                                                                 >  1:381977/1‑143 (MQ=255)
                           cGAGAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACGTTTTCCt                                                                                                             <  1:437736/143‑1 (MQ=255)
                                     tAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACGTTTTCCTTTTATCcgtt                                                                                                   >  1:289288/1‑143 (MQ=255)
                                                aTCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACGTTTTCCTTTTATCCGTTCGTTATAACGc                                                                                        >  1:265625/1‑143 (MQ=255)
                                                        gCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCGGGCGCGCATGAGACCAAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACGTTTTCCTTTTATCCGTTCGTTATAACGCTGATAGtt                                                                                <  1:262669/143‑1 (MQ=255)
                                                               tCATGACTCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACGTTTTCCTTTTATCCGTTCGTTATAACGCTGATAGTTCTTCGAt                                                                         <  2:246265/143‑1 (MQ=255)
                                                                cATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACGTTTTCCTTTTATCCGTTCGTTATAACGCTGATAGTTCTTCGAtt                                                                        >  1:559586/1‑143 (MQ=255)
                                                                         cTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGGTGTTGCATAGCATCCCCCATTCCTGACGCTTCCCTTTTATCCGTTCGTTATAACGCTGACAGTTCTCCGATCTTGCgggcc                                                               >  2:146885/1‑141 (MQ=255)
                                                                                 cGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACGTTTTCCTTTTATCCGTTCGTTATAACGCTGATAGTTCTTCGATTTTGTGGGGCTAAATGAt                                                       <  2:433962/143‑1 (MQ=255)
                                                                                          tCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACGTTTTCCTTTTCTCCGTTCGTTATAACCCCGATAGTTCTTCGATTTTGTGGGGCTAAGTGATACTCCCCGa                                              >  2:525327/1‑143 (MQ=255)
                                                                                                      gCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACGTTTTCCTTTTATCCGTTCGTTATAACGCTGATAGTTCTTCGATTTTGTGGGGCTACATGATAATGCCCGACACCAGTACTCg                                  >  2:269061/1‑143 (MQ=255)
                                                                                                            aaTTGCGTTAAGTCCTGCGTCAATCGTTGTTGCATAGTATCCCCCATTCCTGACGTTTTCCTTTTATCCGTTCGTTATAACGCTGATAGTTCTTCGATTTTGTGGGGCTAAATGATAATGCCCGACACCAGTAATCGTTATTg                            <  2:188098/143‑1 (MQ=255)
                                                                                                                 cGTTAAATCCTGCGTTAATCGTTGTTACATAGTATCCCCCATTCCTGACGTTTTCCTTTTATCCGTTCGTTATAACGCTGATAGTTCTTCGATTTTGTGGGGCTAAATGATAATGcc                                                 <  1:318382/117‑1 (MQ=255)
                                                                                                                 cGTTAAATCCTGCGTTAATCGTTGTTACATAGTATCCCCCATTCCTGACGTTTTCCTTTTATCCGTTCGTTATAACGCTGATAGTTCTTCGATTTTGTGGGGCTAAATGATAATGcc                                                 >  2:318382/1‑117 (MQ=255)
                                                                                                                 cGTTAAATCCTGCGTTAATCGTTGTTACATAGTATCCCCCATTCCTGACGTTTTCCTTTTATCCGTTCGTTATAACGCTGATAGTTCTTCGATTTTGTGGGGCTAAATGATAATGcc                                                 >  2:314233/1‑117 (MQ=255)
                                                                                                                 cGTTAAATCCTGCGTTAATCGTTGTTACATAGTATCCCCCATTCCTGACGTTTTCCTTTTATCCGTTCGTTATAACGCTGATAGTTCTTCGATTTTGTGGGGCTAAATGATAATGcc                                                 <  1:314233/117‑1 (MQ=255)
                                                                                                                   ttAAATCCTGCGTTAATCGTTGTTGCATAGTATCCACCATAACTAACGTTTGTCTTTTATCCGTTCGTTATAACGCTGATAGTTCTTCGATTTTGTGGGGCTAAATGATAATGCCCGACACCAGTAATCGTTATTGAGTaaaa                     <  1:106409/143‑1 (MQ=255)
                                                                                                                       agcccgCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACGTTTTCCTTTTATCCGTTCGTTATAACGCTGATAGTTCTTCGATTTTGTGGGGCTAAATGATAATGCCCGACACCAGTAATCGTTATTGAGTAAAACttt                 <  1:225213/138‑1 (MQ=255)
                                                                                                                                  aTCGTTGTTGCATAGTATCCCCCATTCCTGACGTTTTCCTTTTATCCGTTCGTTATAACGCTGATAGTTCTTCGATTTTGTGGGGCTAAATGATAATGCCCGACACCAGTAATCGTTATTGAGTAAAACTTTGCATATccgcc      <  2:265625/143‑1 (MQ=255)
                                                                                                                                      ttgttgCATAGTATCCCCCATTCCTGACGTTTTCCTTTTATCCGTTCGTTATAACGCTGATAGTTCTTCGATTTTGTGGGGCTAACTGATAATGCCCGACACCAGTCATCGTTATTGAGTAAAACTTTGCATATCCGCCCCCa  >  1:65566/1‑143 (MQ=255)
                                                                                                                                           |                                                                                                                                         
CACAGAACGCAATCCACTGGCTCTTCACGAAAGGGGCTAACCTGATGGATCGCCATGCGGATTTCATTAATCGCTTTGATGCGATCATCTTCCGGCAAGCTGGCGAGAAATTGCGTTAAATCCTGCGTTAATCGTTGTTGCATAGTATCCCCCATTCCTGACGTTTTCCTTTTATCCGTTCGTTATAACGCTGATAGTTCTTCGATTTTGTGGGGCTAAATGATAATGCCCGACACCAGTAATCGTTATTGAGTAAAACTTTGCATATCCGCCGCCA  >  NC_000913/635222‑635498

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