Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,545,596 T→C E181G (GAG→GGG)  yjiC ← YjfZ/YjiC family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,545,5960TC96.6% 211.4 / ‑3.4 59E181G (GAG→GGG) yjiCYjfZ/YjiC family protein
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (27/30);  total (28/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCTTCTATACCAATTATTGGGTCATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATC  >  NZ_CP009273/4545460‑4545726
                                                                                                                                        |                                                                                                                                  
ggCTTCTATACCAATTATTGGGTCATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCa                                                                                                                                  >  2:135363/1‑139 (MQ=255)
ggCTTCTATACCAATTATTGGGTCATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCa                                                                                                                                  >  2:503304/1‑139 (MQ=255)
  cTTCTATACCAATTATTGGGTCATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCAtt                                                                                                                                >  1:140556/1‑139 (MQ=255)
      tataCCAATTATTGGGTCATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTc                                                                                                                            >  2:764011/1‑139 (MQ=255)
                attGGGTCATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTg                                                                                                                  >  2:835666/1‑139 (MQ=255)
                     gTCATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCt                                                                                                               >  2:406197/1‑137 (MQ=255)
                        aTTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTgag                                                                                                            >  1:863905/1‑137 (MQ=255)
                             ccAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTgg                                                                                                     <  2:863905/139‑1 (MQ=255)
                              cAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGa                                                                                                    <  2:35626/139‑1 (MQ=255)
                                    aaTATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTg                                                                                                                               >  1:699656/1‑106 (MQ=255)
                                    aaTATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTg                                                                                                                               <  2:699656/106‑1 (MQ=255)
                                    aaTATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGtt                                                                                              <  2:412246/139‑1 (MQ=255)
                                    aaTATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGt                                                                                               <  2:421652/138‑1 (MQ=255)
                                     atatTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTg                                                                                             <  1:503304/139‑1 (MQ=255)
                                            cATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAatata                                                                                      >  2:792660/1‑139 (MQ=255)
                                              ttAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATc                                                                                    <  1:782562/139‑1 (MQ=255)
                                                         ggCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGtt                                                                         <  1:221294/139‑1 (MQ=255)
                                                               ggCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCa                                                                   <  2:140556/139‑1 (MQ=255)
                                                               ggCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCa                                                                   <  1:135363/139‑1 (MQ=255)
                                                                gCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTgag                                                                                                            <  1:10502/97‑1 (MQ=255)
                                                                gCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTgag                                                                                                            >  2:10502/1‑97 (MQ=255)
                                                                       cGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAgg                                                           <  1:835666/139‑1 (MQ=255)
                                                                                   cTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAAcc                                               <  1:771165/139‑1 (MQ=255)
                                                                                                cGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGa                                                                                                    <  1:750902/73‑1 (MQ=255)
                                                                                                cGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGa                                                                                                    >  2:750902/1‑73 (MQ=255)
                                                                                                  tCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTAt                                >  2:582046/1‑139 (MQ=255)
                                                                                                  tCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTAt                                >  2:150718/1‑139 (MQ=255)
                                                                                                    aTTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTAttt                              >  1:452721/1‑139 (MQ=255)
                                                                                                     ttATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTc                             <  1:609201/139‑1 (MQ=255)
                                                                                                        tCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCAt                           <  1:150718/138‑1 (MQ=255)
                                                                                                            aTCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCGATCGTGGTTCGGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGaa                        <  1:830803/137‑1 (MQ=255)
                                                                                                            aTCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGtt                                                                         <  1:609174/88‑1 (MQ=255)
                                                                                                            aTCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGtt                                                                         >  2:609174/1‑88 (MQ=255)
                                                                                                            aTCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGaa                        <  1:406197/137‑1 (MQ=255)
                                                                                                            aTCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGaa                        <  1:582046/137‑1 (MQ=255)
                                                                                                            aTCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGaa                        <  1:792660/137‑1 (MQ=255)
                                                                                                               aaTATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGta                   <  1:207890/139‑1 (MQ=255)
                                                                                                                      atatCAATCCTTGGCTGCTCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCt            >  2:490622/1‑139 (MQ=255)
                                                                                                                      atatCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTgaga                                                                                                           <  2:266367/44‑1 (MQ=255)
                                                                                                                      atatCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTgaga                                                                                                           >  1:266367/1‑44 (MQ=255)
                                                                                                                      atatCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGaa                                                                                           <  1:417304/60‑1 (MQ=255)
                                                                                                                      atatCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGaa                                                                                           >  2:417304/1‑60 (MQ=255)
                                                                                                                      atatCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTg                                                    >  1:622013/1‑99 (MQ=255)
                                                                                                                      atatCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTg                                                    <  2:622013/99‑1 (MQ=255)
                                                                                                                      atatCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGta                   >  1:775368/1‑132 (MQ=255)
                                                                                                                      atatCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGta                   <  2:775368/132‑1 (MQ=255)
                                                                                                                        atCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAggg                                                                                                <  2:418886/53‑1 (MQ=255)
                                                                                                                        atCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAggg                                                                                                >  1:418886/1‑53 (MQ=255)
                                                                                                                         tCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGaa                        <  2:175989/124‑1 (MQ=255)
                                                                                                                         tCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGaa                        >  1:175989/1‑124 (MQ=255)
                                                                                                                           aaTCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAAcc                                               <  1:401006/99‑1 (MQ=255)
                                                                                                                           aaTCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAAcc                                               >  2:401006/1‑99 (MQ=255)
                                                                                                                           aaTCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCaa         >  1:247439/1‑137 (MQ=255)
                                                                                                                                ttGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATg                          <  2:494431/115‑1 (MQ=255)
                                                                                                                                ttGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATg                          >  1:494431/1‑115 (MQ=255)
                                                                                                                                ttGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGta                   <  1:628428/122‑1 (MQ=255)
                                                                                                                                ttGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGta                   >  2:628428/1‑122 (MQ=255)
                                                                                                                                ttGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATc  >  2:639748/1‑139 (MQ=255)
                                                                                                                                 tGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATc  >  1:691289/1‑138 (MQ=255)
                                                                                                                                        |                                                                                                                                  
GGCTTCTATACCAATTATTGGGTCATTACCCAACTCAATATTAGCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGTTATTTCATGAACTGTATAGTCCTCAACAGGATC  >  NZ_CP009273/4545460‑4545726

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