Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,684,183 A→G V479A (GTG→GCG)  yhjK ← cyclic‑di‑GMP phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,684,1830AG97.4% 132.0 / ‑3.8 39V479A (GTG→GCG) yhjKcyclic‑di‑GMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (22/16);  total (22/17)
Fisher's exact test for biased strand distribution p-value = 4.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.68e-01

TAACAGTTCCAGCATATCCGCCACCATATTCGGGTGCATCAGTTGCAGCGCAGAGAGGTTTACCGACAAGGGCAGCATAATGCCGCGCTCTTGCCAGGCTGCAAGCAATCGACAGGACTCTTCCAGCACCCAGTGACCGACGGTAACCATCAGCCCACAGCACTCAATGCGATCGATTAAGCCATCCGGCAGGTCCCAACTGCCATCCGGTTGCTGGATACGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGG  >  NC_000913/3684056‑3684328
                                                                                                                               |                                                                                                                                                 
tAACAGTTCCAGCATATCCGCCACCATATTCGGGTGCATCAGTTGCAGCGCAGAGAGGTTTACCGACAAGGGCAGCATAATGCCGCGCTCTTGCCAGGCTGCAAGCAATCGACAGGACTCTTCCAGCGCCCAGTGACCGACGGTAACCAt                                                                                                                             <  1:405953/150‑1 (MQ=255)
          aGCATATCCGCCACCATATTCGGGTGCATCGGTTGCAGCGCAGAGAGGTTTACCGACAAGGGCAGCATAATGCCGCGCTCTTGCCAGGCTGCAAGCAATCGACAGGACTCTTCCAGCGCCCAGTGACCGACGGTAACCATCAGCCCACAg                                                                                                                   >  1:588210/1‑150 (MQ=255)
          aGCATATCCGCCACCATATTCGGGTGCATCAGTTGCAGCGCAGAGAGGTTTACCGACAAGGGCAGCATAATGCCGCGCTCTTGCCAGGCTGCAAGCAATCGACAGGACTCTTCCAGCGCCCAGTGACCGACGGTAACCATCAGCCCACAGc                                                                                                                  <  1:269981/151‑1 (MQ=255)
                  cGCCACCATATTCGGGTGCATCCGTTGCAGCGCAGAGAGGTTTACCGACAAGGGCAGCATAATGCCGCGCTCTTGCCAGGCTGCAAGCAATCGACAGGACTCTTCCAGCGCCCAGTGACCGACGGTAACCATCAGCCCACAGCACTCAAt                                                                                                           >  1:573467/1‑150 (MQ=255)
                  cGCCACCATATTCGGGTGCATCAGTTGCAGCGCAGAGAGGTTTACCGACAAGGGCAGCATAATGCCGCGCTCTTGCCAGGCTGCAAGCAATCGACAGGACTCTTCCAGCGCCCAGTGACCGACGGTAACCATCAGCCCACAGCACTCAATg                                                                                                          >  1:588486/1‑151 (MQ=255)
                         atatTCGGGTGCATCAGTTGCAGCGCAGAGAGGTTTACCGACAAGGGCAGCATAATGCCGCGCTCTTGCCAGGCTGCAAGCAATCGACAGGACTCTTCCAGCGCCCAGTGACCGACGGTAACCATCAGCCCACAGCACTCAATGcgatcg                                                                                                    >  1:652817/1‑150 (MQ=255)
                           atTCGGGTGCATCAGTTGCAGCGCAGAGAGGTTTACCGACAAGGGCAGCATAATGCCGCGCTCTTGCCAGGCTGCAAGCAATCGACAGGACTCTTCCAGCGCCCAGTGACCGACGGTAACCATCAGCCCACAGCACTCAATGCGATCGAtt                                                                                                 >  2:162715/1‑151 (MQ=255)
                           atTCGGGTGCATCAGTTGCAGCGCAGAGAGGTTTACCGACAAGGGCAGCATAATGCCGCGCTCTTGCCAGGCTGCAAGCAATCGACAGGACTCTTCCAGCGCCCAGTGACCGACGGTAACCATCAGCCCACAGCACTCAATGCGATCGAtt                                                                                                 <  2:389966/151‑1 (MQ=255)
                             tCGGGTGCATCAGTTGCAGCGCAGAGAGGTTTACCGACAAGGGCAGCATAATGCCGCGCTCTTGCCAGGCTGCAAGCAATCGACAGGACTCTTCCAGCGCCCAGTGACCGACGGTAACCATCAGCCCACAGCACTCAATGCGATCGATTa                                                                                                <  1:395892/150‑1 (MQ=255)
                                   gCATCAGTTGCAGCGCAGAGAGGTTTACCGACAAGGGCAGCATAATGCCGCGCTCTTGCCAGGCTGCAAGCAATCGACAGGACTCTTCCAGCGCCCAGTGACCGACGGTAACCATCAGCCCACAGCACTCAATGCGATCGATTAAGCCATc                                                                                         >  2:198323/1‑151 (MQ=255)
                                     aTCAGTTGCAGCGCAGAAAGGTTTACCGACAAGGGCAGCATAATGCCGCGCTCTTGCCAGGCTGCAAGCAATCGACAGGACTCTTCCAGCGCCCAGTGACCGACGGTAACCATCAGCCCACAGCACTCAATGCGATCGATTAAGCCATCCg                                                                                       >  2:495565/1‑151 (MQ=255)
                                      tCAGTTGCAGCGCAGAGAGGTTTACCGACAAGGGCAGCATAATGCCGCGCTCTTGCCAGGCTGCAAGCAATCGACAGGACTCTTCCAGCGCCCAGTGACCGACGGTAACCATCAGCCCACAGAACTCAATGCGATCGATTAAGCCATCCgg                                                                                      >  2:467666/1‑151 (MQ=255)
                                        aGTTGCAGCGCAGAGAGGTTTACCGACAAGGGCAGCATAATGCCGCGCTCTTGCCAGGCTGCAAGCAATCGACAGGACTCTTCCAGCGCCCAGTGACCGACGGTAACCATCAGCCCACAGCACTCAATGCGATCGATTAAGCCATCCGGc                                                                                     >  2:327986/1‑150 (MQ=255)
                                               gcgcAGAGAGGTTTACCGACAAGGGCAGCATAATGCCGCGCTCTTGCCAGGCTGCAAGCAATCGACAGGACTCTTCCAGCGCCCAGTGACCGACGGTAACCATCAGCCCACAGCACTCAATGCGATCGATTAAGCCATCCGGCAGGTCCCa                                                                             <  2:22957/151‑1 (MQ=255)
                                                cgcAGAGAGGTTTACCGACAAGGGCAGCATAATGCCGCGCTCTTGCCAGGCTGCAAGCAATCGACAGGACTCTTCCAGCGCCCAGTGACCGACGGTAACCATCAGCCCACAGCACTCAATGCGATCGATTAAGCCATCCGGCAGGTCCCa                                                                             >  2:144713/1‑150 (MQ=255)
                                                   agagagGTTTACCGACAAGGGCAGCATAATGCCGCGCTCTTGCCAGGCTGCAAGCAATCGACAGGACTCTTCCAGCGCCCAGTGACCGACGGTAACCATCAGCCCACAGCACTCAATGCGATCGATTAAGCCATCCGGCAGGTCCCAACTg                                                                         <  1:607907/151‑1 (MQ=255)
                                                   agagagGTTTACCGACAAGGGCAGCATAATGCCGCGCTCTTGCCAGGCTGCAAGCAATCGACAGGACTCTTCCAGCGCCCAGTGACCGACGGTAACCATCAGCCCACAGCACTCAATGCGATCGATTAAGCCATCCGGCAGGTCCCAACTg                                                                         >  2:208577/1‑151 (MQ=255)
                                                       agGTTTACCGACAAGGGCAGCATAATGCCGCGCTCTTGCCAGGCTGCAAGCAATCGACAGGACTCTTCCAGCGCCCAGTGACCGACGGTAACCATCAGCCCACAGCACTCAATGCGATCGATTAAGCCATCCGGCAGGTCCCAACTGCCAt                                                                     >  1:277669/1‑151 (MQ=255)
                                                        ggTTTACCGACAAGGGCAGCATAATGCCGCGCTCTTGCCAGGCTGCAAGCAATCGACAGGACTCTTCCAGCGCCCAGTGACCGACGGTAACCATCAGCCCACAGCACTCAATGCGATCGATTAAGCCATCCGGCAGGTCCCAACTGCCATc                                                                    <  1:261773/151‑1 (MQ=255)
                                                               cGACAAGGGCAGCATAATGCCGCGCTCTTGCCAGGCTGCAAGCAATCGACAGGACTCTTCCAGCGCCCAGTGACCGACGGTAACCATCAGCCCACAGCACTCAATGCGATCGATTAAGCCATCCGGCAGGTCCCAACTGCCATCCGGTTGc                                                             >  2:462459/1‑151 (MQ=255)
                                                                gACAAGGGCAGCATAATGCCGCGCTCTTGCCAGGCTGCAAGCAATCGACAGGACTCTTCCAGCGCCCAGTGACCGACGGTAACCATCAGCCCACAGCACTCAATGCGATCGATTAAGCCATCCGGCAGGTCCCAACTGCCATCCGGTTGCt                                                            <  1:291486/151‑1 (MQ=255)
                                                                gACAAGGGCAGCATAATGCCGCGCTCTTGCCAGGCTGCAAGCAATCGACAGGACTCTTCCAGCGCCCAGTGACCGACGGTAACCATCAGCCCACAGCACTCAATGCGATCGATTAAGCCATCCGGCAGGTCCCAACTGCCATCCGGTTGCt                                                            <  2:377601/151‑1 (MQ=255)
                                                                          gcaTAATGCCGCGCTCTTGCCAGGCTGCAAGCAATCGACAGGACTCTTCCAGCGCCCAGTGACCGACGGTAACCATCAGCCCACAGCACTCAATGCGATCGATTAAGCCATCCGGCAGGTCCCAACTGCCATCCGGTTGCTGGATACGc                                                    >  1:677957/1‑149 (MQ=255)
                                                                            aTACTGCCGCGCTCTTGCCAGGCTGCAAGCAATCGACAGTACTCTTCCAGCGCCCAGTGACCGACGGTAACCATCAGACCACAGCACACAATGCGATCGATTAAGCCATCCGGCAGGTCCCAACTGCCATCCGGTTGCTGGATACGCAGt                                                 <  2:136164/150‑1 (MQ=255)
                                                                                 gCCGCGCTCTTGCCAGGCTGCAAGCAATCGACAGGACTCTTCCAGCGCCCAGTGACCGACGGTAACCATCAGCCCACAGCACTCAATGCGATCGATTAAGCCATCCGGCAGGTCCCAACTGCCATCCGGTTGCTGGATACGCAGTAacac                                            <  1:323021/150‑1 (MQ=255)
                                                                                      gcTCTTGCCAGGCTGCAAGCAATCGACAGGACTCTTCCAGCGCCCAGTGACCGACGGTAACCATCAGCCCACAGCACTCAATGCGATCGATTAAGCCATCCGGCAGGTCCCAACTGCCATCCGGTTGCTGGATACGCAGTAACACTTCCGc                                      >  2:483795/1‑151 (MQ=255)
                                                                                           tGCCAGGCTGCAAGCAATCGACAGGACTCTTCCAGCGCCCAGTGACCGACGGTAACCATCAGCCCACAGCACTCAATGCGATCGATTAAGCCATCCGGCAGGTCCCAACTGCCATCCGGTTGCTGGATACGCAGTAACACTTCCGCACTg                                  <  1:372198/150‑1 (MQ=255)
                                                                                                   tGCAAGCAATCGACAGGACTCTTCCAGCGCCCAGTGACCGACGGTAACCATCAGCCCACAGCACTCAATGCGATCGATTAAGCCATCCGGCAGGTCCCAACTGCCATCCGGTTGCTGGATACGCAGTAACACTTCCGCACTGACCAGttt                          >  2:94335/1‑150 (MQ=255)
                                                                                                     caagcaaTCGACAGGACTCTTCCAGCGCCCAGTGACCGACGGTAACCATCAGCCCACAGCACTCAATGCGATCGATTAAGCCATCCGGCAGGTCCCAACTGCCATCCGGTTGCTGGATACGCAGTAACACTTCCGCACTGACCAGTTTAcc                       >  2:131593/1‑151 (MQ=255)
                                                                                                      aagcaaTCGACAGGACTCTTCCAGCGCCCAGTGACCGACGGTAACCATCAGCCCACAGCACTCAATGCGATCGATTAAGCCATCCGGCAGGTCCCAACTGCCATCCGGTTGCTGGATACGCAGTAACACTTCCGCACTGACCAGTTTACCg                      >  2:229063/1‑151 (MQ=255)
                                                                                                      aagcaaTCGACAGGACTCTTCCAGCGCCCAGTGACCGACGGTAACCATCAGCCCACAGCACTCAATGCGATCGATTAAGCCATCCGGCAGGTCCCAACTGCCATCCGGTTGCTGGATACGCAGTAACACTTCCGCACTGACCAGTTTACCg                      >  2:34933/1‑151 (MQ=255)
                                                                                                       agcaaTCGACAGGACTCTTCCAGCGGCCAGTGACCGACGGTAACCATCAGCCCACAGCACTCAATGCGATCGATTAAGCCATCCGGCAGGTCCCAACTGCCATCCGGTTGCTGTAGACGCAGTAACACTTCCGCACTGACCATTTTACCGc                     <  1:401916/151‑1 (MQ=255)
                                                                                                           aTCGACAGGACTCTTCCAGCGCCCAGTGACCGACGGTAACCATCAGCCCACAGCACTCAATGCGATCGATTAAGCCATCCGGCAGGTCCCAACTGCCATCCGGTTGCTGGATACGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGt                 >  1:687547/1‑151 (MQ=255)
                                                                                                             cGACAGGACTCTTCCAGCGCCCAGTGACCGACGGTAACCATCAGCCCACAGCACTCAATGCGATCGATTAAGCCATCCGGCAGGTCCCAACTGCCATCCGGTTGCTGGATACGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCa               <  1:78545/151‑1 (MQ=255)
                                                                                                                       ctTCCAGCGCCCAGTGACCGACGGTAACCATCAGCCCACAGCACTCAATGCGATCGATTAAGCCATCCGGCAGGTCCCAACTGCCATCCGGTTGCTGGATACGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCtg     >  1:44868/1‑151 (MQ=255)
                                                                                                                       ctTCCAGCACCCAGTGACCGACGGTAACCATCAGCCCACAGCACTCAATGCGATCGATTAAGCCATCCGGCAGGTCCCAACTGCCATCCGGTTGCTGGATACGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGAcc       <  1:101804/149‑1 (MQ=255)
                                                                                                                         tCCAGCGCCCAGTGACCGACGGTAACCATCAGCCCACAGCACTCAATGCGATCGATTAAGCCATCCGGCAGGTCCCAACTGCCATCCGGTTGCTGGATACGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCtgt    <  2:394974/150‑1 (MQ=255)
                                                                                                                          ccAGCGCCCAGTGACCGACGGTAACCATCAGCCCACAGCACTCAATGCGATCGATTAAGCCATCCGGCAGGTCCCAACTGCCATCCGGTTGCTGGATACGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCtgtg   <  1:184956/150‑1 (MQ=255)
                                                                                                                           cAGCGCCCAGTGACCGACGGTAACCATCAGCCCACAGCACTCAATGCGATCGATTAAGCCATCCGGCAGGTCCCAACTGCCATCCGGTTGCTGGATACGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTgg  >  1:10824/1‑150 (MQ=255)
                                                                                                                               |                                                                                                                                                 
TAACAGTTCCAGCATATCCGCCACCATATTCGGGTGCATCAGTTGCAGCGCAGAGAGGTTTACCGACAAGGGCAGCATAATGCCGCGCTCTTGCCAGGCTGCAAGCAATCGACAGGACTCTTCCAGCACCCAGTGACCGACGGTAACCATCAGCCCACAGCACTCAATGCGATCGATTAAGCCATCCGGCAGGTCCCAACTGCCATCCGGTTGCTGGATACGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGG  >  NC_000913/3684056‑3684328

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