Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 874,183 G→A 12.5% C402Y (TGT→TAT)  yliE → putative membrane‑anchored cyclic‑di‑GMP phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913874,1830GA12.5% 140.4 / 10.9 48C402Y (TGT→TAT) yliEputative membrane‑anchored cyclic‑di‑GMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base G (20/22);  new base A (3/3);  total (23/25)
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.09e-01

TGAAAGATGCCGAAATTGAACGGTCAAATAAATTGCTCTCACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGCGATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAAACTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTGTGTCAGCGAATTGCCGAGAAATATGCGCAATGCTGCGATATCGTGACATTCAATGCCGATCTCTATTTACTTCT  >  NC_000913/874045‑874321
                                                                                                                                          |                                                                                                                                          
tGAAAGATGCCGAAATTGAACGGTCAAATAAATTGCTCTCACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGCGATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAAACtgtgt                                                                                                                                           <  1:242407/140‑1 (MQ=255)
     gATGCCGAAATTGAACGGTCAAATAAATTGCTCTCACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGCGATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAAACGGTGTGGTCt                                                                                                                                      >  2:1006477/1‑140 (MQ=255)
             aaTTGAACGGTCAAATAAATTGCTCTCACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGCGATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAAACTGTGTGGTCTGGAAGCCg                                                                                                                              <  2:442818/140‑1 (MQ=255)
                     ggTCAAATAAATTGCTCTCACTGATCTCTTAAGATCATGAAAGTGGTTTTATTAAAAATTTGTCGATTTTTGAGTCTTACAATTATCTGTATCTGGATGTGGGGATCATCAAACTGTGTGGTCTGGAAGCCGTGGAAGCg                                                                                                                      <  2:957218/140‑1 (MQ=255)
                          aataaatTGCTCTCACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGCGATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAAACTGTGTGGTCTGGAAGCCGTGGAAGCGGTGtt                                                                                                                 <  1:48715/140‑1 (MQ=255)
                            taaatTGCTCTCACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGCGATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAAACTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTg                                                                                                               <  1:578294/140‑1 (MQ=255)
                            taaatTGCTCTCACTGATCTCTTACGATCAGGAAAGTGGATTTATTAAAAATATGGCGATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAAACTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTg                                                                                                               <  1:1206809/140‑1 (MQ=255)
                                 tGCTCTCACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGCGATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAAACTGTGTGGTCTGGAAGCCGTGGAAGCGGGGTTTGGTGtt                                                                                                          >  1:813257/1‑140 (MQ=255)
                                    tctcACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGCGATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAAACTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTtgat                                                                                                       <  1:700486/140‑1 (MQ=255)
                                       cACTGATCTCTTTCGATCAGGATAGTGGTTTTTTTAAAAATATGGCGATTTTTGAGTCTTACACTAATCAGTATCTGGCTGTGGGGATCATCAAACTGTGTGGTCTGGAAGGCGTGGAAGCGGTTTTTGGTGTTGTTGAg                                                                                                    <  2:716738/140‑1 (MQ=255)
                                       cACTGATCTCTTACTTTCAGGAAAGTGGTTTTATTAAAAATATGGCTATTATTTAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAATCTGTGTTGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAg                                                                                                    <  2:716271/140‑1 (MQ=255)
                                       cACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGCGATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAAACTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAg                                                                                                    <  1:188794/140‑1 (MQ=255)
                                       cACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGCGATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAAACTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAg                                                                                                    <  1:468034/140‑1 (MQ=255)
                                                  tACGATCAGGAAAGTGGTTTTATTAAAAATATGGCGATTATTTAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAAACTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAAt                                                                                         >  2:1160138/1‑140 (MQ=255)
                                                  tACGATCAGGAAAGTGGTTTTATTAAAAATATGGCGATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAAACTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAAt                                                                                         <  2:527429/140‑1 (MQ=255)
                                                  tACGATCAGGAAAGTGGTTTTATTAAAAATATGGCGATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAAACTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAAt                                                                                         >  1:561207/1‑140 (MQ=255)
                                                  tACGATCAGGAAAGTGGTTTTATTAAAAATATGGCGATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAAACTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAAt                                                                                         >  1:561192/1‑140 (MQ=255)
                                                  tACGATCAGGAAAGTGGTTTTATTAAAAATATGGCGATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAAACTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAAt                                                                                         >  1:484558/1‑140 (MQ=255)
                                                      aTCAGGAAAGTGGTTTTATTAAAAATATGGCGATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAAACTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTc                                                                                     >  1:1060996/1‑140 (MQ=255)
                                                      aTCAGGAAAGTGGTTTTATTAAAAATATGGCGATTATTGAGTCTAACAATAATCAGTATCTGGCTGGGGGGATCATCAAACTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTc                                                                                     >  2:1118949/1‑140 (MQ=255)
                                                               gTGGTTTTATTAAAAATATGGCGATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAAACTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATtg                                                                            >  1:31455/1‑140 (MQ=255)
                                                                                  ggCGATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAAACTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTGTGTCAGCGAATTGCCgaga                                                         >  2:756732/1‑140 (MQ=255)
                                                                                  ggCGATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAAACTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTGTGTCAGCGAATTGCCgaga                                                         >  2:756167/1‑140 (MQ=255)
                                                                                   gCGATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAAACTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTGTGTCAGCGAATTGCCGAGaa                                                        <  2:31455/140‑1 (MQ=255)
                                                                                     gATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAAACTGTATGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCa                                                                                                >  1:315674/1‑98 (MQ=255)
                                                                                     gATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAAACTGTATGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCa                                                                                                <  2:315674/98‑1 (MQ=255)
                                                                                        tattGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAAACTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTGTGTCAGCGAATTGCCGAGAAATATg                                                   <  2:465136/140‑1 (MQ=255)
                                                                                          ttGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAAACTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTGTGTCAGCGAATTGCCGAGAAATATgcg                                                 <  1:232277/140‑1 (MQ=255)
                                                                                                            aGTATCTGGCTGTGGGGATCATCAAACTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTGTGTCAGCGAATTGCCGAGAAATATGCGCAATGCTGCGATATCGTg                               <  2:940873/140‑1 (MQ=255)
                                                                                                            aGTATCTGGCTGTGGGGATCATCAAACTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTGTGTCAGCGAATTGCCGAGAAATATGCGCAATGCTGCGATATCGTg                               <  2:1150628/140‑1 (MQ=255)
                                                                                                                 cTGGCTGTGGGGATCATCAAACTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTGTGTCAGCGAATTGCCGAGAAATATGCGCAATGCTGCGATATCGTGACAtt                          >  1:855899/1‑140 (MQ=255)
                                                                                                                   ggCTGTGGGGATCATCAAACTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTGTGTCAGCGAATTGCCGAGAAATATGCGCAATGCTGCGATATCGTGACATTCa                        >  2:687358/1‑140 (MQ=255)
                                                                                                                        tgGGGATCATCAAACTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTGTGTCAGCGAATTGCCGAGAAATATGCGCAATGCTGCGATATCGTGACATTCAATGcc                   >  1:895878/1‑140 (MQ=255)
                                                                                                                        tgGGGATCATCAAACTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTGTGTCAGCGAATTGCCGAGAAATATGCGCAATGCTGCGATATCGTGACATTCAATGcc                   >  1:1174344/1‑140 (MQ=255)
                                                                                                                        tgGGGATCATCAAACTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTGTGTCAGCGAATTGCCGAGAAATATGCGCAATGCTGCGATATCGTGACATTCAATGcc                   >  1:895551/1‑140 (MQ=255)
                                                                                                                        tgGGGATCATCAAACTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTGTGTCAGCGAATTGCCGAGAAATATGCGCAATGCTGCGATATCGTGACATTCAATGcc                   >  2:853323/1‑140 (MQ=255)
                                                                                                                        tgGGGATCATCAAACTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTGTGTCAGCGAATTGCCGAGAAATATGCGCAATGCTGCGATATCGTGACATTCAATGcc                   >  2:864633/1‑140 (MQ=255)
                                                                                                                        tgGGGATCATCAAACTGTATGGTCTGGAAGCCGTGGAAGc                                                                                                                       <  2:180998/40‑1 (MQ=255)
                                                                                                                        tgGGGATCATCAAACTGTATGGTCTGGAAGCCGTGGAAGc                                                                                                                       <  2:181108/40‑1 (MQ=255)
                                                                                                                        tgGGGATCATCAAACTGTATGGTCTGGAAGCCGTGGAAGc                                                                                                                       >  1:180998/1‑40 (MQ=255)
                                                                                                                        tgGGGATCATCAAACTGTATGGTCTGGAAGCCGTGGAAGc                                                                                                                       >  1:181108/1‑40 (MQ=255)
                                                                                                                           ggATCATCAAACTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTGTGTCAGCGAATTGCCGAGAAATATGCGCAATGCTGCGATATCGTGACATTCAATGCCGAt                <  1:199436/140‑1 (MQ=255)
                                                                                                                             atcatcAAACTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTGTGTCAGCGAATTGCCGAGAAATATGCGCAATGCTGCGATATCGTGACATTCAATGCCGAtct              <  2:400997/140‑1 (MQ=255)
                                                                                                                               catcaAACTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTGTGTCAGCGAATTGCCGAGAAATATGCGCAATGCTGCGATATCGTGACATTCAATGCCGATctct            >  1:159725/1‑140 (MQ=255)
                                                                                                                                  caAACTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTGTGTCAGCGAATTGCCGAGAAATATGCGCAATGCTGCGATATCGTGACATTCAATGCCGATCTCTAtt         >  1:1133332/1‑140 (MQ=255)
                                                                                                                                   acaCTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAAACGTCAGGAAATTGTGTCAGCGAATTGCCGAGAAATATGCGCAATGCTGCGATATCGTGACAATCAATGCCGATCTCTAttt        <  2:352014/138‑1 (MQ=255)
                                                                                                                                   aaaCTGTGTGGTCTGGAAGCCGTGGAATCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTGTGTCAGCGAATTGCCGAGAAATATGCGCAATGCTGCGATATCGTGACATTCAATGCCGATCTCTAttt        <  2:352005/140‑1 (MQ=255)
                                                                                                                                   aaaCTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTGTGTCAGCGAATTGCCGAGAAATATGCGCAATGCTGCGATATCGTGACATTCAATGCCGATCTCTAttt        <  2:1034045/140‑1 (MQ=255)
                                                                                                                                   aaaCTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTGTGTCAGCGAATTGCCGAGAAATATGCGCAATGCTGCGATATCGTGACATTCAATGCCGATCTCTAttt        <  2:642462/140‑1 (MQ=255)
                                                                                                                                         tgtgGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTGTGTCAGCGAATTGCCGAGAAATATGCGCAATGCTGCGATATCGTGACATTCAATGCCGATCTCTATTTACTTCt  >  1:1016205/1‑140 (MQ=255)
                                                                                                                                         tgtgGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTGTGTCAGCGAATTGCCGAGAAATATGCGCAATGCTGCGATATCGTGACATTCAATGCCGATCTCTATTTACTTCt  >  1:346349/1‑140 (MQ=255)
                                                                                                                                         tgtgGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTGTGTCAGCGAATTGCCGAGAAATATGCGCAATGCTGCGATATCGTGACATTCAATGCCGATCTCTATTTACTTCt  >  1:598738/1‑140 (MQ=255)
                                                                                                                                         tgtgGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTGTGTCAGCGAATTGCCGAGAAATATGCGCAATGCTGCGATATCGTGACATTCAATGCCGATCTCTATTTACTTCt  >  1:872165/1‑140 (MQ=255)
                                                                                                                                         tgtgGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTGTGTCAGCGAATTGCCGAGAAATATGCGCAATGCTGCGATATCGTGACATTCAATGCCGATCTCTATTTACTTCt  >  1:925879/1‑140 (MQ=255)
                                                                                                                                         tgtgGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTGTGTCAGCGAATTGCCGAGAAATATGCGCAATGCTGCGATATCGTGACATTCAATGCCGATCTCTATTTACTTCt  >  1:926200/1‑140 (MQ=255)
                                                                                                                                         tgtgGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTGTGTCAGCGAATTGCCGAGAAATATGCGCAATGCTGCGATATCGTGACATTCAATGCCGATCTCTATTTACTTCt  >  1:1206216/1‑140 (MQ=255)
                                                                                                                                         tgtgGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTGTGTCAGCGAATTGCCGAGAAATATGCGCAATGCTGCGATATCGTGACATTCAATGCCGATCTCTATTTACTTCt  >  2:176423/1‑140 (MQ=255)
                                                                                                                                         tgtgGTCTGGAAGCCGTGGAAGCGGAGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTGTGTCAGCGAATAGCCGAGAAATATGCGCAATGCTGCGATATCGTGACATTCAATGCCGATCTCTATTTACTTCt  >  1:874834/1‑140 (MQ=255)
                                                                                                                                          |                                                                                                                                          
TGAAAGATGCCGAAATTGAACGGTCAAATAAATTGCTCTCACTGATCTCTTACGATCAGGAAAGTGGTTTTATTAAAAATATGGCGATTATTGAGTCTAACAATAATCAGTATCTGGCTGTGGGGATCATCAAACTGTGTGGTCTGGAAGCCGTGGAAGCGGTGTTTGGTGTTGATGAGCGCAATAAAATCGTCAGGAAATTGTGTCAGCGAATTGCCGAGAAATATGCGCAATGCTGCGATATCGTGACATTCAATGCCGATCTCTATTTACTTCT  >  NC_000913/874045‑874321

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