Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,860,426 Δ1 bp intergenic (+10/+41) pphB → / ← ygbI serine/threonine‑specific protein phosphatase 2/DeoR family putative transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,860,4250T.95.7% 84.8 / ‑2.7 23intergenic (+9/+42)pphB/ygbIserine/threonine‑specific protein phosphatase 2/DeoR family putative transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base . (11/11);  total (11/12)
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

ATAAACGGCGCTGATTATTTTATATTTGGACATATGATGTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAGTCTCATGCTTCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAGCGCGCGACTGGCACTCTCCGGCAAACCGTCGTCT  >  NC_000913/2860288‑2860545
                                                                                                                                         |                                                                                                                        
aTAAACGGCGCTGATTATTTTATATTTGGACATATGATGTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGctc                                                                                                                         >  2:118770/1‑138 (MQ=255)
   aaCGGCGCTGATTATTTTATATTTGGACATATGATGTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGctcttc                                                                                                                      <  1:118770/139‑5 (MQ=255)
       gcgcTGATTATTTTATATTTGGACATATGATGTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGTTCTTCtgtg                                                                                                                 <  2:536988/139‑1 (MQ=255)
          cTGATTAATTTATATTTGAACATATGATGTTTGATAACATTCAGACGTTCGCTAACTAGATTTATATTGATACCGGAGCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAg                                                                                                              <  1:352470/139‑1 (MQ=255)
             attattTTATATTTGGACATATGATGTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGCAt                                                                                                           <  2:261976/139‑1 (MQ=255)
                     atatTTGGACATATGATGTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGCATGAGTAAcc                                                                                                   <  2:320904/139‑1 (MQ=255)
                       atTTGGACATATGATGTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGCATGAGTAACCCg                                                                                                 >  1:631931/1‑139 (MQ=255)
                          tGGACATATGATGTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGCATGAGTAACCCGgtg                                                                                              <  2:510951/139‑1 (MQ=255)
                           ggACATATGATGTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGCATGAGTAACCCGgtgt                                                                                             >  1:672713/1‑139 (MQ=255)
                                tatGATGTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGCATGAGTAACCCGGTGTTATTg                                                                                        >  2:134581/1‑139 (MQ=255)
                                           aTAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCattat                                                                             >  2:773358/1‑139 (MQ=255)
                                                 ttCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCAttt                                                                       >  1:512450/1‑139 (MQ=255)
                                                      aCGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGc                                                                  >  1:677584/1‑139 (MQ=255)
                                                        gTTCGCTAACCAGATTTATCTTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTa                                                                <  1:523209/139‑1 (MQ=255)
                                                             cTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGc                                                           <  2:327564/139‑1 (MQ=255)
                                                                aCCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCaaa                                                        <  1:134581/139‑1 (MQ=255)
                                                                               aTACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGCATGAGt                                                                                                       >  1:335719/1‑77 (MQ=255)
                                                                               aTACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGCATGAGt                                                                                                       <  2:335719/77‑1 (MQ=255)
                                                                                tACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGct                                                                                                                          <  2:144535/58‑1 (MQ=255)
                                                                                tACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGct                                                                                                                          >  1:144535/1‑58 (MQ=255)
                                                                                      aTCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAgcg                                  >  2:54885/1‑139 (MQ=255)
                                                                                                  cgggcggCTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAGCGCGCGACTGGCAc                      >  1:12813/1‑139 (MQ=255)
                                                                                                  cgggcggCTGTCATTTTATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAGCGCGCGACTGGCAc                      >  1:282667/1‑139 (MQ=255)
                                                                                                                tttATAAAATAAAGTAATCTCATGCTCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAGCGCGCGACTGGCACTCTCCGGCAAACcg        <  2:12813/139‑1 (MQ=255)
                                                                                                                  tataAAATAAAGTAATCTCATGCTCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAGCGCGCGACTGGCACTCTCCGGCAAACcgtc      <  2:672713/139‑1 (MQ=255)
                                                                                                                    taaaataaaGTAATCTCATGCTCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAGCGCGCGACTGGCACTCTCCGGCAAACcgtcgt    <  1:732418/139‑1 (MQ=255)
                                                                                                                      aaataaaGTAATCTCATGCTCTTCTGTGAAGCATGAGTAACCCGGCGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAGCGCGCGACTGGCACTCTCCGGCAAACCGTCGTCt  >  1:427420/1‑139 (MQ=255)
                                                                                                                                         |                                                                                                                        
ATAAACGGCGCTGATTATTTTATATTTGGACATATGATGTTTGATAACATTCAGACGTTCGCTAACCAGATTTATATTGATACCGGATCGCCGAACAGCGGGCGGCTGTCATTTTATAAAATAAAGTAGTCTCATGCTTCTTCTGTGAAGCATGAGTAACCCGGTGTTATTGCAGGCCATTATTCATTTTTCGCTACCAGCAAAGAGAGATCCTGCTTCACCAGCGCGCGACTGGCACTCTCCGGCAAACCGTCGTCT  >  NC_000913/2860288‑2860545

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