Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,893,928 G→A *424* (TAG→TAA
A4T (GCC→ACC) 
ygcN →
ygcO →
putative oxidoreductase
putative 4Fe‑4S cluster‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,893,9280GA96.4% 162.1 / ‑1.8 56*424* (TAG→TAA
A4T (GCC→ACC) 
ygcN
ygcO
putative oxidoreductase
putative 4Fe‑4S cluster‑containing protein
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (31/23);  total (32/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.51e-01

GGATATTTCCCGCGATTTATGGGATCAGGGTGATAAACCTGTTCCACCGCTGCGCCAGTTATTCTGGCATCATTTACGTCGTCACGGCCTGTGGCATCTGGCGGGCGATGTTATCAGGAGTCTGCGATGTCTGTAGCCCGTAATCTCTGGCGCGTTGCTGATGCGCCGCACATTGTTCCGGCTGACTCCGTTGAGCGCCAGACGGCAGAACGGTTGATTAACGCCTGTCCGGCAGGTCTTTTTTCGCTCACACCGGAAGGTAA  >  NC_000913/2893793‑2894055
                                                                                                                                       |                                                                                                                               
ggATATTTCCCGCGATTTATGGGATCAGGGTGATAAACCTGTTCCACCGCTGCGCCAGTTATTCTGGCATCATTTACGTCGTCACGGCCTGTGGCATCTGGCGGGCGATGTTATCAGGAGTCTGCGATGTCTGTAAccc                                                                                                                              >  1:197855/1‑139 (MQ=255)
ggATATTTCCCGCGATTTATGGGATCAGGGTGATAAACCTGTTCCACCGCTGCGCCAGTTATTCTGGCATCATTTACGTCGTCACGGCCTGTGGCATCTGGCGGGCGATGTTATCAGGAGTCTGCGATGTCTGTAAccc                                                                                                                              >  1:336860/1‑139 (MQ=255)
 gATATTTCCCGCGATTTATGGGATCAGGGTGATAAACCTGTTCCACCGCTGCGCCAGTTATTCTGGCATCATTTACGTCGTCACGGCCTGTGGCATCTGGCGGGCGATGTTATCAGGAGTCTGCGATGTCTGTAACCCg                                                                                                                             >  1:1213556/1‑139 (MQ=255)
        cccGCGATTTATGGGATCAGGGTGATAAACCTGTTCCACCGCTGCGCCAGTTATTCTGTCATCATTTCCGTCGTCACGGCCTGTGTCTTCTGCCGGGCGATGTTATCAGGAGTCTGCGATGTCTGTAACCCGTAAtctc                                                                                                                      >  2:1023379/1‑139 (MQ=255)
           gcgATTTATGGGATCAGGGTGATAAACCTGTTCCACCGCTGCGCCAGTTATTCTGGCATCATTTACGTCGTCACGGCCTGTGGCATCTGGCGGGCGATGTTATCAGGAGTCTGCGATGTCTGTAACCCGTAATCTCTgg                                                                                                                   >  2:171101/1‑139 (MQ=255)
             gATTTATGGGATCAGGGTGATAAACCTGTTCCACCGCTGCGCCAGTTATTCTGGCATCATTTACGTCGTCACGGCCTGTGGCATCTGGCGGGCGATGTTATCAGGAGTCTGCGATGTCTGTAACCCGTAATCTCTGgcg                                                                                                                 <  2:197855/139‑1 (MQ=255)
                           gggTGATAAACCTGTTCCACCGCTGCGCCAGTTATTCTGGCATCATTTACGTCGTCACGGCCTGTGGCATCTGGCGGGCGATGTTATCAGGAGTCTGCGATGTCTGTAACCCGTAATCTCTGGCGCGt                                                                                                              >  1:1054685/1‑128 (MQ=255)
                           gggTGATAAACCTGTTCCACCGCTGCGCCAGTTATTCTGGCATCATTTACGTCGTCACGGCCTGTGGCATCTGGCGGGCGATGTTATCAGGAGTCTGCGATGTCTGTAACCCGTAATCTCTGGCGCGt                                                                                                              <  2:1054685/128‑1 (MQ=255)
                               gATAAACCTGTTCCACCGCTGCGCCAGTTATTCTGGCATCATTTACGTCGTCACGGCCTGTGGCATCTGGCGGGCGATGTTATCAGGAGTCTGCGATGTCTGTAGCCCGTAATCTCTGGCGCGTTGCTGATGCGCCGca                                                                                               >  1:729946/1‑139 (MQ=255)
                                          tCCACCGCTGCGCCAGTTATTCTGGCATCATTTACGTCGTCACGGCCTGTGGCATCTGGCGGGCGATGTTATCAGGAGTCTGCGATGTCTGTAACCCGTAATCTCTGGCGCGTTGCTGATGCGCCGattgttccggctg                                                                                    <  2:386298‑M1/139‑14 (MQ=255)
                                          tCCACCGCTGCGCCAGTTATTCTGGCATCATTTACGTCGTCACGGCCTGTGGCATCTGGCGGGCGATGTTATCAGGAGTCTGCGATGTCTGTAACCCGTAATCTCTGGCGCGTTGCTGATGCGCCGattgttccggc                                                                                      <  2:450543‑M1/137‑12 (MQ=255)
                                          tCCACCGCTGCGCCAGTTATTCTGGCATCATTTACGTCGTCACGGCCTGTGGCATCTGGCGGGCGATGTTATCAGGAGTCTGCGATGTCTGTAACCCGTAATCTCTGGCGCGTTGCTGATGCGCCGattgttccggc                                                                                      >  1:450543‑M1/1‑126 (MQ=255)
                                                    cgcCAGTTATTCTGGCATCATTTACGTCGTCACGGCCTGTGGCATCTGGCGGGCGATGTTATCAGGAGTCTGCGATGTCTGTAACCCGTAATCTCTGGCGCGTTGCTGATGCGCCGattgttccggctgactccgttga                                                                          >  1:964816‑M1/1‑116 (MQ=255)
                                                            aTTCTGGCATCATTTACGTCGTCACGGCCTGTGGCATCTGGCGGGCGATGTTATCAGGAGTCTGCGATGTCTGTAACCCGTAATCTCTGGCGCGTTGCTGATGCGCCGattgttccggctgactccgttgagcgccaga                                                                  >  2:1002141‑M1/1‑108 (MQ=255)
                                                                        tttACGTCGTCACGGCCTGTGGCATCTGGCGGGCGATGTTATCAGGAGTCTGCGATGTCTGTAACCCGTAATCTCTGGCGCGTTGCTGATGCGCCGattgttccggctgactccgttgag                                                                         <  1:221879‑M1/120‑25 (MQ=255)
                                                                        tttACGTCGTCACGGCCTGTGGCATCTGGCGGGCGATGTTATCAGGAGTCTGCGATGTCTGTAACCCGTAATCTCTGGCGCGTTGCTGATGCGCCGattgttccggctgactccgttgag                                                                         >  2:221879‑M1/1‑96 (MQ=255)
                                                                            cgtcgtCACGGCCTGTGGCATCTGGCGGGCGATGTTATCAGGAGTCTGCGATGTCTGTAAccc                                                                                                                              <  1:25328/63‑1 (MQ=255)
                                                                            cgtcgtCACGGCCTGTGGCATCTGGCGGGCGATGTTATCAGGAGTCTGCGATGTCTGTAAccc                                                                                                                              >  2:25328/1‑63 (MQ=255)
                                                                                gtcACGGCCTGTGGCATCTGGCGGGCGATGTTATCAGGAGTCTGCGATGTCTGTAACCCGTAATCTCTGGCGCGTTGCTGATGCGCCGattgttccggctgactccgttgagcgccagacggcagaacggttgattaac                                              >  1:346695‑M1/1‑88 (MQ=255)
                                                                                gtcACGGCCTGTGGCATCTGGCGGGCGATGTTATCAGGAGTCTGCGATGTCTGTAACCCGTAATCTCTGGCGCGTTGCTGATGCGCCGattgttccggctgactccgttgagcgccagacggcagaacggttgattaac                                              >  1:485914‑M1/1‑88 (MQ=255)
                                                                                gtcACGGCCTGTGGCATCTGGCGGGCGATGTTATCAGGAGTCTGCGATGTCTGTAACCCGTAATCTCTGGCGCGTTGCTGATGCGCCGattgttccggctgactccgttgagcgccagacggcagaacggttgattaac                                              >  1:1210304‑M1/1‑88 (MQ=255)
                                                                                       ccTGTGGCATCTGGCGGGCGATGTTATCAGGAGTCTGCGATGTCTGTAACCCGTAATCTCTGGCGCGTTGCTGATGCGCCGattgttccggctgactccgttgagcgccagacggcagaacggttgattaacgcctgtc                                       >  2:172545‑M1/1‑81 (MQ=255)
                                                                                         tgtgGCATCTGGCGGGCGATGTTATCAGGAGTCTGCGATGTCTGTAACCCGTAATCTCTGGCGCGTTGCTGATGCGCCGattgttccggctgactccgttgagcgccagacggcagaacggttgattaacgcctgtccg                                     >  1:1006297‑M1/1‑79 (MQ=255)
                                                                                         tgtgGCATCTGGCGGGCGATGTTATCAGGAGTCTGCGATGTCTGTAACCCGTAATCTCTGGCGCGTTGCTGATGCGCCGattgttccggctgactccgttgagcgccagacggcag                                                            >  2:1233724‑M1/1‑79 (MQ=255)
                                                                                         tgtgGCATCTGGCGGGCGATGTTATCAGGAGTCTGCGATGTCTGTAACCCGTAATCTCTGGCGCGTTGCTGATGCGCCGattgttccggctgactccgttgagcgccagacggcag                                                            <  1:1233724‑M1/116‑38 (MQ=255)
                                                                                             gCATCTGGCGGGCGATGTTATCAGGAGTCTGCGATGTCTGTAACCCGTAATCTCTGGCGCGTTGCTGATGCGCCGattgttccggctgactccgttgagcgccagacggcagaacggttgattaacgcctgtccggcag                                 <  2:1006297‑M1/139‑65 (MQ=255)
                                                                                                tCTGGCGGGCGATGTTATCAGGAGTCTGCGATGTCTGTAACCCGTAATCTCTGGCGCGTTGCTGATGCGCCGattgttccggctgactccgttgagcgccagacggcagaa                                                          >  1:1079168‑M1/1‑72 (MQ=255)
                                                                                                tCTGGCGGGCGATGTTATCAGGAGTCTGCGATGTCTGTAACCCGTAATCTCTGGCGCGTTGCTGATGCGCCGattgttccggctgactccgttgagcgccagacggcagaa                                                          <  2:1079168‑M1/111‑40 (MQ=255)
                                                                                                   ggcgggcgATGTTATCAGGAGTCTGCGATGTCTGTAACCCGTAATCTCTGGCGCGTTGCTGATGCGCCGattgttccggctgactccgttgagcgccagacggcagaacggttgattaacgcctgtccggcaggtcttt                           >  1:976818‑M1/1‑69 (MQ=255)
                                                                                                   ggcgggcgATGTTATCAGGAGTCTGCGATGTCTGTAACCCGTAATCTCTGGCGCGTTGCTGATGCGCCGattgttccggctgactccgttgagcgccagacggcagaacgg                                                       <  1:940833‑M1/111‑43 (MQ=255)
                                                                                                   ggcgggcgATGTTATCAGGAGTCTGCGATGTCTGTAACCCGTAATCTCTGGCGCGTTGCTGATGCGCCGattgttccggctgactccgttgagcgccagacggcagaacgg                                                       >  2:940833‑M1/1‑69 (MQ=255)
                                                                                                     cgggcgATGTTATCAGGAGTCTGCGATGTCTGTAACCCGTAATCTCTGGCGCGTTGCTGATGCGCCGattgttccggctgactccgttgagcgccagacggcagaacggttgattaacgcctgtccggcaggtcttttt                         >  1:208711‑M1/1‑67 (MQ=255)
                                                                                                     cgggcgATGTTATCAGGAGTCTGCGATGTCTGTAACCCGTAATCTCTGGCGCGTTGCTGATGCGCCGattgttccggctgactccgttgagcgccagacggcagaacggttgattaacgcctgtccggcaggtcttttt                         >  2:192470‑M1/1‑67 (MQ=255)
                                                                                                     cgggcgATGTTATCAGGAGTCTGCGATGTCTGTAACCCGTAATCTCTGGCGCGTTGCTGATGCGCCGattgttccggctgactccgttgagcgccagacggcagaacggttgattaacgcctgtccggcaggtcttttt                         >  1:281784‑M1/1‑67 (MQ=255)
                                                                                                     cgggcgATGTTATCAGGAGTCTGCGATGTCTGTAACCCGTAATCTCTGGCGCGTTGCTGATGCGCCGattgttccggctgactccgttgagcgccagacggcagaacggttgattaacgcctgtccggcaggtctttt                          <  1:299890‑M1/138‑72 (MQ=255)
                                                                                                     cgggcgATGTTATCAGGAGTCTGCGATGTCTGTAACCCGTAATCTCTGGCGCGTTGCTGATGCGCCGattgttccggctgactccgttgagcgccagacggcagaacggttgattaacgcctgtccggcaggtctttt                          >  2:299890‑M1/1‑67 (MQ=255)
                                                                                                     cgggcgATGTTATCAGGAGTCTGCGATGTCTGTAACCCGTAATCTCTGGCGCGTTGCTGATGCGCCGattgttccggctgactccgttgagcgccagacggcagaacggttgattaacgcctgtccgg                                    <  1:420123‑M1/128‑62 (MQ=255)
                                                                                                     cgggcgATGTTATCAGGAGTCTGCGATGTCTGTAACCCGTAATCTCTGGCGCGTTGCTGATGCGCCGattgttccggctgactccgttgagcgccagacggcagaacggttgattaacgcctgtccgg                                    >  2:420123‑M1/1‑67 (MQ=255)
                                                                                                       ggcgATGTTATCAGGAGTCTGCGATGTCTGTAACCCGTAATCTCTGGCGCGTTGCTGATGCGCCGattgttccggctgactccgttgagcgccagacggcagaacggttgattaacgcctgtccggcaggtcttttttc                       >  2:1159721‑M1/1‑65 (MQ=255)
                                                                                                        gcgATGTTATCAGGAGTCTGCGATGTCTGTAACCCGTAATCTCTGGCGCGTTGCTGATGCGCCGattgttccggctgactccgttgagcgccagacggcagaacggttgattaacgcctgtccggcaggtcttttttcg                      <  1:172545‑M1/139‑76 (MQ=255)
                                                                                                         cgATGTTATCAGGAGTCTGCGATGTCTGTAACCCGTAATCTCTGGCGCGTTGCTGATGCGCCGattgttccggctgactccgttgagcgccagacggcagaacggttgattaacgcctgtccggcaggtcttttttcgc                     <  2:485914‑M1/139‑77 (MQ=255)
                                                                                                             gTTATCAGGAGTCTGCGATGTCTGTAGCCCGTAATCTCTGGCGCGTTGCTGATGCGCCGCACATTGTTCCGGCTGACTCCGTTGAGCGCCAGACGGCAGAACGGTTGATTAACGCCTGTCCGGCAGGTCTTTTTTCGCt                 <  2:729946/139‑1 (MQ=255)
                                                                                                              ttATCAGGAGTCTGCGATGTCTGTAACCCGTAATCTCTGGCGCGTTGCTGATGCGCCGattgttccggctgactccgttgagcgccagacggcagaacggttgattaacgcctgtccggcaggtcttttttcgctcaca                <  2:208711‑M1/139‑82 (MQ=255)
                                                                                                               tATCAGGAGTCTGCGATGTCTGTAACCCGTAATCTCTGGCGCGTTGCTGATGCGCCGattgttccggctgactccgttgagcgccagacggcagaacggttgattaacgcctgtccggcaggtcttttttcgctcacac               <  2:964816‑M1/139‑83 (MQ=255)
                                                                                                               tATCAGGAGTCTGCGATGTCTGTAACCCGTAATCTCTGGCGCGTTGCTGATGCGCCGattgttccggctgactccgttgagcgccagacggcagaacggttgattaacgcctgtccggcaggtcttttttcgctcacac               <  1:171101‑M1/139‑83 (MQ=255)
                                                                                                                  cAGGAGTCTGCGATGTCTGTAACCCGTAATCTCTGGCGCGTTGCTGATGCGCCGattgttccggctgactccgttgagcgccagacggcagaacggttgattaacgcctgtccggcaggtcttttttcgctcacaccgg            <  1:132350‑M1/139‑86 (MQ=255)
                                                                                                                  cAGGAGTCTGCGATGTCTGTAACCCGTAATCTCTGGCGCGTTGCTGATGCGCCGattgttccggctgactccgttgagcgccagacggcagaacggttgattaacgcctgtccggcaggtcttttttcgctcacaccgg            <  2:976818‑M1/139‑86 (MQ=255)
                                                                                                                       gTCTGCGATGTCTGTAACCCGTAATCTCTGGCGCGTTGCTGATGCGCCGattgttccggctgactccgttgagcgccagacggcagaa                                                          <  1:15103‑M1/88‑40 (MQ=255)
                                                                                                                       gTCTGCGATGTCTGTAACCCGTAATCTCTGGCGCGTTGCTGATGCGCCGattgttccggctgactccgttgagcgccagacggcagaa                                                          >  2:15103‑M1/1‑49 (MQ=255)
                                                                                                                        tCTGCGATGTCTGTAACCCGTAATCTCTGGCGCGTTGCTGATGCGCCGattgttccggctgactccgttgagcgccagacggcagaacggttgattaacgcctgtccggcaggtcttttttcgctcacaccggaaggta      >  1:555839‑M1/1‑48 (MQ=255)
                                                                                                                        tCTGCGATGTCTGTAACCCGTAATCTCTGGCGCGTTGCTGATGCGCCGattgttccggctgactccgttgagcgccagacggcagaacggttgattaacgcctgtccggcaggtcttttttcgctcacaccggaaggta      >  2:647797‑M1/1‑48 (MQ=255)
                                                                                                                        tCTGCGATGTCTGTAACCCGTAATCTCTGGCGCGTTGCTGATGCGCCGattgttccggctgactccgttgagcgccagacggcagaacggttgattaacgcctgtccggcaggtcttttttcgctcacaccgg            >  2:395434‑M1/1‑48 (MQ=255)
                                                                                                                        tCTGCGATGTCTGTAACCCGTAATCTCTGGCGCGTTGCTGATGCGCCGattgttccggctgactccgttgagcgccagacggcagaacggttgattaacgcctgtccggcaggtcttttttcgctcacaccgg            <  1:395434‑M1/133‑86 (MQ=255)
                                                                                                                          tGCGATGTCTGTAACCCGTAATCTCTGGCGCGTTGCTGATGCGCCGattgttccggctgactccgttgagcgccagacggcagaacggttgattaacgcctgtccggcaggtcttttttcgctcacaccggaaggtaac    >  2:1141284‑M1/1‑46 (MQ=255)
                                                                                                                            cGATGTCTGTAACCCGTAATCTCTGGCGCGTTGCTGATGCGCCGattgttccggctgactccgttgagcgccagacggcagaacggttgattaacgcctgtccggcaggtcttttttcgctcacaccggaaggtaactt  <  1:647797‑M1/139‑96 (MQ=255)
                                                                                                                            cGATGTCTGTAACCCGTAATCTCTGGCGCGTTGCTGATGCGCCGattgttccggctgactccgttgagcgccagacggcagaacggttgattaacgcctgtccggcaggtcttttttcgctcacaccggaaggtaactt  <  2:281784‑M1/139‑96 (MQ=255)
                                                                                                                                       |                                                                                                                               
GGATATTTCCCGCGATTTATGGGATCAGGGTGATAAACCTGTTCCACCGCTGCGCCAGTTATTCTGGCATCATTTACGTCGTCACGGCCTGTGGCATCTGGCGGGCGATGTTATCAGGAGTCTGCGATGTCTGTAGCCCGTAATCTCTGGCGCGTTGCTGATGCGCCGCACATTGTTCCGGCTGACTCCGTTGAGCGCCAGACGGCAGAACGGTTGATTAACGCCTGTCCGGCAGGTCTTTTTTCGCTCACACCGGAAGGTAA  >  NC_000913/2893793‑2894055

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