Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,420,727 G→A G36D (GGC→GAC)  yfcG → GSH‑dependent disulfide bond oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,420,7270GA97.3% 248.4 / ‑4.1 73G36D (GGC→GAC) yfcGGSH‑dependent disulfide bond oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base A (41/30);  total (41/32)
Fisher's exact test for biased strand distribution p-value = 1.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCTATTTCGCCCCGACACCCAATGGTCACAAAATTACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGGCGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTCCTGCCGATGGCGGCGAACCGCTA  >  NC_000913/2420631‑2420821
                                                                                                |                                                                                              
tctATTTCGCCCCGACACCCAATGGTCACAAAATTACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACgg                                                                                             <  2:329849/100‑1 (MQ=255)
  tATTTCGCCCCGACACCCAATGGTCACAAAATTACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGGCGGTc                                                                                           <  1:58800/100‑1 (MQ=255)
  tATTTCGCCCCGACACCCAATGGTCACAAAATTACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTc                                                                                           <  2:1129324/100‑1 (MQ=255)
       cGCCCCGACACCCAATGGTCACAAAATTACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGttt                                                                                      <  2:598186/100‑1 (MQ=255)
        gCCCCGACACCCAATGGTCACAAAATTACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGttt                                                                                      <  2:247403/99‑1 (MQ=255)
        gCCCCGACACCCAATGGTCACAAAATTACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTc                                                                                     >  1:503161/1‑100 (MQ=255)
            cGACACCCAATGGTCACAAAATTACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGAAAAGACGGTCAGTTTCGTcc                                                                                 <  2:1356364/100‑1 (MQ=255)
                 cccAATGGTCACAAAATTACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGaa                                                                             >  2:943924/1‑99 (MQ=255)
                 cccAATGGTCACAAAATTACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAAt                                                                            >  2:784031/1‑100 (MQ=255)
                 cccAATGGTCACAAAATTACACTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAAt                                                                            >  2:567139/1‑100 (MQ=255)
                  ccAATGGTCACAAAATTACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAAtt                                                                           >  2:613093/1‑100 (MQ=255)
                        gTCACAAAATTACACTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTgcg                                                                     <  1:567140/100‑1 (MQ=255)
                          cacaAAATTACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCa                                                                   >  2:588592/1‑100 (MQ=255)
                            caAAATTACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCAtt                                                                 >  2:1412434/1‑100 (MQ=255)
                            caAAATTACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCAtt                                                                 >  1:192987/1‑100 (MQ=255)
                            caAAATTACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCAtt                                                                 >  2:1417482/1‑100 (MQ=255)
                             aaaaTTACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCAttt                                                                <  1:473189/100‑1 (MQ=255)
                                aTTACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGc                                                             >  1:1566062/1‑100 (MQ=255)
                                  tACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGccc                                                           <  1:1262449/100‑1 (MQ=255)
                                  tACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGccc                                                           >  2:1206995/1‑100 (MQ=255)
                                   aCGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGccc                                                           <  1:1422738/99‑1 (MQ=255)
                                    cGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCcaa                                                         >  2:637126/1‑100 (MQ=255)
                                    cGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCcaa                                                         >  2:45108/1‑100 (MQ=255)
                                    cGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCca                                                          <  2:1048771/99‑1 (MQ=255)
                                    cGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCACCca                                                          <  2:1137957/99‑1 (MQ=255)
                                     gCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCcaac                                                        <  2:103755/100‑1 (MQ=255)
                                      cTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCcaaca                                                       <  1:1417490/100‑1 (MQ=255)
                                       tGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCcaacaa                                                      >  1:1135911/1‑100 (MQ=255)
                                       tGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCcaacaa                                                      >  1:232619/1‑100 (MQ=255)
                                            ctcGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACaaaa                                                 >  1:842983/1‑100 (MQ=255)
                                             tcGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAAt                                                <  1:640691/100‑1 (MQ=255)
                                             tcGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAAt                                                >  2:1022158/1‑100 (MQ=255)
                                              cGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAAtt                                               >  2:1143728/1‑100 (MQ=255)
                                              cGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAAtt                                               >  1:731846/1‑100 (MQ=255)
                                              cGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAAtt                                               >  2:1296998/1‑100 (MQ=255)
                                                 agaagCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCg                                            <  1:1027447/100‑1 (MQ=255)
                                                 agaagCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCg                                            >  2:604013/1‑100 (MQ=255)
                                                 agaagCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCg                                            >  2:894983/1‑100 (MQ=255)
                                                  gaagCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCgg                                           >  2:1162375/1‑100 (MQ=255)
                                                    agCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCa                                         >  1:214942/1‑100 (MQ=255)
                                                     gCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCaa                                        >  1:1294697/1‑100 (MQ=255)
                                                      cAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAAt                                       >  1:1411460/1‑100 (MQ=255)
                                                          gCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAAttgtt                                   >  1:620239/1‑100 (MQ=255)
                                                           cTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAAttgttg                                  <  1:1451273/100‑1 (MQ=255)
                                                               aTTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCa                              >  1:1372127/1‑100 (MQ=255)
                                                                ttATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCAt                             >  2:340462/1‑100 (MQ=255)
                                                                 tATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCAtt                            >  1:1144687/1‑100 (MQ=255)
                                                                  aTCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATtc                           <  2:115165/100‑1 (MQ=255)
                                                                  aTCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATtc                           <  2:214942/100‑1 (MQ=255)
                                                                  aTCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATtc                           >  1:347813/1‑100 (MQ=255)
                                                                    cGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATtctc                         >  2:673227/1‑100 (MQ=255)
                                                                    cGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTACGGCAATTGTTGATCATtctc                         >  2:1405799/1‑100 (MQ=255)
                                                                     gCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAATATTCCGGCAATTGTTGATCATTCTcc                        <  1:960891/100‑1 (MQ=255)
                                                                     gCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTcc                        >  1:1414445/1‑100 (MQ=255)
                                                                       ttGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTCCTg                      <  2:231112/100‑1 (MQ=255)
                                                                        tGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTCCTGc                     >  2:1058204/1‑100 (MQ=255)
                                                                        tGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTCCTGc                     <  1:309264/100‑1 (MQ=255)
                                                                         gATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTCCTGcc                    <  1:894988/100‑1 (MQ=255)
                                                                          aTTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTCCTGCCg                   <  2:75113/100‑1 (MQ=255)
                                                                           ttAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTCCTGCCGa                  >  1:1265488/1‑100 (MQ=255)
                                                                             aaGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTCCTGCCGATg                <  2:59199/100‑1 (MQ=255)
                                                                              aGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTCCTGCCGATgg               <  1:1162920/100‑1 (MQ=255)
                                                                                gTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTCCTGCCGATggcg             <  2:1144682/100‑1 (MQ=255)
                                                                                gTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTCCTGCCGATggcg             <  2:503160/100‑1 (MQ=255)
                                                                                 tAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTCCTGCCGATggcgg            >  1:1448290/1‑100 (MQ=255)
                                                                                  aGACCTGGGGAAAGGCGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAGAATTCCGGCAATTGTTGATCAGTCTCCTGCCGATggcggc           <  2:1502000/100‑1 (MQ=255)
                                                                                   gACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTCCTGCCGATGgcggcg          >  1:864742/1‑100 (MQ=255)
                                                                                     ccTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTCCTGCCGATGGCGGCGaa        >  2:58893/1‑100 (MQ=255)
                                                                                       tGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTCCTGCCGATGGCGGCGAAcc      <  2:395680/100‑1 (MQ=255)
                                                                                       tGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTCCTGCCGATGGCGGCGAAcc      >  2:432717/1‑100 (MQ=255)
                                                                                        ggggAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTCCTGCCGATGGCGGCGAACCg     >  2:274136/1‑100 (MQ=255)
                                                                                          ggAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTCCTGCCGATGGCGGCGAACCGCt   <  2:242881/100‑1 (MQ=255)
                                                                                           gAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTCCTGCCGATGGCGGCGAACCGCTa  <  2:1029747/100‑1 (MQ=255)
                                                                                                |                                                                                              
TCTATTTCGCCCCGACACCCAATGGTCACAAAATTACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGGCGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTCCTGCCGATGGCGGCGAACCGCTA  >  NC_000913/2420631‑2420821

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: ‑