Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,488,579 (C)7→8 coding (1485/1503 nt) yjgR ← DUF853 domain‑containing protein YjgR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,488,5721.C97.8% 164.8 / ‑3.2 46coding (1492/1503 nt)yjgRDUF853 domain‑containing protein YjgR
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (27/18);  total (27/19)
Fisher's exact test for biased strand distribution p-value = 4.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

TGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCTCCCCCCCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCTTTTGGCCATTGTTTGCACCACACCATCTTTCTTCCCGCCGCGTGGTCCGGTAGTGCCGAACA  >  NC_000913/4488442‑4488697
                                                                                                                                    |                                                                                                                             
tGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCTcccccccca                                                                                                                      >  1:650835/1‑138 (MQ=255)
   tAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCTccccccccagca                                                                                                                   <  2:372786/140‑6 (MQ=255)
    aGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCGATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCTCCCCCCCCAGCaa                                                                                                                   >  2:93357/1‑140 (MQ=255)
    aGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCTCCCCCCCCAGCaa                                                                                                                   >  1:266581/1‑140 (MQ=255)
    aGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCTCCCCCCCCAGCaa                                                                                                                   >  2:300343/1‑140 (MQ=255)
           tgtgCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCTCCCCCCCCAGCAAACTcccc                                                                                                            >  1:678059/1‑140 (MQ=255)
              gCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCTCCCCCCCCAGCAAACTCCCCAAc                                                                                                         <  2:266581/140‑1 (MQ=255)
              gCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCTCCCCCCCCAGCAAACTCCCCAAc                                                                                                         >  1:464655/1‑140 (MQ=255)
                  cGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCTccccccccagca                                                                                                                   <  1:736766/125‑6 (MQ=255)
                  cGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCTccccccccagca                                                                                                                   >  2:736766/1‑120 (MQ=255)
                   gTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCTCCCCCCCCAGCAAACTCCCCAACATccc                                                                                                    <  1:745231/140‑1 (MQ=255)
                      cATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCTCCCCCCCCAGCAAACTCCCCAACATCCcgcg                                                                                                 <  2:347946/140‑1 (MQ=255)
                                               aaTATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCTCCCCCCCCAGCAAACTCCCCAACATCCCGCGTACAATCTTATTCGTCACTTGGCGa                                                                        >  1:490671/1‑140 (MQ=255)
                                                      tGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCTCCCCCCCCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGgcgc                                                                 >  1:449065/1‑140 (MQ=255)
                                                      tGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCTCCCCCCCCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGgcgc                                                                 >  2:427108/1‑140 (MQ=255)
                                                      tGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCTCCCCCCCCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGgcgc                                                                 >  2:518969/1‑140 (MQ=255)
                                                      tGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCTCCCCCCCCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGgcgc                                                                 >  1:511427/1‑140 (MQ=255)
                                                      tGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCTCCCCCCCCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGgcgc                                                                 >  2:715865/1‑140 (MQ=255)
                                                        cAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCTCCCCCCCCAGCAAACTCCCCAACATCCCGAGTACAATCTGATTCGTCACTTGGCGAGCGGCGCtt                                                               >  2:635930/1‑140 (MQ=255)
                                                         aGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCTCCCCCCCCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCttt                                                              >  1:630547/1‑140 (MQ=255)
                                                           tggggtgggCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCTCCCCCCCCAGCAAACTCCCCAACATCCCGCGTACAATCTGCTTCGTCACTTGGCGAGCGGCGCTTTTg                                                            >  2:398235/1‑140 (MQ=255)
                                                           tggggtgggCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCTCCCCCCCCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCTTTTg                                                            >  1:656168/1‑140 (MQ=255)
                                                           tggggtgggCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTATTCTCCCCCCCCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCTTTTg                                                            >  1:815894/1‑140 (MQ=255)
                                                             gggtgggCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCTCCCCCCCCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCTTTTGGc                                                          >  1:356693/1‑140 (MQ=255)
                                                             gggtgggCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCTCCCCCCCCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCTTTTGGc                                                          >  1:346936/1‑140 (MQ=255)
                                                             gggtgggCATGATAAGACGCGCCAGCGTCGCATCAGACAAAGACAAATCTCGCCTGACTTACCTTCTTCTCCCCCCCCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCTTTTGGc                                                          >  1:499939/1‑140 (MQ=255)
                                                              ggtgggCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCTCCCCCCCCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCTTTTGGcc                                                         <  2:630547/140‑1 (MQ=255)
                                                               gtgggCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCTCCCCCCCCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCTTTTGGCCa                                                        <  1:427108/140‑1 (MQ=255)
                                                               gtgggCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCTCCCCCCCCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCTTTTGGCCa                                                        >  2:10435/1‑140 (MQ=255)
                                                                tgggCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCTCCCCCCCCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCTTTTGGcc                                                         <  2:290495/138‑1 (MQ=255)
                                                                tgggCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCTCCCCCCCCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCTTTTGGcc                                                         >  1:290495/1‑138 (MQ=255)
                                                                   gCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCTCCCCCCCCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCTTTTGGCCATTGt                                                    <  1:635930/140‑1 (MQ=255)
                                                                       gATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCTCCCCCCCCAGCAAACTCCCCAACATCCCGCGTACAATCTGGTTCGTCACTTGGCGAGCGGCGCTTTTGGCCATTGt                                                    >  1:488763/1‑136 (MQ=255)
                                                                       gATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCTCCCCCCCCAGCAAACTCCCCAACATCCCGCGTACAATCTGGTTCGTCACTTGGCGAGCGGCGCTTTTGGCCATTGt                                                    <  2:488763/136‑1 (MQ=255)
                                                                       gATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCTCCCCCCCCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCTTTTGGCCATTGTTTGc                                                <  1:518969/140‑1 (MQ=255)
                                                                        aTAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCTCCCCCCCCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCTTTTGGCCATTGTTTGca                                               >  1:209781/1‑140 (MQ=255)
                                                                           aGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCTCCCCCCCCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCTTTTGGCCATTGTTTGcacca                                            <  2:656168/140‑1 (MQ=255)
                                                                                 gcCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCTCCCCCCCCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCTTTTGGCCATTGTTTGCACCACACCAt                                      >  1:230468/1‑140 (MQ=255)
                                                                                  ccAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCTCCCCCCCCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCACTTTTGGCCATTGTTTGCACCACACCAtc                                     >  2:267458/1‑140 (MQ=255)
                                                                                      cGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCTCCCCCCCCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCTTTTGGCCATTGtt                                                   >  2:865285/1‑122 (MQ=255)
                                                                                      cGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCTCCCCCCCCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCTTTTGGCCATTGtt                                                   <  1:865285/122‑1 (MQ=255)
                                                                                                 gCAAAGACAAATCTCGCCTGACTTACCTTCTTCTCCCCCCCCAGCAAACTCCCCAACATCCCGCGTACAATCTTATTCGTCACTTGGCGAGCGGCGCTTTTGGCCATTGTTTGCACCACACCATCTTTCTTCCCGCcgcg                      <  2:490671/140‑1 (MQ=255)
                                                                                                 gCAAAGACAAATCTCGCCTGACTTACCTTCTTCTCCCCCCCCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCTTTTGGCCATTGTTTGCACCACACCATCTTTCTTCCCGCcgcg                      <  1:535699/140‑1 (MQ=255)
                                                                                                 gCAAAGACAAATCTCGCCTGACTTACCTTCTTCTCCCCCCCCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCACTTTTGGCCATTGTTTGCACCACACCATCTTTCTTCCCGCcgcg                      <  1:267458/140‑1 (MQ=255)
                                                                                                    aaGACAAATCTCGCCTGACTTCCCTTCTTCCCCCCCCCCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCTTTTGGCCATTGTTTGCACCACCCCATCTTTCTTCCCGCCGCGTgg                   <  1:321905/140‑1 (MQ=255)
                                                                                                                gCCTGACTTACCTTCTTCTCCCCCCCCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCTTTTGGCCATTGTTTGCACCACACCATCTTTCTTCCCGCCGCGTGGTCCGGTAGTGcc       <  2:449065/140‑1 (MQ=255)
                                                                                                                gCCTGACTTACCTTCTTCTCCCCCCCCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCTTTTGGCCATTGTTTGCACCACACCATCTTTCTTCCCGCCGCGTGGTCCGGTAGTGcc       >  2:345489/1‑140 (MQ=255)
                                                                                                                  cTGACTTACCTTCTTCTCCCCCCCCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCTTTTGGCCATTGTTTGCACCACACCATCTTTCTTCCCGCCGCGTGGTCCGGTAGTGCCGa     <  1:634003/140‑1 (MQ=255)
                                                                                                                   tGACTTACCTTCTTCTCCCCCCCCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCACTTTTGGCCATTGTTTGCACCACACCATCTTTCTTCCCGCCGCGTGGTCCGGTAGTGCCGaa    <  1:393665/140‑1 (MQ=255)
                                                                                                                     aCTTACCTTCTTCTCCCCCCCCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCTTTTGGCCATTGTTTGCACCACACCATCTTTCTTCCCGCCGCGTGGTCCGGTAGTGCCGAACa  <  1:715865/140‑1 (MQ=255)
                                                                                                                                    |                                                                                                                             
TGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCTCCCCCCCAGCAAACTCCCCAACATCCCGCGTACAATCTGATTCGTCACTTGGCGAGCGGCGCTTTTGGCCATTGTTTGCACCACACCATCTTTCTTCCCGCCGCGTGGTCCGGTAGTGCCGAACA  >  NC_000913/4488442‑4488697

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