Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,490,055 T→G E3D (GAA→GAC yjgR ← DUF853 family protein with NTPase fold

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,490,0550TG98.3% 203.9 / ‑3.8 58E3D (GAA→GACyjgRDUF853 family protein with NTPase fold
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base G (28/29);  total (29/29)
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.85e-01

CCCGTCCCCGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGTTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAACCTAAAGCGTG  >  NC_000913/4489960‑4490150
                                                                                               |                                                                                               
cccGTCCCCGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCAc                                                                                             >  2:892148/1‑100 (MQ=255)
      cccGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATga                                                                                       <  2:528295/100‑1 (MQ=255)
        cGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATgaga                                                                                     <  1:364565/100‑1 (MQ=255)
        cGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATgaga                                                                                     <  2:438401/100‑1 (MQ=255)
          ttGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGata                                                                                   <  1:1366314/100‑1 (MQ=255)
             cgcCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATAtcc                                                                                >  1:562550/1‑100 (MQ=255)
              gcCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATAtcct                                                                               <  1:1404379/100‑1 (MQ=255)
                cGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATAtcctcc                                                                             <  1:444882/100‑1 (MQ=255)
                cGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATAtcctcc                                                                             <  2:274304/100‑1 (MQ=255)
                 ggTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCa                                                                            >  1:1632458/1‑100 (MQ=255)
                  gTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAg                                                                           <  2:1779233/100‑1 (MQ=255)
                  gTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAg                                                                           <  2:1458438/100‑1 (MQ=255)
                  gTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAg                                                                           <  2:968034/100‑1 (MQ=255)
                     aTCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGttt                                                                        <  2:193808/100‑1 (MQ=255)
                         gCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTcctg                                                                    >  2:1348584/1‑100 (MQ=255)
                              tGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTgg                                                               <  2:51459/100‑1 (MQ=255)
                                 cGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGt                                                            >  1:1785269/1‑100 (MQ=255)
                                 cGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGt                                                            <  2:336978/100‑1 (MQ=255)
                                  ggTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGta                                                           <  1:1907987/100‑1 (MQ=255)
                                   gTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGtat                                                          >  1:314985/1‑100 (MQ=255)
                                    ttAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGtata                                                         <  2:404635/100‑1 (MQ=255)
                                    ttAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCAGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGtata                                                         <  1:650084/100‑1 (MQ=255)
                                     tAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAg                                                        <  1:452418/100‑1 (MQ=255)
                                       gCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGc                                                      >  2:266076/1‑100 (MQ=255)
                                       gCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGc                                                      >  2:586346/1‑100 (MQ=255)
                                       gCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGc                                                      >  1:1012514/1‑100 (MQ=255)
                                            ncccGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGt                                                 >  2:1523809/2‑100 (MQ=255)
                                              ccGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTcc                                               >  2:1554833/1‑100 (MQ=255)
                                                 gCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCgg                                             <  1:101841/99‑1 (MQ=255)
                                                    gTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCaaa                                         <  1:610266/100‑1 (MQ=255)
                                                    gTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCaaa                                         >  1:699949/1‑100 (MQ=255)
                                                     tAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAag                                        <  2:910625/100‑1 (MQ=255)
                                                      aaaaaCAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAaga                                       >  1:1586874/1‑100 (MQ=255)
                                                      aaaaaCAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAaga                                       >  1:332715/1‑100 (MQ=255)
                                                       aaaaCAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAagag                                      >  2:1271943/1‑100 (MQ=255)
                                                        aaaCAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAgaga                                     >  1:1290525/1‑100 (MQ=255)
                                                         aaCAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGaa                                    >  1:1150341/1‑100 (MQ=255)
                                                          aCAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAAt                                   <  2:1035183/100‑1 (MQ=255)
                                                          aCAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAAt                                   >  2:1413512/1‑100 (MQ=255)
                                                                   tgtCCGGCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAAt                          >  2:49867/1‑100 (MQ=255)
                                                                        ggCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAaga                     <  2:7344/100‑1 (MQ=255)
                                                                        ggCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAaga                     >  2:1231627/1‑100 (MQ=255)
                                                                        ggCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAaga                     <  2:381015/100‑1 (MQ=255)
                                                                        ggCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCCGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAaga                     >  1:650147/1‑100 (MQ=255)
                                                                         gCGTGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAagag                    <  1:977345/100‑1 (MQ=255)
                                                                            tGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGaaa                 <  1:892147/100‑1 (MQ=255)
                                                                            tGCGGGCAATTAACAGGGGGTCACTCATGAGATATCCGCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGaaa                 >  1:1045432/1‑100 (MQ=255)
                                                                             gCGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGaaaa                <  2:1685697/100‑1 (MQ=255)
                                                                              cGGGCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAAt               >  1:44141/1‑100 (MQ=255)
                                                                               gggCAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATa              <  2:699950/100‑1 (MQ=255)
                                                                                  cAATTAACAGGGGTTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAAcc           >  2:1673861/1‑100 (MQ=255)
                                                                                  cAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAAcc           >  1:1886372/1‑100 (MQ=255)
                                                                                  cAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAAcc           <  2:325024/100‑1 (MQ=255)
                                                                                  cAATTAACAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAAcc           <  2:1045433/100‑1 (MQ=255)
                                                                                       aaCAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAACCTAAAg      >  1:105936/1‑100 (MQ=255)
                                                                                        aCAGGGGGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAACCTAAAGc     >  1:1473539/1‑100 (MQ=255)
                                                                                           ggggGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAACCTAAAGCgtg  >  2:946168/1‑100 (MQ=255)
                                                                                           ggggGTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAACCTAAAGCgtg  >  2:1304717/1‑100 (MQ=255)
                                                                                               |                                                                                               
CCCGTCCCCGTTGCGCCGGTAATCAGCCCATGACGGTTAGCCATTCCCGGCAGTAAAAACAGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGTTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAACCTAAAGCGTG  >  NC_000913/4489960‑4490150

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