Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,054,234 A→G E232G (GAA→GGA)  torC → pentaheme c‑type cytochrome TorC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,054,2340AG98.6% 253.2 / ‑3.9 69E232G (GAA→GGA) torCpentaheme c‑type cytochrome TorC
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (47/21);  total (48/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.16e-01

CACGCTGTACTCTATTGATATTAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGAAGTGAAAGTCCTTAAACGTGACGGCGACTGGCTGCAAATTGAAATTACCGGCTGGACGGAAAGCGCCGGACGTCAGCGTGTACTCA  >  NZ_CP009273/1054151‑1054320
                                                                                   |                                                                                      
cacGCTGTACTCTATTGATATTAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGGAGTGaa                                                                                  <  2:2374408/90‑1 (MQ=255)
     tGTACTCTATTGATATTAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGGAGTGAAAGTcc                                                                             >  2:1128494/1‑90 (MQ=255)
        aCTCTATTGATATTAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGGAGTGAAAGTCCTTa                                                                          <  1:2579108/90‑1 (MQ=255)
          tctATTGATATTAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGGAGTGAAAGTCCTTaaa                                                                        >  2:401928/1‑90 (MQ=255)
          tctATTGATATTAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGGAGTGAAAGTCCTTaaa                                                                        >  1:1216575/1‑90 (MQ=255)
              ttGATATTAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGGAGTGAAAGTCCTTAAACGTg                                                                    >  2:1142881/1‑90 (MQ=255)
              ttGATATTAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGGAGTGAAAGTCCTTAAACGTg                                                                    >  2:451673/1‑90 (MQ=255)
              ttGATATTAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGGAGTGAAAGTCCTTAAACGTg                                                                    >  1:2142730/1‑90 (MQ=255)
                gATATTAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGGAGTGAAAg                                                                                >  1:3022359/1‑76 (MQ=255)
                gATATTAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGGAGTGAAAg                                                                                <  2:3022359/76‑1 (MQ=255)
                gATATTAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGGAGTGAAAGTCCTTAAACGTGAc                                                                  >  2:726767/1‑90 (MQ=255)
                gATATTAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGGAGTGAAAGTCCTTAAACGTGAc                                                                  >  1:577502/1‑90 (MQ=255)
                gATATTAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGGAGTGAAAGTCCTTAAACGTGAc                                                                  >  1:1136333/1‑90 (MQ=255)
                gATATTAAGCCAATTTATGCGGCGAAAGGCGATAAAAAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGGGGGGAAAGTCCTTAAACGTgcc                                                                  >  1:1029186/1‑88 (MQ=255)
                  tatTAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGGAGTGAAAGTCCTTAAACGTGACgg                                                                >  2:227049/1‑90 (MQ=255)
                     tAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGGAGTGAAAGTCCTTAAACGTGACGGCGa                                                             >  1:217526/1‑90 (MQ=255)
                     tAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGGAGTGAAAGTCCTTAAACGTGACGGCGa                                                             >  1:415512/1‑90 (MQ=255)
                     tAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGGAGTGAAAGTCCTTAAACGTGACGGCGa                                                             >  1:2964187/1‑90 (MQ=255)
                     tAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGGAGTGAAAGTCCTTAAACGTGACGGCGa                                                             >  2:1074075/1‑90 (MQ=255)
                     tAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGGAGTGAAAGTCCTTAAACGTGACGGCGa                                                             >  1:1986855/1‑90 (MQ=255)
                     tAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGGAGTGAAAGTCCTTAAACGTGACGGCGa                                                             >  2:2134180/1‑90 (MQ=255)
                     tAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGGAGTGAAAGTCCTTAAACGTGACGGCGa                                                             >  2:1106805/1‑90 (MQ=255)
                     tAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGGAGTGAAAGTCCTTAAACGTGACGGCGa                                                             >  2:842536/1‑90 (MQ=255)
                     tAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGGAGTGAAAGTCCTTAAACGTGACGGCGa                                                             >  2:923022/1‑90 (MQ=255)
                     tAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCACTGCTGCCTGCTTCGGGAGTGAAAGTCCTTAAACGTGACGGCGa                                                             >  2:423569/1‑90 (MQ=255)
                       aGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGGAGTGAAAGTCCTTAAACGTGACGGCGACt                                                           >  2:1068955/1‑90 (MQ=255)
                       aGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGGAGTGAAAGTCCTTAAACGTGACGGCGACt                                                           >  1:1095592/1‑90 (MQ=255)
                       aGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGGAGTGAAAGTCCTTAAACGTGACGGCGACt                                                           >  2:2277050/1‑90 (MQ=255)
                       aGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGGAGTGAAAGTCCTTAAACGTGACGGCGACt                                                           >  2:1909588/1‑90 (MQ=255)
                       aGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGGAGTGAAAGTCCTTAAACGTGACGGCGACt                                                           >  2:2706932/1‑90 (MQ=255)
                          cAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGGAGTGAAAGTCCTTAAACGTGACGGCGACTGGc                                                        <  2:882564/90‑1 (MQ=255)
                          cAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGGAGTGAAAGTCCTTAAACGTGACGGCGACTGGc                                                        >  1:1300240/1‑90 (MQ=255)
                           aaTTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGGAGTGAAAGTCCTTAAACGTGACGGCGACTGGCt                                                       >  1:2540885/1‑90 (MQ=255)
                            aTTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGGAGTGAAAGTCCTTAAACGTGACGGCGACTGGCTg                                                      >  1:2488180/1‑90 (MQ=255)
                                          aGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGGAGTGAAAGTCCTTAAACGTGACGGCGACTGGCTGCAAATTGAAATTAc                                        <  1:451673/90‑1 (MQ=255)
                                              gATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGGAGTGAAAGTCCTTAAACGTGACGGCGACTGGCTGCAAATTGAAATTACCGGc                                    >  2:245492/1‑90 (MQ=255)
                                               aTAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGGAGTGAAAGTCCTTAAACGTGACGGCGACTGGCTGCAAATTGAAATTACCGGCt                                   <  1:2277050/90‑1 (MQ=255)
                                                    gaagCCTCTGGTTCTCTGCTGCCTGCTTCGGGAGTGAAAGTCCTTAAACGTGACGGCGACTGGCTGCAAATTGAAATTACCGGCTGGACg                              <  1:2706932/90‑1 (MQ=255)
                                                    gaagCCTCTGGTTCTCTGCTGCCTGCTTCGGGAGTGAAAGTCCTTAAACGTGACGGCGACTGGCTGCAAATTGAAATTACCGGCTGGACg                              >  2:1278191/1‑90 (MQ=255)
                                                     aagCCTCTGGTTCTCTGCTGCCTGCTTCGGGAGTGAAAGTCCTTAAACGTGACGGCGACTGGCTGCAAATTGAAATTACCGGCTGGACgg                             >  1:2810315/1‑90 (MQ=255)
                                                      agCCTCTGGTTCTCTGCTGCCTGCTTCGGGAGTGAAAGTCCTTAAACGTGACGGCGACTGGCTGCAAATTGAAATTACCGGCTGGACGGa                            >  1:1911321/1‑90 (MQ=255)
                                                      agCCTCTGGTTCTCTGCTGCCTGCTTCGGGAGTGAAAGTCCTTAAACGTGACGGCGACTGGCTGCAAATTGAAATTACCGGCTGGACGGa                            >  2:791538/1‑90 (MQ=255)
                                                         ctctGGTTCTCTGCTGCCTGCTTCGGGAGTGAAAGTCCTTAAACGTGACGGCGACTGGCTGCAAATTGAAATTACCGGCTGGACGGAAAg                         >  2:1321539/1‑90 (MQ=255)
                                                           ctGGTTCTCTGCTGCCTGCTTCGGGAGTGAAAGTCCTTAGACGTGACGGCGACTGGCTGCAAATTGAAATTACCGGCTGGACGGAAAgcg                       >  2:2710153/1‑90 (MQ=255)
                                                           ctGGTTCTCTGCTGCCTGCTTCGGGAGTGAAAGTCCTTAAACGTGACGGCGACTGGCTGCAAATTGAAATTACCGGCTGGACGGAAAgcg                       >  1:716951/1‑90 (MQ=255)
                                                            tGGTTCTCTGCTGCCTGCTTCGGGAGTGAAAGTCCTTAAACGTGACGGCGACTGGCTGCAAATTGAAATTACCGGCTGGACGGAAAgcgc                      <  1:227049/90‑1 (MQ=255)
                                                             ggTTCTCTGCTGCCTGCTTCGGGAGTGAAAGTCCTTAAACGTGACGGCGACTGGCTGCAAATTGAAATTACCGGCTGGACGGAAAGCGcc                     <  1:791538/90‑1 (MQ=255)
                                                                 ctctGCTGCCTGCTTCGGGAGTGAAAGTCCTTAAACGTGACGGCGACTGGCTGCAAATTGAAATTACCGGCTGGACGGAAAGCGCCGGAc                 >  2:1777719/1‑90 (MQ=255)
                                                                  tctGCTGCCTGCTTCGGGAGTGAAAGTCCTTAAACGTGACGGCGACTGGCTGCAAATTGAAATTACCGGCTGGACGGAAAGCGCCGGACg                >  2:2242314/1‑90 (MQ=255)
                                                                  tctGCTGCCTGCTTCGGAAGTGAAAGTCCTTAAACGTGACGGCGACTGGCTGCAAATTGAAATTACCGGCTGGACGGAAAGCGCCGGACg                >  2:1317829/1‑90 (MQ=255)
                                                                   ctgctgCCTGCTTCGGGAGTGAAAGTCCTTAAACGTGACGGCGACTGGCTGCAAATTGAAATTACCGGCTGGACGGAAAgcgc                      >  2:1783287/1‑83 (MQ=255)
                                                                   ctgctgCCTGCTTCGGGAGTGAAAGTCCTTAAACGTGACGGCGACTGGCTGCAAATTGAAATTACCGGCTGGACGGAAAgcgc                      <  1:1783287/83‑1 (MQ=255)
                                                                   ctgctgCCTGCTTCGGGAGTGAAAGTCCTTAAACGTGACGGCGACTGGCTGCAAATTGAAATTACCGGCTGGACGGAAAGCGCCGGACGt               <  2:1864151/90‑1 (MQ=255)
                                                                   ctgctgCCTGCTTCGGGAGTGAAAGTCCTTAAACGTGACGGCGACTGGCTGCAAATTGAAATTACCGGCTGGACGGAAAGCGCCGGACGt               <  1:2122141/90‑1 (MQ=255)
                                                                   ctgctgCCTGCTTCGGGAGTGAAAGTCCTTAAACGTGACGGCGACTGGCTGCAAATTGAAATTACCGGCTGGACGGAAAGCGCCGGACGt               >  1:38465/1‑90 (MQ=255)
                                                                    tgctgcCTGCTTCGGGAGTGAAAGTCCTTAAACGTGACGGCGACTGGCTGCAAATTGAAATTACCGGCTGGACGGAAAGCGCCGGACGTc              <  1:1909588/90‑1 (MQ=255)
                                                                    tgctgcCTGCTTCGGGAGTGAAAGTCCTTAAACGTGACGGCGACTGGCTGCAAATTGAAATTACCGGCTGGACGGAAAGCGCCGGACGTc              <  1:1142881/90‑1 (MQ=255)
                                                                    tgctgcCTGCTTCGGGAGTGAAAGTCCTTAAACGTGACGGCGACTGGCTGCAAATTGAAATTACCGGCTGGACGGAAAGCGCCGGACGTc              <  1:2134180/90‑1 (MQ=255)
                                                                     gctgcCTGCTTCGGGAGGGAAAGTCCTTAAACGTGACGGCGACTGGCTGCAAATTGAAATTACCGGCTGGACGGAAAGCGCCGGACGTCa             >  1:1649393/1‑90 (MQ=255)
                                                                       tgcctgcTTCGGGAGTGAAAGTCCTTAAACGTGACGGCGACTGGCTGCAAATTGAAATTACCGGCTGGACGGAAAGCGCCGGACGTCAGc           >  2:1351268/1‑90 (MQ=255)
                                                                       tgcctgcTTCGGGAGTGAAAGTCCTTAAACGTGACGGCGACTGGCTACAAATTGAAATTACCGGCTGGGCGGAAAGCGCCGGACGGCAGc           >  1:1811224/1‑90 (MQ=255)
                                                                        gcctgcTTCGGGAGTGAAAGTCCTTAAACGTGACGTCGACTGGCTGCAAATTGAAATTACCGGCTGGACGGAAAGCGCCGGACGTCAGCg          >  2:96170/1‑90 (MQ=255)
                                                                        gcctgcTTCGGGAGTGAAAGTCCTTAAACGTGACGGCGACTGGCTGCAAATTGAAATTACCGGCTGGACGGAAAGCGCCGGACGTCAGCg          >  1:3010282/1‑90 (MQ=255)
                                                                        gcctgcTTCGGGAGTGAAAGTCCTTAAACGTGACGGCGACTGGCTGCAAATTGAAATTACCGGCTGGACGGAAAGCGCCGGACGTCAGCg          <  1:245492/90‑1 (MQ=255)
                                                                        gcctgcTTCGGGAGTGAAAGTCCTTAAACGTGACGGCGACTGGCTGCAAATTGAAATTACCGGCTGGACGGAAAGCGCCGGACGTCAGCg          <  1:423569/90‑1 (MQ=255)
                                                                              ttCGGGAGTGAAAGTCCTTAAACGTGACGGCGACTGGCTGCAAATTGAAATTACCGGCTGGACGGAAAGCGCCGGACGTCAGCGTGTACt    <  2:1095592/90‑1 (MQ=255)
                                                                               tCGGGAGTGAAAGTCCTTAAACGTGACGGCGACTGGCTGCAAATTGAAATTACCGGCTGGACGGAAAGCGTCGGACGTCAGCGTGTACTc   <  1:1977960/90‑1 (MQ=255)
                                                                               tCGGGAGTGAAAGTCCTTAAACGTGACGGCGACTGGCTGCAAATTGAAATTACCGGCTGGACGGAAAGCGCCGGACGTCAGCGTGTACTc   <  1:401928/90‑1 (MQ=255)
                                                                                cGGGAGTGAAAGTCCTTAAACGTGACGGCGACTGGCTGCAAATTGAAATTACCGGCTGGACGGAAAGCGCCGGACGTCAGCGTGTACTCa  <  2:1771883/90‑1 (MQ=255)
                                                                                   |                                                                                      
CACGCTGTACTCTATTGATATTAAGCCAATTTATGCGGCGAAAGGCGATAAAGAAGCCTCTGGTTCTCTGCTGCCTGCTTCGGAAGTGAAAGTCCTTAAACGTGACGGCGACTGGCTGCAAATTGAAATTACCGGCTGGACGGAAAGCGCCGGACGTCAGCGTGTACTCA  >  NZ_CP009273/1054151‑1054320

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