Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_pae_pgi 4,015,435 Δ1 bp coding (735/816 nt) ysgA ← putative carboxymethylenebutenolidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_pgi4,015,4350G.94.6% 147.1 / 0.6 37coding (735/816 nt)ysgAputative carboxymethylenebutenolidase
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base . (22/13);  total (24/13)
Fisher's exact test for biased strand distribution p-value = 5.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTACACGATAATCTCTGCTTTCGCATTAGCAGCCCGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAA  >  NC_000913_3_pae_pgi/4015294‑4015575
                                                                                                                                             |                                                                                                                                            
cAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATA‑CTCGGGCg                                                                                                                                      >  1:100512/1‑149 (MQ=255)
cAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGAGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGTTA‑CTCGGGCg                                                                                                                                      >  1:258768/1‑149 (MQ=255)
     cTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCGATCTTCTGCAGATGCGGCATGATA‑CTCGGGCGATAAt                                                                                                                                 >  1:252821/1‑149 (MQ=255)
                ggTAAAATCCAAACCGGGCGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCACACCATACTGCTTAAACAATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGGGGCATGCT‑GCTCGGGGGATAATCAGCGTTGAat                                                                                                                      >  1:334922/1‑149 (MQ=255)
                 gTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGAATACGCTACCAGCCATCTTCTGCAGATGCGGCATGATA‑CTCGGGCAATAATCAGCGTTGAatg                                                                                                                     >  1:402280/1‑149 (MQ=255)
                     aaTCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATA‑CTCGGCCGATACTCACCGTTGAatgcctg                                                                                                                 >  2:156793/1‑149 (MQ=255)
                       tCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACCCGACCAGCTACTTTCTTCGGATGCGGCCTGATA‑CTCGGGCGATAATCCGCGTTGAATGCATGcc                                                                                                               >  1:183348/1‑149 (MQ=255)
                        ccAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATA‑CTCGGGCGATAATCAGCGTTGAATGCATGccc                                                                                                              <  2:100512/149‑1 (MQ=255)
                         cAAACCGGGTGGTAATACCACCCGGTCTTGTCTCATTACAGCCACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATA‑CTCGGGCGATAATCAGCGTTGAATGCATGCCCg                                                                                                             <  2:101154/149‑1 (MQ=255)
                         cAAACCGGGCGGTAACACCACCCGGTCCTCGCTCATTACAGCTACTTCTTCCCACCATACTGCTTAAACCATTCCAGTATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATA‑CTCGGGCGATAATCAGCGTTGAATGCATGCCCg                                                                                                             <  2:221724/149‑1 (MQ=255)
                          agaCCGGGCGGTATTACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCCGATGCGACATGATA‑CTCGGGCGCTAATCAGCGTTGAATGCATGccccc                                                                                                            >  1:248900/3‑147 (MQ=255)
                          aaaCCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATA‑CTCGGGCGATAATCAGCGTTGAATGCATCCCCgg                                                                                                            >  2:293714/1‑149 (MQ=255)
                              cGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATA‑CTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTc                                                                                                        >  2:282105/1‑149 (MQ=255)
                              cGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATA‑CTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTc                                                                                                        >  2:235859/1‑149 (MQ=255)
                                  tggtAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATA‑CTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCggg                                                                                                    >  1:142563/1‑149 (MQ=255)
                                                    ttttCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATATGCTGCCAGCCATCTTCTGCAGATGCGGCATGATA‑CTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTACACGATAATCTCTGCt                                                                                  >  2:409755/1‑149 (MQ=255)
                                                    ttttCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATA‑CTCGGGCGCTAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTACACGATAATCTCTGCt                                                                                  >  1:254852/1‑149 (MQ=255)
                                                        ctcATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATA‑CTCGGGCGATAATCAGCGt                                                                                                                           >  1:181713/1‑104 (MQ=255)
                                                        ctcATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATA‑CTCGGGCGATAATCAGCGt                                                                                                                           <  2:181713/104‑1 (MQ=255)
                                                             tACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATA‑CTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTACACGATAATCTCTGCTTTCGCATTa                                                                         <  1:282105/149‑1 (MQ=255)
                                                                   gACTTCCTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATA‑CTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTACACGATAATCTCTGCTTTCGCATTAGCAGcc                                                                   <  2:337769/149‑1 (MQ=255)
                                                                          tcCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATA‑CTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTACACGATAATCTCTGCTTTCGCATTAGCAGCCCGCAgcg                                                            >  1:82810/1‑149 (MQ=255)
                                                                             caccaTACTGCTTAAACCATTTCAGCATACGCTGCCAGCCATCTGCTGCAGATGCGGCATGATA‑CTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTACACGATAATCTCTGCTTTCGCATTAGCAGCCCGCAGCGCCt                                                         <  1:152508/149‑1 (MQ=255)
                                                                                    cTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATA‑CTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTACACGATAATCTCTGCTTTCGCATTAGCAGCCCGCAGCGCCTGGCGCAt                                                  >  1:5108/1‑149 (MQ=255)
                                                                                     tGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATA‑CTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTACACGATAATCTATGCTTTGGAATTAGCAGCCCGCAGCGCCTGGCGCATg                                                 >  1:70537/1‑149 (MQ=255)
                                                                                      gCTTAAGCCATTCCAGCATACGCGGCCAGCCATCTTCTGCAGATGCGGCATGATA‑CTCGGGCGATAATCAGCGTCGAATGCATGCCCGGCGTCCGGGTACACGATAATCTCTGCTTTCGCATTAGCAGCCCGCAGCGCCTGGCGCATgg                                                <  2:161311/149‑1 (MQ=255)
                                                                                      gCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATA‑CTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTACACGATAATCTCTGCTTTCGCATTAGCAGCCCGCAGCGCCTGGCGCATgg                                                >  1:242008/1‑149 (MQ=255)
                                                                                          aaaCCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATA‑CTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTACACGATAATCTCTGCTTTCGCATTAGCAGCCCGCAGCGCCTGGCGCCTGGTTTc                                            >  1:239888/1‑149 (MQ=255)
                                                                                            aCCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGAAGCGGCATGATA‑CTCGGCAGATAATCAGCGTTGAAAGCATGCCCGGCGTCCGGGTACACGATAATCTCTGCTCTCGCATTCGCAACCCGCAGCGCCTGGCGCATCGTTTCaa                                          >  2:78434/1‑149 (MQ=255)
                                                                                             ccATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATA‑CTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTACACGATAATCTCTGCTTTCGCATTAGCAGCCCGCAGCGCCTGGCGCATGGTTTCAAc                                         <  2:142563/149‑1 (MQ=255)
                                                                                                   cAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATA‑CTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTACACGATAATCTCTGCTTTCGCATTAGCAGCCCGCAGCGCCTGGCGCATGGTTTCAACGCtctc                                   >  2:210769/1‑149 (MQ=255)
                                                                                                      cATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATA‑CTCGAGCGATAAGGAGCGTTTAATGCATGCCCGGCGTCCGGGTACACGATAATCTGTGCTTTCGCATTAGCAGCCCGCAGCGCCTGGCCACGGGTTTAAAAGCtctcgtc                                >  2:2173/1‑146 (MQ=255)
                                                                                                           gCTGCCAGCCATCTTCTGCAGATGCGGCATGATA‑CTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTACACGATAATCTCTGCTTTCGCATTAGCAGCCCGCAGCGCCTGGCGCATGGTTTCCACGCTCTCCTGCGGaa                           <  1:359330/149‑1 (MQ=255)
                                                                                                           gCTGCCAGCCATCTTCTGCAGATGCGGCATGATA‑CTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTACACGATAATCTCTGCTTTCGCATTAGCAGCCCGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGaa                           <  2:5108/149‑1 (MQ=255)
                                                                                                                     aTCTTCTGCGGAGGGGGCATGATA‑CTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTACACGATAATCTCTGCTTTCGCATTAGCAGCCCGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATc                 <  2:251344/149‑1 (MQ=255)
                                                                                                                                    ggCATGATA‑CTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTACACGATAATCTCTGCTTTCGCATTAGCAGCCCGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATaa  <  1:249225/149‑1 (MQ=255)
                                                                                                                                     gCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTACACGATAATCTCTGCTTTCGCATTAGCAGCCCGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATaa  >  1:153436/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                            
CAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTACACGATAATCTCTGCTTTCGCATTAGCAGCCCGCAGCGCCTGGCGCATGGTTTCAACGCTCTCCTGCGGAATGCTGTTATCCTGACCACCATATAA  >  NC_000913_3_pae_pgi/4015294‑4015575

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