Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,616,073 T→C E331E (GAA→GAG)  yjjI ← DUF3029 domain‑containing protein YjjI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,616,0730TC90.9% 86.4 / ‑3.2 33E331E (GAA→GAG) yjjIDUF3029 domain‑containing protein YjjI
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base C (11/19);  total (13/20)
Fisher's exact test for biased strand distribution p-value = 5.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

AAGCATGATGCGGGGCGACAGTTTGCAGATGGGTGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCTACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGG‑CAATCCCTTCTTTTTCACAC‑AGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGATGATGGCGATCTGCTGCTGGCAGTAGTGCGGTAGAGTGCGCGTAAAGAAGTCATCCAGCGATTCGCTGCGCTCGGCAATGGCTTTCAGGTTAAGGCGT  >  NC_000913/4615826‑4616316
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    
aaGCATGAGGCGGGGCGACAGTTTGCAGATGGGTGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGTCTGTGCTTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGG‑CAATCCCTTCTTTCTca                                                                                                                                                                                                                                                   <  1:361343/251‑1 (MQ=255)
    atgatgCGGTGCGACCGTTTGCAGATGGGTGATCGAATATGGCTCATCGCCATAAGGCAAACGCGCGCCCGGCGTGGTGCCGCTATAGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTCCGGCGGCATCTTTACCGTAGCGCGCGGACAA‑CCCTTCTTTCTcacacc                                                                                                                                                                                                                                               >  2:336162/1‑250 (MQ=255)
     tgatgCGGGGCGACAGTTTGCAGATGGGTGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGG‑CAATCCCTTCTTTCTcacac                                                                                                                                                                                                                                                >  1:324739/1‑249 (MQ=255)
                gACAGTTTGCAGATGGGTGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGG‑CAATCCCTTCTTTCTCACAC‑AGCAAGTTAACCg                                                                                                                                                                                                                                  <  1:284946/251‑1 (MQ=255)
                    gtgtGCAGATGGGTGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGAGGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGG‑CAATCCCTTCTTTCTCACAC‑AGCAAGTTAACCGCTTc                                                                                                                                                                                                                              <  1:196977/248‑1 (MQ=255)
                               ggTGCTCGGATCTGGCGCATCGCCATACGGCAAACGCGCGTCCGGCGTGGTGCCGATATCGGAACGGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGG‑CAATCCCTTCTTTCTCACAC‑AGCAAGTTAACCGCTTCCGCCAGCCCAt                                                                                                                                                                                                                   <  1:357043/251‑1 (MQ=255)
                                gTGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGG‑CAATCCCTTCTTTCTCACAC‑AGCAAGTTAACCGCTTCCGCCAGCCCAt                                                                                                                                                                                                                   >  1:336462/1‑250 (MQ=255)
                                  gATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGG‑CAATCCCTTCTTTCTCACAC‑AGCAAGTTAACCGCTTCCGCCAGccccgatc                                                                                                                                                                                                                >  1:317510/1‑246 (MQ=255)
                                     cGGATCTAGCTCATCGCCATATGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGCCTGTGCGTGTAACATGGCGCGTTTTTTCCAGACATATTTCACGGGGGCATTGGCGACAAACACCGCCAGTTGCGCGCTTATGCGATAACCCACTCCCTTGGCGGCGGCCTCTTTATCTGAGCTCGCGG‑CAATCCCTTCTTTT‑CACACCATCCAGTTTACCGCTTCCGCCAGcacccacctac                                                                                                                                                                                                             >  2:92630/1‑241 (MQ=255)
                                      ggATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATCTCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACCAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGG‑CAATCCCTTCTTTCTCACAC‑AGCAAGTTAACCGCTTCCGCCAGCCCATACATGcc                                                                                                                                                                                                            >  1:118090/1‑251 (MQ=255)
                                                                    gcgcCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGG‑CAATCCCTTCTTTCTCACAC‑AGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAACGTTCAGGGTTAAt                                                                                                                                                                              <  1:41826/251‑1 (MQ=255)
                                                                           ggggggtgCCGATAGCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGG‑CAATCCCTTCTTTCTCACAC‑AGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCct                                                                                                                                                                       <  1:360456/247‑1 (MQ=255)
                                                                                    gtatactggaatctacccccgtcggtgCGTGTAACAAGCCTCGGTTTGGCCCGCCATATTTCACAGGGTTATTGGCGCCAAACGCCGCCAGTTGCGCGCTGAGGCGATACCCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGG‑CAACCCCTGCTTGCTCACCC‑AGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAACGGTCAGTGTTAATCAGCCCTTCTTTCAcc                                                                                                                                                              <  2:323790/227‑1 (MQ=255)
                                                                                             aaCTGATCCCCGACTGTGCTTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGG‑CAATCCCTTCTTTCTCACAC‑AGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTa                                                                                                                                                     <  1:246840/251‑1 (MQ=255)
                                                                                              ccTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGG‑CAATCCCTTCTCTCTCACAC‑AGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTAt                                                                                                                                                    <  1:238226/250‑1 (MQ=255)
                                                                                                                          gcgcGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGG‑CAATCCCTTCTTTCTCACAC‑AGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATagag                                                                                                                        <  1:123406/251‑1 (MQ=255)
                                                                                                                                    gccagccaTATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATACCCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGG‑CAATCCCTTCTTTCTCACAC‑AGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTcaca                                                                                                               <  1:54288/250‑1 (MQ=255)
                                                                                                                                                  aCGGGGGTATTGGCGACAAACTCCACCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGG‑CAATCCCTTCTTTCTCACAC‑AGCAAGTTAACCGCTTCCGCCACCCCATACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCgatgat                                                                                                >  1:146170/1‑251 (MQ=255)
                                                                                                                                                    gggggTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGG‑CAATCCCTTCTTTCTCACAC‑AGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAACGTTCAGGGt                                                                                                                                                                                  >  1:285327/1‑167 (MQ=255)
                                                                                                                                                    gggggTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGG‑CAATCCCTTCTTTCTCACAC‑AGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAACGTTCAGGGt                                                                                                                                                                                  <  2:285327/167‑1 (MQ=255)
                                                                                                                                                       ggTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCCCGCGG‑CAATCCCTTCTTTCTCACAC‑AGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCCCAAATCGTTCAGGGTTAATCAGCCCTTCTTTCACCCGGAAGCTATTCTCCAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGATGATGGca                                                                                            >  1:206599/1‑249 (MQ=255)
                                                                                                                                                                gACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCTACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGG‑CAATCCCTTCTTTT‑CACAC‑AGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGATGATGGCGATctgctgctg                                                                                 >  1:89390/1‑251 (MQ=255)
                                                                                                                                                                        ccGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGG‑CAATCCCTTCTTTCTCACAC‑AGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGATGATGGCGATCTGCTGCTGGCagt                                                                            <  1:37753/249‑1 (MQ=255)
                                                                                                                                                                        ccGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGGAGCGCGCGG‑CAATCCCTTCTTTCTCACAC‑AGCAAGTTAACCGCTGCCGCCAGCCCATACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGTAATTCACACCGCGCATCGATGATGGCGATCTGCTGCTGGCagt                                                                            <  1:349147/249‑1 (MQ=255)
                                                                                                                                                                                   gtccTGAAGCGAAAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGG‑CAATCCCTTCTTTCGCACAC‑AGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATGGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGTATCGATGATGGCGATCTGCTGCTGGCAGTAGTGCGGTAGAGt                                                               <  2:123834/248‑1 (MQ=255)
                                                                                                                                                                                   gcgcTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGG‑CAATCCCTTCTTTCTCACAC‑AGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAACGATCAGGGTTAATCAGCCCTTCTTTCCACAGGAAGCTATTCTCAAGGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGATAATGGCGATCTGCTGCTAGAAGTAGTGCGGTagagg                                                               >  1:85679/1‑250 (MQ=255)
                                                                                                                                                                                       tGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGG‑CAATCCCTTCTTTCTCACAC‑AGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGATGATGGCGATCTGCTGCTGGCAGTAGTGCGGTAGAGTgcgc                                                           >  1:32587/1‑251 (MQ=255)
                                                                                                                                                                                                     ccTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGG‑CAATCCCTTCTTTCTCACAC‑AGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGATGATGGCGATCTGCTGCTGGCAGTAGTGCGGTAGAGTGCGCGTAAAGAAGTCATc                                             <  1:345396/250‑1 (MQ=255)
                                                                                                                                                                                                              gcggcgGCTTCTTTACCGTAGCGCGCGG‑CAATCCCTTCTTTCTCACAC‑AGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGATGATGGCGATCTGCTGCTGGCAGTAGTGCGGTAGAGTGCGCGTAAAGAAGTCATCCAGCGATTc                                    <  1:262450/251‑1 (MQ=255)
                                                                                                                                                                                                                          ttACCGTAGCGCGCGG‑CAATCCCTTCTTTCTCACAC‑AGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGATGATGGCGATCTGCTGCTGGCAGTAGTGCGGTAGAGTGCGCGTAAAGAAGTCATCCAGCGATTCGCTGCGCTCgg                         <  1:365466/250‑1 (MQ=255)
                                                                                                                                                                                                                                 aGCGCGCGG‑CAATCCCTTCTTTCTCACAC‑AGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGATGATGGCGATCTGCTGCTGGCAGTAGTGCGGTAGAGTGCGCGTAAAGAAGTCATCCAGCGATTCGCTGCGCTCGGCAATGGCt                 >  1:158533/1‑251 (MQ=255)
                                                                                                                                                                                                                                             aTCCCTTCTTTT‑CACAC‑AGAAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACAGTGGCACAAAACGTTCAGGGGTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGATGATGGCGATCTGCTGCTGGCAGTAGTGCGGTAGAGTGCGCGTAAAGAAGTCATCCAGCGATTCGCTGCGCTCGGCAATGGCTTTCAGGTTAAgg     <  2:89390/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 cgtcttTCTCACAC‑AGCAAGTTAACCGTTTCCGCCAGCCCAGACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGTGCATCGATGATGGCGATCTGCTGCTGGCAGTAGCGCGGTAGAGTGCGCGTAAAGAAGTCATCCAGCGATTCGCTGCGCTCGGCAATGGCTTTCAGGTTAAGGCGt  <  1:207942/249‑1 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    
AAGCATGATGCGGGGCGACAGTTTGCAGATGGGTGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCTACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGG‑CAATCCCTTCTTTTTCACAC‑AGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGATGATGGCGATCTGCTGCTGGCAGTAGTGCGGTAGAGTGCGCGTAAAGAAGTCATCCAGCGATTCGCTGCGCTCGGCAATGGCTTTCAGGTTAAGGCGT  >  NC_000913/4615826‑4616316

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