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.89.7% 202.7 / 22.7 59intergenic (+33/‑180)pspE/ycjMthiosulfate:cyanide sulfurtransferase (rhodanese)/alpha amylase catalytic domain family protein
Reads supporting (aligned to +/- strand):  ref base T (4/2);  new base . (28/24);  total (33/26)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTG  >  NC_000913/1369937‑1370131
                                                                                                   |                                                                                                
gggATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGt                                                                                                  >  2:1303238/1‑100 (MQ=255)
 ggATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAAAATTCACCATGAGCGGCGCTTATGCGCCGtt                                                                                                 >  2:515838/1‑100 (MQ=255)
        cccACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGttttttttc                                                                                         <  2:1326777/100‑2 (MQ=255)
         ccACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGttttttttct                                                                                        <  1:920444/100‑3 (MQ=255)
           aCGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGttttttttctgt                                                                                      <  2:940913/100‑5 (MQ=255)
            cGTTGAGCACTCCGGTGGCATGAAAGACATCGCAATGCCGAAAGTCAAAAGTTAAAAGAATGCACCATGAGCGGCGCTTATGTGCCGTTTTTTTTCTGTc                                                                                    <  1:444155/100‑1 (MQ=255)
            cGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGttttttttct                                                                                        >  2:66798/1‑95 (MQ=255)
            cGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTc                                                                                    >  2:915866/1‑100 (MQ=255)
            cGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTc                                                                                    <  2:480004/100‑1 (MQ=255)
                gagaACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATTCGTCGTTTTTTTTCTGTCacac                                                                                <  1:393343/100‑1 (MQ=255)
                 agaACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACAcc                                                                               <  2:279595/100‑1 (MQ=255)
                 agaACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACAcc                                                                               >  2:1352790/1‑100 (MQ=255)
                     cGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTAt                                                                           <  1:1303228/100‑1 (MQ=255)
                     atctggTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTa                                                                            <  1:295777/96‑1 (MQ=255)
                      gCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTAtt                                                                          <  1:128268/100‑1 (MQ=255)
                        cGGTGGCCCGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATTAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTAtttn                                                                        <  1:132951/100‑2 (MQ=255)
                         gtttgCCTTAAAGATATTGGAATGCCGAAGGTCAATGGTTATAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTAtt                                                                          <  2:269326/93‑1 (MQ=255)
                         ggTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTaaga                                                                         >  2:295779/1‑96 (MQ=255)
                          gTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACa                                                                      <  2:924608/100‑1 (MQ=255)
                             gCCTGCAAGACATAGCCATGCCGAAGGTCAAAGGTTCAAAGAATTCAACATGCGCGGCGCTTATGCGCCGTTTTTTTTTCTGTTAAAAATTAtttttatc                                                                    >  2:360109/1‑95 (MQ=255)
                              ccTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATcc                                                                   >  1:240793/1‑99 (MQ=255)
                               cTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATcccc                                                                 >  1:1614513/1‑100 (MQ=255)
                                      aCATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGatt                                                          <  1:770482/100‑1 (MQ=255)
                                        aTCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGcttat                                                        >  1:702795/1‑100 (MQ=255)
                                        aTCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGattat                                                        >  2:53251/1‑100 (MQ=255)
                                         tCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGattatt                                                       >  1:1793519/1‑100 (MQ=255)
                                            cAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGattatt                                                       >  2:440991/1‑97 (MQ=255)
                                               tGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTc                                                 >  2:1462245/1‑100 (MQ=255)
                                               tGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTc                                                 <  1:255307/100‑1 (MQ=255)
                                               tGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCAAACCTTATTTACATCCCCATAGATTATTTGCGTc                                                 >  1:63286/1‑100 (MQ=255)
                                                gCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCa                                                >  2:167495/1‑100 (MQ=255)
                                                 ccGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAg                                               >  2:929539/1‑100 (MQ=255)
                                                  cGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTGATTTACATGCCCATAGATTATTTGCGTCAGc                                              >  1:1223029/1‑100 (MQ=255)
                                                   gAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCt                                             <  2:813477/100‑1 (MQ=255)
                                                     aGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCGCACCTTATTTACATCCCCATAGATTATTTGCGTCAGc                                              >  1:1340316/1‑97 (MQ=255)
                                                      ggTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTcac                                          >  1:457679/1‑100 (MQ=255)
                                                      ggTCAAAGGTTAAAAGAATTCACCATGAGCGGCGATTATGCGCCGTTTTTTTTCTGTCACACCTTATTTCCCTCCCCATAGATAATTTGCGTCAGCTcac                                          >  1:1710971/1‑100 (MQ=255)
                                                       gTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTcaca                                         <  2:80641/100‑1 (MQ=255)
                                                       gTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGATTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTcaca                                         >  1:634890/1‑100 (MQ=255)
                                                         cAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACaaa                                       >  2:609484/1‑100 (MQ=255)
                                                         cAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACaaa                                       >  2:1455910/1‑100 (MQ=255)
                                                           aaGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATa                                     <  2:791886/100‑1 (MQ=255)
                                                             ggTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACg                                   >  2:762750/1‑100 (MQ=255)
                                                               ttAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCt                                 <  1:304476/100‑1 (MQ=255)
                                                                 aaaaGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCttt                               <  1:736288/100‑1 (MQ=255)
                                                                  aaaGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCtttt                              >  2:184148/1‑100 (MQ=255)
                                                                  aaaGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCtttt                              >  2:967602/1‑100 (MQ=255)
                                                                   aaGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCttttt                             <  2:11204/100‑1 (MQ=255)
                                                                    aGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTc                            <  2:1793532/100‑1 (MQ=255)
                                                                         agaaaaaaaacggcgCATAAGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATCATTTGCGTCAGCTCACAAATACGCTTTTTCCCTgg                       >  1:556410/11‑100 (MQ=255)
                                                                           aCCATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACTCCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTa                     >  1:927476/1‑100 (MQ=255)
                                                                            ccATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTaa                    >  1:752572/1‑100 (MQ=255)
                                                                            ccATGAGCGGCGCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTaa                    >  2:1506368/1‑100 (MQ=255)
                                                                             cATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTaa                    >  2:1257516/1‑100 (MQ=255)
                                                                                 aacggcgCATAAGCGCCGCTTTTTTTTTTTTTTTTCCTTATTTACATCCCCATATATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTTGTaaaaa                 >  1:1775962/3‑99 (MQ=255)
                                                                                 aGCGGCGCTTATGCGCCGTTTTTCTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTaaaaaa                <  1:1064925/99‑1 (MQ=255)
                                                                                 aGCGGCGCATAAGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTTCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAAt               >  2:1685545/1‑100 (MQ=255)
                                                                                  gcggcgCTTATGCGCCGTTTTTTTTCTGTTACACCTTATTTTCATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAAt               <  2:136957/99‑1 (MQ=255)
                                                                                   cggcgCTTATGCGCCGTTTACTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAAt               <  2:800630/100‑1 (MQ=255)
                                                                                   cggcgCATAAGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAAt               >  2:1638037/1‑99 (MQ=255)
                                                                                   ccgcgcTTATGCGCCGTTTTTTTTCTTTCACGACTTATTTACATCCCCATCGATTATTTGCGTCAGCTCACATATGCGCTTTTTCCCTGGTAAAAAAtt              <  1:677891/97‑2 (MQ=255)
                                                                                    ggcgCTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGAt            <  1:1363298/100‑1 (MQ=255)
                                                                                     gtgcTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGAt            <  1:1353654/97‑1 (MQ=255)
                                                                                       gcTTATGCGCCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGAttt          >  2:1261518/1‑99 (MQ=255)
                                                                                            atgcgcttttttttttCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGc     <  2:1804353/93‑1 (MQ=255)
                                                                                              cgcCGTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTg  <  2:36630/100‑1 (MQ=255)
                                                                                                   |                                                                                                
GGGATATACCCACGTTGAGAACGCCGGTGGCCTGAAAGACATCGCAATGCCGAAGGTCAAAGGTTAAAAGAATTCACCATGAGCGGCGCTTATGCGCCGTTTTTTTTTCTGTCACACCTTATTTACATCCCCATAGATTATTTGCGTCAGCTCACAAATACGCTTTTTCCCTGGTAAAAAATGATTTCCTGCGTG  >  NC_000913/1369937‑1370131

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