Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,370,044 (T)9→8 intergenic (+41/‑172) pspE → / → ycjM thiosulfate sulfurtransferase PspE/glucosylglycerate phosphorylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,370,0360T.92.3% 87.0 / 2.6 26intergenic (+33/‑180)pspE/ycjMthiosulfate sulfurtransferase PspE/glucosylglycerate phosphorylase
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base . (9/15);  total (10/16)
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

CAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAAT  >  NC_000913/1369915‑1370152
                                                                                                                         |                                                                                                                    
caaAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGttttttttctgt                                                                                                           <  2:449538/133‑5 (MQ=255)
   aggagaTCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTcac                                                                                                       <  1:302613/131‑1 (MQ=255)
   aaGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTcac                                                                                                       <  2:74440/133‑1 (MQ=255)
   aaGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTcac                                                                                                       <  1:31529/133‑1 (MQ=255)
     gagaTCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCacac                                                                                                     <  1:332173/133‑1 (MQ=255)
          ccTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTAt                                                                                                >  2:392527/1‑133 (MQ=255)
               gCGAGATGGGTTATAACCACGTTGAGAACGCCCGTGGCCTGCAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACa                                                                                           <  2:147940/133‑1 (MQ=255)
                     tGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCa                                                                                     >  1:32774/1‑133 (MQ=255)
                      gggATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCAt                                                                                    >  2:363423/1‑133 (MQ=255)
                              cccACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGattatt                                                                            >  1:301716/1‑133 (MQ=255)
                               ccACGTTGAGAACGCCGGTGGCCTGAAAGACATCACAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTAttt                                                                           >  2:445896/1‑133 (MQ=255)
                                 aCGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGc                                                                         >  1:337396/1‑133 (MQ=255)
                                   gTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGt                                                                       <  1:450359/133‑1 (MQ=255)
                                             ccGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACaa                                                             >  1:402188/1‑133 (MQ=255)
                                              cGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACAACCCCATAGATTATTCGCGTCAGCTCACaaa                                                            >  1:195534/1‑133 (MQ=255)
                                                        aaaGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCttttt                                                  >  2:414900/1‑133 (MQ=255)
                                                                  cAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTaaa                                        <  1:445896/133‑1 (MQ=255)
                                                                    aTGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTaaaaa                                      <  2:330385/133‑1 (MQ=255)
                                                                                 aaGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCg                         <  1:392527/133‑1 (MQ=255)
                                                                                  aGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGt                        <  1:266696/133‑1 (MQ=255)
                                                                                       aaaaGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACt                    >  1:394133/1‑133 (MQ=255)
                                                                                               tCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCtt           <  2:301716/133‑1 (MQ=255)
                                                                                                   cATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGc       <  1:391783/133‑1 (MQ=255)
                                                                                                    aTGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGc       <  2:394133/133‑1 (MQ=255)
                                                                                                    aTGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCt      <  1:414900/133‑1 (MQ=255)
                                                                                                        gcggcgCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAAt  <  2:337396/133‑1 (MQ=255)
                                                                                                        gcggcgCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAAt  <  2:402188/133‑1 (MQ=255)
                                                                                                                         |                                                                                                                    
CAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAAT  >  NC_000913/1369915‑1370152

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