Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 297,827 T→C E300E (GAA→GAG)  paoD ← molybdenum cofactor insertion chaperone for PaoABC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913297,8270TC96.5% 85.7 / NA 29E300E (GAA→GAG) paoDmolybdenum cofactor insertion chaperone for PaoABC
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (11/17);  minor base G (1/0);  total (12/17)
Fisher's exact test for biased strand distribution p-value = 4.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGACTGACGCCAGCCGATACACTTATGGGTATTTCCCCCGGAGGCGAGAAAGCACTCTCCACGCCCGGCCGCAAGGATCAGGACGACGGGGGCAGGCATGAATCCTCCTCCTGATGGAGACGTACAGAGGCGACTTCTGCCAGCACGGAGAGTGCCAGAGTATGCGCATCCCGGGCTTTGGGGAATATCCCGACGGGTGCCCGGATTTGCGTTGTTTCCTCCCTGGACCATCCCAGCTCGTGGAGCTTTTGCAGACGTAACGTGTGGGTT  >  NC_000913/297692‑297962
                                                                                                                                       |                                                                                                                                       
gggACTGACGCCAGCCGATACACTTATGGGTATTTCCCCCGGAGGCGAGAAAGCACTCTCCACGCCCGGCCGCAAGGATCAGGACGACGGGGGCAGGCATGAATCCTCCTCCTGATGGAGACGTACAGAGGCGACCTCTg                                                                                                                                     <  1:501281/140‑1 (MQ=255)
  gACTGACGCCAGCCGATACACTTATGGGTATTTCCCCCGGAGGCGAGAAAGCACTCTCCACGCCCGGCCGCAAGGATCAGGACGACGGGGGCAGGCATGAATCCTCCTCCTGATGGAGACGTACAGAGGCGACCTCTGcc                                                                                                                                   <  2:846424/140‑1 (MQ=255)
          ccAGCCGATACACTTATGGGTATTTCCCCCGGAGGCGAGAAAGCACTCTCCACGCCCGGCCGCAAGGATCAGGACGACGGGGGCAGGCATGAATCCTCCTCCTGATGGAGACGTACAGAGGCGACCTCTGCCAGCACGga                                                                                                                           <  1:658045/140‑1 (MQ=255)
             gCCGATACACTTATGGGTATTTCCCCCGGAGGCGAGAAAGCACTCTCCACGCCCGGCCGCAAGGATCAGGACGACGGGGGCAGGCATGAATCCTCCTCCTTATGGAGAGGTACAGAGACGACGTCTGCCAGCACGGagag                                                                                                                        >  2:269936/1‑140 (MQ=255)
                gATACACTTATGGGTCTTTCCCCCGGAGGCGAGAAAGAACTCTCCACGCCCGGCCGCAAGGATCAGGACGACGGGGGCAGGCATGAATCCTCCTCCTGATGGAGACGTACAGAGGCGACCTCTGCCAGCACGGAGAGTGc                                                                                                                     <  1:269936/140‑1 (MQ=255)
                 aTACACTTATGGGTATTTCCCCCGGAGGCGAGAAAGCACTCTCCACGCCCGGCCGCAAGGATCAGGACGACGGGGGCAGGCATGAATCCTCCTCCTGATGGATACGTACAAAGGGGACCTCTGCCAGCACGGAGAGTGcc                                                                                                                    >  1:281506/1‑140 (MQ=255)
                 aTACACTTATGGGTATTTCCCCCGGAGGCGAGAAAGCACTCTCCACGCCCGGCCGCAAGGATCAGGACGACGGGGGCAGGCATGAATCCTCCTCCTGATGGAGACGTACAGAGGCGACCTCTGCCAGCACGGAGAGTGcc                                                                                                                    >  1:629747/1‑140 (MQ=255)
                            ggTATTTCCCCCGGAGGCGAGAAAGCACTCTCCACGCCCGGCCGCAAGGATCAGGACGACGGGGGCAGGCATGAATCCTCCTCCTGATGGAGACGTACAGAGGCGACCTCTGCCAGCACGGAGAGTGCCAGAGTATgcgc                                                                                                         >  1:550495/1‑140 (MQ=255)
                                    ccccGGAGGCGAGAAAGCACTCTCCACGCCCGGCCGCAAGGATCAGGACGACGGGGGCAGGCATGAATCCTCCTCCTGATGGAGACGTACAGAGGCGACCTCTGCCAGCACGGAGAGTGCCAGAGTATGCGCATCCCggg                                                                                                 <  2:508223/140‑1 (MQ=255)
                                              gagaAAGCACTCTCCACGCCCGGCCGCAAGGATCAGGACGACGGGGGCAGGCATGAATCCTCCTCCTGATGGAGACGTACAGAGGCGACCTCTGCCAGCACGGAGAGTGCCAGAGTATGCGCATCCCGGGCTTTGGGGaa                                                                                       <  1:507934/140‑1 (MQ=255)
                                                              cGCCCGGCCGCAAGGCTCAGGACGATGGGGGCAGGCATGAATCCTCCTCCTGATGGAGACGTACAGAGGCGACCTCTGCCAGCACGGAGAGTGCCAGAGTATGCGCATCCCGGGCTTTGGGGAATATCCCGACGGGTGcc                                                                       <  1:367104/140‑1 (MQ=255)
                                                                  cGGCCGCAAGGATCAGGACGACGGGGGCAGGCATGAATCCTCCTCCTGATGGAGACGTACAGAGGCGACCTCTGCCAGCACGGAGAGTGCCAGAGTATGCGCATCCCGGGCTTTGGGGAATATCCCGACGGGTGCCCGGa                                                                   >  2:830473/1‑140 (MQ=255)
                                                                    gCCGCAAGGATCAGGACGACGGGGGCAGGCATGAATCCTCCTCCTGATGGAGACGTACAGAGGCGACCTCTGCCAGCACGGAGAGTGCCAGAGTATGCGCATCCCGGGCTTTGGGGAATATCCCGACGGGTGCCCGGAtt                                                                 >  1:176440/1‑140 (MQ=255)
                                                                         aaGGATCAGGACGACGGGGGCAGGCATGAATCCTCCTCCTGATGGAGACGTACAGAGGCGACCTCTGCCAGCACGGAGAGTGCCAGAGTATGCGCATCCCGGGCTTTGGGGAATATCCCGACGGGTGCCCGGATTTGCgt                                                            <  2:19909/140‑1 (MQ=255)
                                                                         aaGGATCAGGACGACGGGGGCAAGCATGAATCCTCTTCCTGATGGAGACGTACAGAGGCGACCTCTGCCAGCACGGAGAGTGCCAGAGTATGCGCATCCCGGGCTTTGGGGAATATCCCGACGGGTGCCCGGATTTGCgt                                                            <  2:344573/140‑1 (MQ=255)
                                                                              tCAGGACGACGGGGGCAGGCATGAATCCTCCTCCTGATGGAGACGTACAGAGGCGACCTCTGCCAGCACGGAGAGTGCCAGAGTATGCGCATCCCGGGCTTTGGGGAATATCCCGACGGGTGCCCGGATTTGCGTTGttt                                                       <  2:664216/140‑1 (MQ=255)
                                                                                       cgGGGGCAGGCATGAATCCTCCTCCTGATGGAGACGTACAGAGGCGACCTCTGCCAGCACGGAGAGTGCCAGAGTATGCGCATCCCGGGCTTTGGGGAATATCCCGCCGGGTGCCCGGATTTGCGTTGTTTCCTCCCTgg                                              >  2:163618/1‑140 (MQ=255)
                                                                                       cgGGGGCAGGCATGAATCCTCCTCCTGATGGAGACGTACAGAGGCGACCTCTGCCAGCACGGAGAGTGCCAGAGTATGCGCATCCCGGGCTTTGGGGAATATCCCGACGGGTGCCCGGATTTGCGTTGTTTCCTCCCTgg                                              >  2:663776/1‑140 (MQ=255)
                                                                                       cgGGGGCAGGCATGAATCCTCCTCCTGATGGAGACGTACAGAGGCGACCTCTGCCAGCACGGAGAGTGCCAGAGTATGCGCATCCCGGGCTTTGGGGAATATCCCGACGGGTGCCCGGATTTGCGTTGTTTCCTCCCTgg                                              >  2:46110/1‑140 (MQ=255)
                                                                                        gggggCAGGCATGAATCCTCCTCCTGATGGAGACGTACAGAGGCGACCTCTGCCAGCACGGAGAGTGCCAGAGTATGCGCATCCCGGGCTTTGGGGAATATCCCGACGGGTGCCCGGATTTGCGTTGTTTCCTCCCTGGa                                             <  1:663776/140‑1 (MQ=255)
                                                                                                         ctcctcctGATGGAGACTTACAAAGGCGACCTCTGGCAGTACGGAGAGCGCCAGAGTATGCGCATCCCGGGCTTTGGGGAATATCCCGACGGGTGCCCGGATTTGCGTTGTTTCCTCCCTGGACCATCCCAGCTCGTGGa                            <  2:281506/140‑1 (MQ=255)
                                                                                                         ctcctcctGATGGAGACGTACAGAGGCGACCTCTGCCAGCACGGAGAGTGCCAGAGTATGCGCATCCCGGGCTTTGGGGAATATCCCGACGGGTGCCCGGATTTGCGTTGTTTCCTCCCTGGACCATCCCAGCTCGTGGa                            <  2:629747/140‑1 (MQ=255)
                                                                                                             tcctGATGGAGACGTACAGAGGCGACCTCTGCCAGCACGGAGAGTGCCAGAGTATGCGCATCCCGGGCTTTGGGGAATATCCCGACGGGTGCCCGGATTTGCGTTGTTTCCTCCCTGGACCATCCCAGCTCGTGGAGCtt                        <  2:550495/140‑1 (MQ=255)
                                                                                                                  gagacagcGTACAGACGCGACCCCTCCCAACACGGAGCGTGCCAGAGTATGCGCATCCCGGGCTTTGGGGAATATCCCGACGGGTGCCCGGATTTGCGTTGTTTCCTCCCTGGACCATCCCAGCTCGTGGAGCTTTTGCa                   <  1:242790/133‑1 (MQ=255)
                                                                                                                  aTGGAGACGTACAGAGGAGACCGCTGCCAGCGCGGAGAGGGCCAGAGTGTGCGCATCCCGGGCTTTGGGGAATATCCCGACGGGTGCCCGGATTTGCGTTGTTTCCTCCCTGGGCCATCCCAGCTCGTGGAGCTTTTGCa                   <  2:298239/140‑1 (MQ=255)
                                                                                                                    ggAGACGTACAGAGGCGACCTCTGCCAGCACGGAGAGTGCCAGAGTATGCGCATCCCGGGCTTTGGGGAATATCCCGACGGGTGCCCGGATTTGCGTTGTTTCCTCCCTGGACCATCCCAGCTCGTGGAGCTTTTg                     <  2:742357/136‑1 (MQ=255)
                                                                                                                    ggAGACGTACAGAGGCGACCTCTGCCAGCACGGAGAGTGCCAGAGTATGCGCATCCCGGGCTTTGGGGAATATCCCGACGGGTGCCCGGATTTGCGTTGTTTCCTCCCTGGACCATCCCAGCTCGTGGAGCTTTTg                     >  1:742357/1‑136 (MQ=255)
                                                                                                                     gagaCGTACAGAGGCGACCTCTGCCAGCACGGAGAGTGCCAGAGTATGCGCATCCCGGGCTTTGGGGAATATCCCGACGGGTGCCCGGATTTGCGTTGTTTCCTCCCTGGACCATCCCAGCTCGTGGAGCTTTTGCAGAc                >  1:690902/1‑140 (MQ=255)
                                                                                                                         cGTACAGAGGCGACCTCTGCCAGCACGGAGAGTGCCAGAGTATGCGCATCCCGGGCTTT‑GGGAATATCCCGACGGGTGCCCGGATTTGCGTTGTTTCCTCCCTGGACCATCCCAGCTCGTGGAGCTTTTGCAGACGTAAc           >  1:757454/1‑140 (MQ=255)
                                                                                                                            aCAGAGGCGACCTCTGCCAGCACGGAGAGTGCCAGAGTATGCGCATCCCGGGCTTTGGGGAATATCCCGACGGGTGCCCGGATTTGCGTTGTTTCCTCCCTGGACCATCCCAGCTCGTGGAGCTTTTGCAGACGTAACgt         <  2:176440/140‑1 (MQ=255)
                                                                                                                                   cGACCTCTGCCAGCACGGAGAGTGCCAGAGTATGGGCATCCCGGGCTTTGGGGAATATCCCGACGGGTGCCCGGATTTGCGTTGTTTCCTCCCTGGACCATCCCAGCTCGTGGAGCTTTTGCAGACGTAACGTGTGGGtt  <  2:649586/140‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                       
GGGACTGACGCCAGCCGATACACTTATGGGTATTTCCCCCGGAGGCGAGAAAGCACTCTCCACGCCCGGCCGCAAGGATCAGGACGACGGGGGCAGGCATGAATCCTCCTCCTGATGGAGACGTACAGAGGCGACTTCTGCCAGCACGGAGAGTGCCAGAGTATGCGCATCCCGGGCTTTGGGGAATATCCCGACGGGTGCCCGGATTTGCGTTGTTTCCTCCCTGGACCATCCCAGCTCGTGGAGCTTTTGCAGACGTAACGTGTGGGTT  >  NC_000913/297692‑297962

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