Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,531,596 (G)6→5 coding (134/762 nt) cysZ → sulfate transporter, sulfite inhibited

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,531,5910G.93.8% 276.8 / 12.7 65coding (129/762 nt)cysZsulfate transporter, sulfite inhibited
Reads supporting (aligned to +/- strand):  ref base G (2/2);  new base . (23/38);  total (25/40)
Fisher's exact test for biased strand distribution p-value = 6.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

CAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAGTTATCT  >  NC_000913/2531492‑2531686
                                                                                                   |                                                                                               
cAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATg                                                                                                 <  2:1683694/100‑1 (MQ=255)
cAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATg                                                                                                 <  1:1904227/100‑1 (MQ=255)
 aGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATgg                                                                                                >  1:60768/1‑100 (MQ=255)
  gCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCCATATTTTGTTGATggg                                                                                               <  1:2085709/100‑1 (MQ=255)
   cGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATgggg                                                                                              >  2:1321173/1‑100 (MQ=255)
     gTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATgggggc                                                                                            <  2:2058217/100‑2 (MQ=255)
            cTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGCGCATTCt                                                                                    >  1:1117713/1‑100 (MQ=255)
            cTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGCGCATTCt                                                                                    >  1:654126/1‑100 (MQ=255)
            cTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGCGCATTCt                                                                                    >  1:1540268/1‑100 (MQ=255)
              aTTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGCGCATTCtgg                                                                                  <  2:674543/100‑1 (MQ=255)
              aTTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGCGCATTCtgg                                                                                  >  2:1025207/1‑100 (MQ=255)
               ttttGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGCGCATTCtggt                                                                                 >  1:357975/1‑100 (MQ=255)
                tttGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGCGCATTCtggtg                                                                                <  1:761750/100‑1 (MQ=255)
                tttGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGCGCATTCtggtg                                                                                <  2:174537/100‑1 (MQ=255)
                tttGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGCGCATTCtggtg                                                                                <  2:572926/100‑1 (MQ=255)
                 ttGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCtggtg                                                                                <  2:412105/100‑1 (MQ=255)
                    cgcAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGCGCATTCTGGTGGctc                                                                            >  1:73642/1‑100 (MQ=255)
                     gcAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGCGCATTCTGGTGGctct                                                                           >  1:1844064/1‑100 (MQ=255)
                      cAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGCGCATTCTGGTGGCTCtt                                                                          <  2:1660138/100‑1 (MQ=255)
                      cAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGCGCATTCTGGTGGCTCtt                                                                          <  2:1144176/100‑1 (MQ=255)
                        aGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGCGCATTCTGGTGGCTCTTTa                                                                        <  1:281917/100‑1 (MQ=255)
                         ggCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGCGCATTCTGGTGGCTCTTTac                                                                       <  1:418943/100‑1 (MQ=255)
                           cTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGCGCATTCTGGTGGCTCTTTacac                                                                     <  1:76864/100‑1 (MQ=255)
                           cTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGCGCATTCTGGTGGCTCTTTacac                                                                     <  1:1943739/100‑1 (MQ=255)
                            tGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTTGTCAATATTTTGTTGATGGGGGCGCATTCTGGTGGCTCTTTAcaca                                                                    <  1:1239230/100‑1 (MQ=255)
                            tGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTacac                                                                     <  1:1172685/100‑1 (MQ=255)
                            tGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGCGCATTCTGGTGGCTCTTTAcaca                                                                    <  2:1353290/100‑1 (MQ=255)
                             ggAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGCGCATTCTGGTGGCTCTTTAcaca                                                                    <  2:949955/99‑1 (MQ=255)
                             ggAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGCGCATTCTGGTGGCTCTTTACACAg                                                                   >  2:796818/1‑100 (MQ=255)
                                aGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGCGCATTCTGGTGGCTCTTTACACAGCTc                                                                <  2:1117363/100‑1 (MQ=255)
                                 gCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCg                                                               >  1:1026870/1‑100 (MQ=255)
                                  cTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGa                                                              >  2:798511/1‑100 (MQ=255)
                                   tCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGAt                                                             >  1:585125/1‑100 (MQ=255)
                                    cGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATg                                                            <  2:156745/100‑1 (MQ=255)
                                    cGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATg                                                            <  2:156850/100‑1 (MQ=255)
                                      tctcGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTc                                                          >  1:1847315/1‑100 (MQ=255)
                                        tcGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTg                                                        <  2:107857/100‑1 (MQ=255)
                                        tcGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTg                                                        <  1:1729362/100‑1 (MQ=255)
                                        tcGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTg                                                        <  2:357976/100‑1 (MQ=255)
                                         cGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTg                                                        >  1:1438451/1‑100 (MQ=255)
                                          gCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGCGCAGTCTGGTGGCTCTTTACACAGCTCGATGTCTGGa                                                      >  2:1928619/1‑100 (MQ=255)
                                           cAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGAt                                                     <  2:804880/100‑1 (MQ=255)
                                           cAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGAt                                                     <  1:1825500/100‑1 (MQ=255)
                                              ccTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATccc                                                  <  1:1505371/100‑1 (MQ=255)
                                              ccTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATccc                                                  <  1:514965/100‑1 (MQ=255)
                                               cTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCg                                                 >  1:1112647/1‑100 (MQ=255)
                                               cTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCg                                                 >  2:1607793/1‑100 (MQ=255)
                                                      tCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCa                                          <  2:1250326/100‑1 (MQ=255)
                                                        ggCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGCGCATTCTGGTGGCTCTTTACTCAGCTCGATGTCTGGATCCCGACTCTCATg                                        <  2:60768/100‑1 (MQ=255)
                                                        ggCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATg                                        <  2:1695028/100‑1 (MQ=255)
                                                        ggCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATg                                        >  1:554907/1‑100 (MQ=255)
                                                            tttCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGtt                                    >  2:1457652/1‑100 (MQ=255)
                                                               cGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACg                                 >  2:287068/1‑100 (MQ=255)
                                                                attATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGt                                <  2:1529214/99‑1 (MQ=255)
                                                                    ttttACCGCTGCTGGTCAATATTTTGTTGATGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCg                            >  1:97928/1‑100 (MQ=255)
                                                                    ttttACCGCTGCTGGTCAATATTTTGTTGATGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCg                            <  2:588928/100‑1 (MQ=255)
                                                                       tACCGCTGCTGGTCAATATTTTGTTGATGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGAc                         >  1:287097/1‑100 (MQ=255)
                                                                           gctgctGGTCAATATTTTGTTGATGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGc                     <  2:899341/100‑1 (MQ=255)
                                                                            ctgctgGTCAATATTTTGTTGATGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCt                    >  2:1470581/1‑100 (MQ=255)
                                                                              gctgGTCAATATTTTGTTGATGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTAc                  <  1:288839/100‑1 (MQ=255)
                                                                              gctgGTCAATATTTTGTTGATGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTAc                  <  2:903519/100‑1 (MQ=255)
                                                                                    cAATATTTTGTTGATGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGc            <  2:939292/100‑1 (MQ=255)
                                                                                     aaTATTTTGTTGATGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCt           <  2:1283253/100‑1 (MQ=255)
                                                                                     aaTATTTTGTTGATGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCt           <  2:301828/100‑1 (MQ=255)
                                                                                      atatTTTGTTGATGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTg          >  1:843963/1‑100 (MQ=255)
                                                                                       tatTTTGTTGATGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGa         <  1:1400300/100‑1 (MQ=255)
                                                                                         ttttGTTGATGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAGt       <  1:992228/100‑1 (MQ=255)
                                                                                           ttgttgATGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAGTTa     >  2:205310/1‑100 (MQ=255)
                                                                                             gttgATGGGGGCGCGTTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAGTTATc   <  1:502921/100‑1 (MQ=255)
                                                                                             gttgATGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAGTTATc   <  2:871140/100‑1 (MQ=255)
                                                                                               tgATGGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCGTGAGTTACGTTCCGGACTGGCTACAATGGCTGAGTTATCt  >  2:723286/1‑100 (MQ=255)
                                                                                                   |                                                                                               
CAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAGTTATCT  >  NC_000913/2531492‑2531686

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