Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,629,481 G→A *64* (TAG→TAA yfgG → uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,629,4810GA97.3% 112.4 / ‑6.0 37*64* (TAG→TAAyfgGuncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (17/19);  minor base T (0/1);  total (17/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TATCGTACTGCTCATCAGCTTTATCTTCTTCTTTGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAGCGGTTACCTTCTCCACTTTCACAGAACATAACGGCACTTCGCTGTCGGATGCTTTTGCTGTTTGGGATTATCAAAGCGGCAGATATTCTTTCATCTTAAATTTTACGTCTTTATCCTGACTGATGTTTA  >  NC_000913/2629346‑2629610
                                                                                                                                       |                                                                                                                                 
tATCGTACTGCTCATCAGCTTTATCTTCTTCTTTGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAACgg                                                                                                                                <  2:6668/139‑1 (MQ=255)
        tGCTCATCAGCTTTATCTTCTTCTTTGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAACGGTTACCTtc                                                                                                                        <  1:1210649/139‑1 (MQ=255)
                 gCTTTATCTTCTTCTTTGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTc                                                                                                               >  2:1059545/1‑139 (MQ=255)
                       tcttcttctTTGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTCACAGaa                                                                                                         >  2:945487/1‑139 (MQ=255)
                                 tGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAACGGTTACCTtctc                                                                                                                      <  2:341969/116‑1 (MQ=255)
                                 tGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAACGGTTACCTtctc                                                                                                                      >  1:341969/1‑116 (MQ=255)
                                 tGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTCACAGAACATAACGGCa                                                                                               >  2:552989/1‑139 (MQ=255)
                                 tGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTCACAGAACATAACGGCa                                                                                               >  2:427126/1‑139 (MQ=255)
                                                  cgtccgtcGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTCACAGAACATAACGGCACTTCGCTGTCGGATGCt                                                                              >  2:1338474/1‑139 (MQ=255)
                                                      cgtcGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTCACAGAACATAACGGCACTTCGCTGTCGGATGCTTTTg                                                                          <  1:1084794/139‑1 (MQ=255)
                                                      cgtcGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTCACAGAACATAACGGCACTTCGCTGTCGGATGCTTTTg                                                                          >  2:1143455/1‑139 (MQ=255)
                                                              ccTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTCACAGAACATAACGGCACTTCGCTGTCGGATGCTTTTGCTGTTTgg                                                                  <  1:427126/139‑1 (MQ=255)
                                                                      caCCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTCACAGAACATAACGGCACTTCGCTGTCGGATGCTTTTGCTGTTTGGGATTATCa                                                          >  2:1149024/1‑139 (MQ=255)
                                                                           tCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTCACAGAACATAACGGCACTTCGCTGTCGGATGCTTTTGCTGTTTGGGATTATCAAAGCg                                                     >  2:547992/1‑139 (MQ=255)
                                                                            cAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTCACAGAACATAACGGCACTTCGCTGTCGGATGCTTTTGCTGTTTGGGATTATCAAAGCgg                                                    <  1:552989/139‑1 (MQ=255)
                                                                            cAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTCACAGAACATAACGGCACTTCGCTGTCGGATGCTTTTGCTGTTTGGGATTATCAAAGCgg                                                    <  1:1149024/139‑1 (MQ=255)
                                                                             agagCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTCACAGAACATAACGGCACTTCGCTGTCGGATGCTTTTGCTGTTTGGGATTATCAAAGCGGc                                                   <  1:958897/139‑1 (MQ=255)
                                                                                 cAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTca                                                                                                              >  1:1225145/1‑76 (MQ=255)
                                                                                 cAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTca                                                                                                              <  2:1225145/76‑1 (MQ=255)
                                                                                 cAAAAAAGAAGCTCAGCAATCCACACTCCCCGTCGAATCACCGGTACAACGTTAACGGGTACCATCGCCAGTTTAACAGAACATTGCGGTACTTCGCTGTCGGATGCTTTTGCTGTTCGGGATTATCCAAGCGGCAGat                                               <  1:630873/139‑1 (MQ=255)
                                                                                               aGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTCACAGAACATAACGGCACTTCGCTGTCGGATGCTTTTGCTGTTTGGGATTATCAAAGCGGCAGATATTCTTTCATCTTa                                 <  2:158440/139‑1 (MQ=255)
                                                                                               aGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTCACAGAACATAACGGCACTTCGCTGTCGGATGCTTTTGCTGTTTGGGATTATCAAAGCGGCAGATATTCTTTCATCTTa                                 >  2:104735/1‑139 (MQ=255)
                                                                                                   aTCCACACTCTCCGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTca                                                                                                              <  1:991446/58‑1 (MQ=255)
                                                                                                   aTCCACACTCTCCGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTca                                                                                                              >  2:991446/1‑58 (MQ=255)
                                                                                                      cacaCTCTCCGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTCACAGAACATAACGGCACTTCGCTGTCGGATGCTTTTGCTGTTTGGGATTATCAAAGCGGCAGATATTCTTTCATCTTAAATTTTa                          <  1:547992/139‑1 (MQ=255)
                                                                                                            ctcCGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTCACAGAACATAACGGCACTTCGCTGTCGGATGCTTTTGCTGTTTGGGATTATCAAAGCGGCAGATATTCTTTCATCTTAAATTTTACGTCtt                    >  2:994722/1‑139 (MQ=255)
                                                                                                             tcCGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTCACAGAACATAACGGCACTTCGCTGTCGGATGCTTTTg                                                                          >  2:1284929/1‑84 (MQ=255)
                                                                                                             tcCGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTCACAGAACATAACGGCACTTCGCTGTCGGATGCTTTTg                                                                          <  1:1284929/84‑1 (MQ=255)
                                                                                                               cGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTCACAGAACATAACGGCACTTCGCTGTCGGa                                                                                  >  2:1057779/1‑74 (MQ=255)
                                                                                                               cGTCGAATCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTCACAGAACATAACGGCACTTCGCTGTCGGa                                                                                  <  1:1057779/74‑1 (MQ=255)
                                                                                                                     aTCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTCACAGAACATAACGGCACTTCGCTGTCGGATGCTTTTGCTGTTTGGGATTATCAAAGCGGCAGATATTCTTTCATCTTAAATTTTACGTCTTTATCctgac           >  2:1322241/1‑139 (MQ=255)
                                                                                                                      tCACCGGTACAACGTTAACGGTTACCTTCTCCACTTTCACAGAACATAACGGCACTTCGCTGTCGGATGCTTTTGCTGTTTGGGATTATCAAAGCGGCAGATATTCTTTCATCTTAAATTTTACGTCTTTATCctgact          <  1:1338474/139‑1 (MQ=255)
                                                                                                                           ggTACAACGTTAACGGTTACCTTCTCCACTTTCACAGAACATAACGGCACTTCGCTGTCGGATGCTTTTGCTGTTTGGGATTATCAAAGCGGCAGATATTCTTTCATCTTAAATTTTACGTCTTTAt                 >  2:1319489/1‑127 (MQ=255)
                                                                                                                           ggTACAACGTTAACGGTTACCTTCTCCACTTTCACAGAACATAACGGCACTTCGCTGTCGGATGCTTTTGCTGTTTGGGATTATCAAAGCGGCAGATATTCTTTCATCTTAAATTTTACGTCTTTAt                 <  1:1319489/127‑1 (MQ=255)
                                                                                                                              aCAACGTTATCGGTTACATTCTCCGCCTTCACAGAACATTACGGCTCTTCGCTGTCGGATGGTTTGGCTGTTTGGGATTATCAAAGCGGCAGATATTCTTTCATCTTAAATTTTACGTCTTTATCCTGACTGATGTTTa  <  1:994722/139‑1 (MQ=255)
                                                                                                                              aCAACGTTAACGGTTACCTTCTCCACTTTCACAGAACATAACGGCACTTCGCTGTCGGATGCTTTTGCTGTTTGGGATTATCAAAGCGGCAGATATTCTTTCATCTTAAATTTTACGTCTTTATCCTGACTGATGTTTa  <  1:1322241/139‑1 (MQ=255)
                                                                                                                              aCAACGTTAACGGTTACCTTCTCCACTTTCACAGAACATAACGGCACTTCGCTGTCGGATGCTTTTGCTGTTTGGGATTATCAAAGCGGCAGATATTCTTTCATCTTAAATTTTACGTCTTTATCCTGACTGATGTTTa  <  1:1143455/139‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                 
TATCGTACTGCTCATCAGCTTTATCTTCTTCTTTGGCCGTTTTATCTACTCGTCCGTCGGTGCCTGGCAGCACCATCAGAGCAAAAAAGAAGCTCAGCAATCCACACTCTCCGTCGAATCACCGGTACAACGTTAGCGGTTACCTTCTCCACTTTCACAGAACATAACGGCACTTCGCTGTCGGATGCTTTTGCTGTTTGGGATTATCAAAGCGGCAGATATTCTTTCATCTTAAATTTTACGTCTTTATCCTGACTGATGTTTA  >  NC_000913/2629346‑2629610

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: