Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,212,233 A→G intergenic (‑230/+470) iraM ← / ← ycgX RpoS stabilzer during Mg starvation, anti‑RssB factor/DUF1398 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,212,2330AG96.4% 94.5 / ‑6.0 28intergenic (‑230/+470)iraM/ycgXRpoS stabilzer during Mg starvation, anti‑RssB factor/DUF1398 family protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (11/16);  minor base C (0/1);  total (11/17)
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

AAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATACATTCAAATGACTGAACATCCTGTAATTAAAACATAGCCTTTATGCTACTTTGTGCCAATTTGCTAAACATTATGGTTGCCTTTTTATATAACGATAATAATGAATATAAGCATGACA  >  NC_000913/1212099‑1212351
                                                                                                                                      |                                                                                                                      
aaatTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATGCAtt                                                                                                                    <  1:231715/139‑1 (MQ=255)
      attaaAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATg                                                                                                              <  2:997008/139‑1 (MQ=255)
      attaaAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATg                                                                                                              <  1:414208/139‑1 (MQ=255)
      attaaAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATg                                                                                                              >  2:138457/1‑139 (MQ=255)
                tgGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATc                                                                                                    <  2:863032/139‑1 (MQ=255)
                    ttATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGt                                                                                                <  2:906198/139‑1 (MQ=255)
                      atataAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTaa                                                                                              >  2:340285/1‑139 (MQ=255)
                            aataaaatTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAATTaaaa                                                                                        >  1:410692/1‑139 (MQ=255)
                            aataaaatTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAATTaaaa                                                                                        >  1:871137/1‑139 (MQ=255)
                                 aatTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAATTAAAACATAg                                                                                   >  2:174954/1‑139 (MQ=255)
                                    ttAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTTTAATTAAAACATAGCCt                                                                                <  2:631228/139‑1 (MQ=255)
                                                            gTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAATTAAAACATAGCCTTTATGCTACTTTGTGCCAATTTGc                                                        >  2:981386/1‑139 (MQ=255)
                                                             ttAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAATTAAAACATAGCCTTTATGCTACTTTGTGCCAATTTGCt                                                       <  1:174954/139‑1 (MQ=255)
                                                              tAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAATTAAAACATAGCCTTTATGCTACTTTGTGCCAATTTGCTa                                                      <  2:265343/139‑1 (MQ=255)
                                                                     aaTGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAATTAAAACATAGCCTTTATGCTACTTTGTGCCAATTTGCTAAACATTa                                               >  1:451792/1‑139 (MQ=255)
                                                                         agaTTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAATTAAAACATAGCCTTTATGCTACTTTGTGCCAATTTGCTAAACATTATGGt                                           <  2:712409/139‑1 (MQ=255)
                                                                                 tgAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAATTAAAACAt                                                                                     >  1:863465/1‑89 (MQ=255)
                                                                                 tgAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAATTAAAACAt                                                                                     <  2:863465/89‑1 (MQ=255)
                                                                                                attGACTTCGTTTTAGATTTGTTTAGCTATATTGTTATCCATGTAAGTGACTTATCAGCCTGTAATTAAAACATAGCCTTTAGACTACGTTGTGCCTATTTGCTACACATTATGGTGGCCTTTCTGTATACCGataata                    <  1:142028/139‑1 (MQ=255)
                                                                                                 ttGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAATTAAAACATAGCCTTTATGCTACTTTGTGCCAATTTGCTAAACATTATGGTTGCCTTTTTATATAACGAtaataa                   >  2:729368/1‑139 (MQ=255)
                                                                                                  tGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAATTAAAACATAGCCTTTATGCTACTTTGTGCCAATTTGCTAAACATTATGGTTGCCTTTTTATATAACGATaataat                  <  1:138457/139‑1 (MQ=255)
                                                                                                         gTTTTAGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAATTAAAACATAGCCTTTATGCTACTTTGTGCCAATTTGCTAAACATTATGGTTGCCTTTTTATATAACGATAATAATGAAtata           <  1:729368/139‑1 (MQ=255)
                                                                                                              aGATTTGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAATTAAAACATAGCCTTTATGCTACTTTGTGCCAATTTGCTAAACATTATGGTTGCCTTTTTATATAACGATAATAATGAATATAAGcat      <  1:695927/139‑1 (MQ=255)
                                                                                                                  ttGTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAATTAAAACATAGCCTTTATGCTACTTTGTGCCAATTTGCTAAACATTATGGTTGCCTTTTTATATAACGATAATAATGAATATAAGcatgaca  <  2:451792/139‑1 (MQ=255)
                                                                                                                    gTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAATTAAAACATAGCCTTTATGCTACTTTGTGCCAATTTGCt                                                       >  2:390867/1‑84 (MQ=255)
                                                                                                                    gTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAATTAAAACATAGCCTTTATGCTACTTTGTGCCAATTTGCt                                                       <  1:390867/84‑1 (MQ=255)
                                                                                                                    gTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAATTAAAACATAGCCTTTATGCTACTTTGTGCCAATTTGCt                                                       <  1:326344/84‑1 (MQ=255)
                                                                                                                    gTTTAGCTATAATGTTATGCATTCAAATGACTGAACATCCTGTAATTAAAACATAGCCTTTATGCTACTTTGTGCCAATTTGCt                                                       >  2:326344/1‑84 (MQ=255)
                                                                                                                                      |                                                                                                                      
AAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATACATTCAAATGACTGAACATCCTGTAATTAAAACATAGCCTTTATGCTACTTTGTGCCAATTTGCTAAACATTATGGTTGCCTTTTTATATAACGATAATAATGAATATAAGCATGACA  >  NC_000913/1212099‑1212351

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: