Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,045,252 C→T P352S (CCG→TCG)  yihF → DUF945 domain‑containing protein YihF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,045,2520CT95.0% 52.2 / NA 20P352S (CCG→TCG) yihFDUF945 domain‑containing protein YihF
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (10/9);  minor base G (0/1);  total (10/10)
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.79e-01

GATAAGTTAACAAAAGGGAATCAGCAACTCGGTAGTGGTGAGTGGTCATTGATTGCTGAATCTATTGATCCCTCAGCAGTGCGCCAATTTATCATCCAGTATAACATTGCGATGCAGAAGCAGCTTGCTGCACACCCTGAGTTAGCAAACGATGAAGTTGCTCTGCAAGAAGTGAATGCTGCATTGTTCAAAGAGTATTTACCGTTATTACAAAAAAGTGAGCCGACCATTAAACAACCGGTAAGATGGAAGAACGCACTCGGCGAACTAAATGCCAATCTGGATATCAGTATTGCCGACCCAGCCAAATCTTCATCATCCACAAACAAAGATATCAAATCGCTCAATTTTGATGTGAAGTTACCGCTTAATGTCGTCACAGAA  >  NC_000913/4045015‑4045398
                                                                                                                                                                                                                                             |                                                                                                                                                  
gATAAGTTAACAAAAGGGAATCAGCAACTCGGTAGTGGTGAGTGGTCATTGATTGCTGAATCTATTGATCCCTCAGCAGTGCGCCAATTTATCATCCAGTATAACATTGCGATGCAGAAGCAGCTTGCTGCACACCCTGAGTTAGCAACCGATGAAGTTGCTCTGCAAGAAGTGAATGCTGCATTGTTCAAAGAGTATTTACCGTTATTACAAAAAAGTGAGCCGACCATTAAACAATCGGTAAGATGga                                                                                                                                        >  1:117345/1‑250 (MQ=255)
  tAAGTTAACAAAAGGGAATCAGCAACTCGGTAGTGGTGCGTTGTCATTGATTGCTGAATCTATTGATCCCTCAGCAGTGCGCCAATTTATCCTCCAGTATAACATTGCGATGCAGAAGCAGATTGCTGCACACCCTGAGTTAGCAAACGATGAAGTTGCTCTGCAAGAAGTGAATGCTGCATTGTTCAAAGAGTATTTACCGTAATTACAATAAAGTGAGCCGACCATTAAACAATCGGTAAGATGga                                                                                                                                        >  1:82892/1‑248 (MQ=255)
     gTTAACAAAAGGGAATCATCAACTCGTTAGTGTTGAGTGGTCATTGATTGCTGAATCTATTGATCCCTCAGCAGCGCGCCAATTTATCATCCAGTATAACATTGCGATGCAGAAGCAGCTTGCTGCACACTCTGATTTAGTAATCGATTAAGGTGCTCTGCAAGAATTTAATGCTGCATTTTGCAAAGAGTAGTTTCCGTTATCACAAAAAAGTGAGCCGACCATTAAACAATCGGTAAGATGga                                                                                                                                        <  1:102366/245‑1 (MQ=255)
                        cAACTCGGTAGTGGTGAGTGGTCATTGATTGCTGACTCACTAGATCCCTCCGCAGTGCGCCAATTTATCATCCAGTATAACATTGCGATGCAGAAGCAGCTTGCTGCACACCCTGAGTTCGCAAACGCTGAAGTTGCTCTGCAAGAAGTGAATGCTGCATTGTTCAAAGAGTATTTACCGTTATTACAAAAAAGTGAGCCGACCATTAAACAATCGGTAAGATGGAAGAACGCACTCGGCGAACTAAATGc                                                                                                               >  1:191371/1‑251 (MQ=255)
                         aaCTCGGTAGTGGTGAGTGGTCATTGATTGCTGAATCTATTGATCCCTCAGCAGTGCGCCAATTTATCATCCAGTATAACATTGCGATGCAGAAGCAGCTTGCTGCACACCCTGAGTTAGCAAACGATGAAGTTGCTCTGCAAGAAGTGAATGCTGCATTGTTCAAAGAGTATTTACCGTTATTACAAAAAAGTGAGCCGACCATTAAACAATCGGTAAGATGGAAGAACGCACTCGGCGAACTAAATGcc                                                                                                              <  1:92353/251‑1 (MQ=255)
                                                      gCTGAATCTATTGATCCCTCAGCAGTGCGCCAATTTATCATCCAGTATAACATTGCGATGCAGAAGCAGCTTGCTGCACACCCTGAGTTAGCAAACGATGAAGTTGCTCTGCAAGAAGTGAATGCTGCATTGTTCAAAGAGTATTTACCGTTATTACAAAAAAGTGAGCCGACCATTAAACAATCGGTAAGATGGAAGAACGCACTCGGCGAACTAAATGCCAATCTGGATATCAGTATTGCCGACCCAGc                                                                                 <  1:172042/251‑1 (MQ=255)
                                                                                   ccAATTTATCATCCAGTATAACATTGCGATGCAGAAGCAGCTTGCTGCACACCCTGAGTTAGTAAACGATGCAGTTGCTCTGCAAGAAGTGAATGCTGCCTTGTTCCAAGAGTATTTACCGTTCTTCCAAAAAAGTGAGCCGACCATTAAACAATCGGTAAGATGGCAGAACGCACTCGGCGCACTAAATTCCAATCTGGATATCAGTATTGCCGCCCCAGCCAAATCTTCATCATCCACAAACAAAGat                                                     >  1:160451/1‑250 (MQ=255)
                                                                                                            gCGATGCAGAAGCAGCTTGCTGCACANCCTGAGNNAGNAAACGATGAAGTTGCTCTGNAANNNGNNNANGCTGCNTTNNTCAAAGAGTNTTTACCGTTATTACAAAAAAGTGAGCCGACCATTAAACAATCGGTAAGATGGANGAACGCACTNNGCGAACTAAATGCCAATCTGGNTATNNGNNTNGCCNAcn                                                                                     >  2:46708/1‑192 (MQ=255)
                                                                                                            gCGATGCAGAAGCAGCTTGCTGCACACCCTGAGTTAGCAAACGATGAAGTTGCTCTGCAAGAAGTGAATGCTGCATTGTTCAAAGAGTATTTACCGTTATTACAAAAAAGTGAGCCGACCATTAAACAATCGGTAAGATGGAAGAACGCACTCGGCGAACTAAATGCCAATCTGGATATCAGTATTGCCGACCCAg                                                                                  <  1:46708/196‑1 (MQ=255)
                                                                                                                                 gcACACCCTGAGTTAGCAAACGATGAAGTTGCTCTGCAAGAAGTGAATGCTGCATTGTTCAAAGAGTATTTACCGTTATTACAAAAAAGTGAGCCGACCATTAAACAATCGGTAAGATGGACGAACGCACTCGGCGAACTAAATGCCAATCTGGATATCAGTATTGCCGACCCAGCCAAATCTTCATCATCCACAAACAAAGATATCAAATCGCTCAATTTTGATGTGAAGTTACCGCTTAATGTCGTcac      >  1:178838/1‑251 (MQ=255)
                                                                                                                                 gcACACCCTGAGTTAGCAAACGATGAAGTTGCTCTGCAAGAAGTGAATGCTGCATTGTTCAAAGAGTATTTACCGTTATTACAAAAAAGTGAGCCGACCATTAAACAATCGGTAAGATGGAAGAACGCACTCGGCGAACTAAATGCCAATCTGGATATCAGTATTGCCGACCCAGCCAAATCTTCATCATCCACAAACAAAGATATCAAATCGCTCAATTTTGATGTGAAGTTACCGCTTAATGTCGTcac      <  1:14616/251‑1 (MQ=255)
                                                                                                                                      cccTGAGTTAGCAAACGATGAAGTTGCTCTGCAAGAAGTGAATGCTGCATTGTTCAAAGAGTATTTACCGTTATTACAAAAAAGTGAGCCGACCATTAAACAATCGGTAAGATGGAAGAACGCACTCGGCGAACTAAATGCCAATCTGGATATCAGTATTGCCGACCCAGCCAAATCTTCATCATCCACAAACAAAGATATCAAATCGCTCAATTTTGATGTGAAGTTACCGCTTAATGTCGTCACAGaa  >  1:125741/1‑250 (MQ=255)
                                                                                                                                      cccTGAGTTAGCAAACGATGAAGCTGCTCTGCAAGAAGTGACTGCTGCATTGTTGCAAGAGTCTTTACCGTTATTACAAAAAAGTGAGCCGCCCATTAAACACTCGGTAAGATGGAAGAACGCACTCGGCGAACTAAATGCCACTCTGGATATCAGTATTGCCGACCCAGCCAACTCTtca                                                                       >  1:89087/1‑181 (MQ=255)
                                                                                                                                                      gATGAAGTTGCTCTTCAAGAAGTGAATGCTGCATTGTCCAAAGAGTATTTACCGTTATTACAAAAAAGTGAGCCGACCATTAAACAAGCGGTAAGATGGAAGAACGCACTCGGCGAACTAAATGc                                                                                                               <  1:211122/125‑1 (MQ=255)
                                                                                                                                                                     cAAGANNTGNNNGCNNCATTGTTCAAAGAGTATTTACCGNTATTACAAAAAAGTGAGCCGACCATTAAANANTNGGTANNATGNAAGAACGCACNCGGCGAACTAAATGCCAATCTGGATATCAGTATTGcc                                                                                         >  2:164091/1‑132 (MQ=255)
                                                                                                                                                                     cAAGAATTGAATGCTGCATTGTTCAAAGAGTATTTACCGTTATTACAAAAAAGTGAGCCGACCATTAAACAATCGGTAAGATGGAAGAACGCACGCGGCGAACTAAATGCCAATCTGGATATCAGTATTGcc                                                                                         <  1:164091/132‑1 (MQ=255)
                                                                                                                                                                                              aaNAGTATNTACCNTTATTACAAAAAAGTNAGCCGACCATTAAACAATCNNTNAGAGGnaagaa                                                                                                                                    <  2:82892/64‑1 (MQ=255)
                                                                                                                                                                                                            ttattaCAAAAAAGTGAGCCGACCATTAAACAATCGGTAAGATGGAAGAACGCACTCGGCGAACTAAATGCCAATCTGGATATCAGTATTg                                                                                           <  1:45703/91‑1 (MQ=255)
                                                                                                                                                                                                                         gTGAGCCGACCATTAAACAATCGGTANGATGGANNAANGCACTCGGCGAACTAAATGNCANNNTNNNTNTCAGTNTTNNCGACCCAGCNAAATCTTCATCATCCACAAACAAAGATATCAAATCGCTCAATTTTGATGTGAANTTACCGCTTNNTg             >  2:46047/1‑156 (MQ=255)
                                                                                                                                                                                                                         gTGAGCCGACCATTAAACAATCGGTAAGATGGAAGAACGCACTCGGCGAACTAAATGCCAATCTGGATATCAGTATTGCCGACCCAGCCAAATCTTCATCATCCACAAACAAAGATATCAAATCGCTCAATTTTGATGTGAAGTTACCGCTTAATg             <  1:46047/156‑1 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                  
GATAAGTTAACAAAAGGGAATCAGCAACTCGGTAGTGGTGAGTGGTCATTGATTGCTGAATCTATTGATCCCTCAGCAGTGCGCCAATTTATCATCCAGTATAACATTGCGATGCAGAAGCAGCTTGCTGCACACCCTGAGTTAGCAAACGATGAAGTTGCTCTGCAAGAAGTGAATGCTGCATTGTTCAAAGAGTATTTACCGTTATTACAAAAAAGTGAGCCGACCATTAAACAACCGGTAAGATGGAAGAACGCACTCGGCGAACTAAATGCCAATCTGGATATCAGTATTGCCGACCCAGCCAAATCTTCATCATCCACAAACAAAGATATCAAATCGCTCAATTTTGATGTGAAGTTACCGCTTAATGTCGTCACAGAA  >  NC_000913/4045015‑4045398

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