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:cyanide sulfurtransferase (rhodanese)/alpha amylase catalytic domain family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,370,0360T.97.1% 127.8 / ‑1.9 34intergenic (+33/‑180)pspE/ycjMthiosulfate:cyanide sulfurtransferase (rhodanese)/alpha amylase catalytic domain family protein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base . (16/17);  total (16/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.38e-01

CAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTG  >  NC_000913/1369902‑1370169
                                                                                                                                      |                                                                                                                                     
cAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGttttt                                                                                                                                   >  1:1228045/1‑139 (MQ=255)
  gTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGttttttt                                                                                                                                 >  1:140671/1‑139 (MQ=255)
  gTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGttttttt                                                                                                                                 >  1:1431761/1‑139 (MQ=255)
          aagcaaAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTca                                                                                                                         <  2:127030/139‑1 (MQ=255)
             caaAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACAcc                                                                                                                     <  2:1431761/139‑1 (MQ=255)
                  gagaTCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTAtt                                                                                                                <  2:1228045/139‑1 (MQ=255)
                                 aTGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGat                                                                                                 >  2:129454/1‑139 (MQ=255)
                                          aCCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGt                                                                                        <  1:129454/139‑1 (MQ=255)
                                          aCCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGt                                                                                        <  2:140671/139‑1 (MQ=255)
                                               cGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGtttttttt                                                                                                                                <  2:791544/95‑1 (MQ=255)
                                               cGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGtttttttt                                                                                                                                >  1:791544/1‑95 (MQ=255)
                                                   gagaACGCCGGTGGCCTGAAAGACTTCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACAcc                                                                                                                     >  1:1094720/1‑101 (MQ=255)
                                                   gagaACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGACG‑TTTTTTTTCTGTCACAcc                                                                                                                     <  2:1094720/101‑1 (MQ=255)
                                                      aaCGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAAt                                                                            >  2:157621/1‑139 (MQ=255)
                                                         gCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACAcc                                                                                                                     >  1:751106/1‑95 (MQ=255)
                                                         gCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACAcc                                                                                                                     <  2:751106/95‑1 (MQ=255)
                                                                            tCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCa                                                                                                      <  1:1471684/91‑1 (MQ=255)
                                                                            tCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCa                                                                                                      >  2:1471684/1‑91 (MQ=255)
                                                                                       aaGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGc                                           >  2:842706/1‑139 (MQ=255)
                                                                                              aaGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATc                                                                                                          <  1:1451872/69‑1 (MQ=255)
                                                                                              aaGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATc                                                                                                          >  2:1451872/1‑69 (MQ=255)
                                                                                              aaGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGa                                                                                                  <  1:7591/77‑1 (MQ=255)
                                                                                              aaGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGa                                                                                                  >  2:7591/1‑77 (MQ=255)
                                                                                                    aaaaGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACa                                                                                                            <  1:1491633/61‑1 (MQ=255)
                                                                                                    aaaaGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACa                                                                                                            >  2:1491633/1‑61 (MQ=255)
                                                                                                    aaaaGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAAccc                              <  1:842706/139‑1 (MQ=255)
                                                                                                     aaaGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTAc                                                                                                             >  1:684524/1‑59 (MQ=255)
                                                                                                     aaaGAATTCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTAc                                                                                                             <  2:684524/59‑1 (MQ=255)
                                                                                                           ttCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTca                                                                                 <  2:420859/81‑1 (MQ=255)
                                                                                                           ttCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTca                                                                                 >  1:420859/1‑81 (MQ=255)
                                                                                                           ttCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAAt                                                                            >  2:470487/1‑86 (MQ=255)
                                                                                                           ttCACCATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAAt                                                                            <  1:470487/86‑1 (MQ=255)
                                                                                                                cATGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAAtt                  >  2:726111/1‑139 (MQ=255)
                                                                                                                  tGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTcac                                                                                                                        <  1:1425563/35‑1 (MQ=255)
                                                                                                                  tGAGCGGCGCTTATGCGCCG‑TTTTTTTTCTGTcac                                                                                                                        >  2:1425563/1‑35 (MQ=255)
                                                                                                                         cgcTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGattat                                                                                              >  2:402823/1‑54 (MQ=255)
                                                                                                                         cgcTTATGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGattat                                                                                              <  1:402823/54‑1 (MQ=255)
                                                                                                                              aTGCGCCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTc    >  2:1227634/1‑139 (MQ=255)
                                                                                                                                gcgcCG‑TTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTg  <  1:1227634/139‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                     
CAGTCAGGGCAAGCAAAAGAGATCCTTAGCGAGATGGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTGACTAAAACCCTTGTGCTCAATTGACAGTTTATTTTCTG  >  NC_000913/1369902‑1370169

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