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,7270GA94.1% 205.6 / 3.2 69G36D (GGC→GAC) yfcGGSH‑dependent disulfide bond oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (2/2);  new base A (35/29);  total (37/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.62e-01

CTCTATTTCGCCCCGACACCCAATGGTCACAAAATTACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGGCGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTCCTGCCGATGGCGGCGAACCGCTAAG  >  NC_000913/2420630‑2420823
                                                                                                 |                                                                                                
ctctATTTCGCCCCGACACCCAATGGTCACAAAATTACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACg                                                                                                <  2:1256958/100‑1 (MQ=255)
 tctATTTCGCCCCGACACCCAATGGTCACAAAATTACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACgg                                                                                               >  2:1655288/1‑100 (MQ=255)
  ctATTTCGCCCCGACACCCAATGGTCACAAAATTACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGt                                                                                              >  1:130787/1‑100 (MQ=255)
     tttCGCCCCGACACCCAATGGTCACAAAATTACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAg                                                                                           >  2:1831714/1‑100 (MQ=255)
      ttCGCCCCGACACCCAATGGTCACAAAATTACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGt                                                                                          >  1:800822/1‑100 (MQ=255)
          ccccGACACCCAATGGTCACAAAATTACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCg                                                                                      <  2:1357937/100‑1 (MQ=255)
          ccccGACACCCAATGGTCACAAAATTACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCg                                                                                      >  2:800036/1‑100 (MQ=255)
          ccccGACACCCAATGGTCACAAAATTACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCg                                                                                      <  1:1802041/100‑1 (MQ=255)
            ccGACACCCAATGGTCACAAAATTACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTc                                                                                    >  1:76087/1‑100 (MQ=255)
            ccGACACCCAATGGTCACAAAATTACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTc                                                                                    <  2:1789335/100‑1 (MQ=255)
                  cccAATGGTCACAAAATTACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAAt                                                                              >  2:606606/1‑100 (MQ=255)
                   ccAATGGTCACAAAATTACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAAtt                                                                             >  2:886179/1‑100 (MQ=255)
                   ccAATGGTCACAAAATTACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAAtt                                                                             >  2:510795/1‑100 (MQ=255)
                             caAAATTACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCAtt                                                                   >  2:29638/1‑100 (MQ=255)
                                  ttACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGGCGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGcc                                                              <  1:1306711/100‑1 (MQ=255)
                                  ttACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGGCGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGcc                                                              >  1:611211/1‑100 (MQ=255)
                                  ttAAGCTTTTTTTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGcc                                                              <  1:1743140/100‑1 (MQ=255)
                                   tACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGccc                                                             <  2:829357/100‑1 (MQ=255)
                                   tACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGccc                                                             <  2:696648/100‑1 (MQ=255)
                                     cGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCcaa                                                           >  2:1040127/1‑100 (MQ=255)
                                       cTGTTTCTCGCAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGCTTCGTCCGGAATTTTTGCGCATTTCGCCcaaca                                                         <  1:1279289/100‑1 (MQ=255)
                                        tGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCcaacaa                                                        >  2:1184428/1‑100 (MQ=255)
                                         gTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCcaacaac                                                       <  2:168846/100‑1 (MQ=255)
                                             ctcGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACaaaa                                                   >  1:1453986/1‑100 (MQ=255)
                                               cGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAAtt                                                 >  1:235193/1‑100 (MQ=255)
                                                gaagaaGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTc                                                >  2:474560/1‑100 (MQ=255)
                                                 aagaagCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTcc                                               <  2:130788/100‑1 (MQ=255)
                                                 aagaagCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTcc                                               <  2:1084761/100‑1 (MQ=255)
                                                 aagaagCAGAGCTGGATTATCGCTTGAGTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTcc                                               <  1:1671278/100‑1 (MQ=255)
                                                  agaagCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCg                                              <  2:871451/100‑1 (MQ=255)
                                                    aagCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGc                                            <  2:800824/100‑1 (MQ=255)
                                                     agCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCa                                           >  1:484176/1‑100 (MQ=255)
                                                        agagCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAAtt                                        <  2:1508357/100‑1 (MQ=255)
                                                        agagCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAAtt                                        >  2:951200/1‑100 (MQ=255)
                                                         gagCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAAttg                                       >  1:1780474/1‑100 (MQ=255)
                                                         gagCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAAttg                                       <  2:744977/100‑1 (MQ=255)
                                                           gCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAAttgtt                                     >  2:1248729/1‑100 (MQ=255)
                                                           gCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAAttgtt                                     >  1:533725/1‑100 (MQ=255)
                                                           gCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAAttgtt                                     >  1:723940/1‑100 (MQ=255)
                                                           gCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAAttgtt                                     >  2:1546658/1‑100 (MQ=255)
                                                            cTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAAttgttg                                    <  2:302992/100‑1 (MQ=255)
                                                             tGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGa                                   >  2:1848663/1‑100 (MQ=255)
                                                              ggATTATCGCTTGATTAAGGTAGACCTGGGGAAAGGCGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGAt                                  >  2:751958/1‑100 (MQ=255)
                                                              ggATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCTACAACAAAATTCCGGCAATTGTTGAt                                  <  2:1183506/100‑1 (MQ=255)
                                                              ggATTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGAt                                  >  2:1923815/1‑100 (MQ=255)
                                                                aTTATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCa                                <  1:886178/100‑1 (MQ=255)
                                                                  tATCGCTTGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCAtt                              >  1:563667/1‑100 (MQ=255)
                                                                   aTCGCTTGATTAAGGTAGACCTGGTGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATtc                             <  1:1142179/100‑1 (MQ=255)
                                                                       attGATTAAGTTAGACCTGGGGAATGTCTGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCTCAACTACAAAATTCCTGCAATTGTTGATTATtc                             <  1:1142221/95‑1 (MQ=255)
                                                                         tGATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTCCTGc                       >  2:775551/1‑100 (MQ=255)
                                                                          gATTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTCCTGcc                      <  2:1584535/100‑1 (MQ=255)
                                                                           aTTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTCCTGCCg                     <  1:1015497/100‑1 (MQ=255)
                                                                           aTTAAGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTCCTGCCg                     <  1:253521/100‑1 (MQ=255)
                                                                              aaGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTCCTGCCGATg                  <  1:863807/100‑1 (MQ=255)
                                                                               aGGTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTCCTGCCGATgg                 <  1:1831702/100‑1 (MQ=255)
                                                                                ggTAGACCTGGGGAAAGGCGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTCCTGCCGATggn                <  2:194206/100‑2 (MQ=255)
                                                                                ggTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTCCTGCCGATggc                >  2:214036/1‑100 (MQ=255)
                                                                                ggTAGACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTCCTGCCGATggc                >  2:607525/1‑100 (MQ=255)
                                                                                    gACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTCCTGCCGATGgcggcg            <  1:1923802/100‑1 (MQ=255)
                                                                                    gACCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTCCTGCCGATGgcggcg            >  1:1597453/1‑100 (MQ=255)
                                                                                     aCCTGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTCCTGCCGATGgcggcg            >  1:11650/1‑99 (MQ=255)
                                                                                        tGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTCCTGCCGATGGCGGCGAAcc        >  1:1874631/1‑100 (MQ=255)
                                                                                        tGGGGAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTCCTGCCGATGGCGGCGAAcc        <  1:1184426/100‑1 (MQ=255)
                                                                                         ggggAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTCCTGCCGATGGCGGCGAACCg       >  1:1484013/1‑100 (MQ=255)
                                                                                         ggggAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTCCTGCCGATGGCGGCGAACCg       >  1:658646/1‑100 (MQ=255)
                                                                                          gggAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTCCTGCCGATGGCGGCGAACCGc      <  2:320189/100‑1 (MQ=255)
                                                                                          gggAAAGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTCCTGCCGATGGCGGCGAACCGc      >  2:869917/1‑100 (MQ=255)
                                                                                             aaaGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATCCTCCTGCCGATGGCGGCGAACCGCTa    <  2:1857220/99‑1 (MQ=255)
                                                                                              aaGACGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTCCTGCCGATGGCGGCGAACCGCTAAg  >  1:821999/1‑100 (MQ=255)
                                                                                                 |                                                                                                
CTCTATTTCGCCCCGACACCCAATGGTCACAAAATTACGCTGTTTCTCGAAGAAGCAGAGCTGGATTATCGCTTGATTAAGGTAGACCTGGGGAAAGGCGGTCAGTTTCGTCCGGAATTTTTGCGCATTTCGCCCAACAACAAAATTCCGGCAATTGTTGATCATTCTCCTGCCGATGGCGGCGAACCGCTAAG  >  NC_000913/2420630‑2420823

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