Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,942,993 G→A G111E (GGG→GAG)  znuC → zinc ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,942,9930GA97.8% 137.0 / ‑5.3 45G111E (GGG→GAG) znuCzinc ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (20/24);  total (21/24)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

TGCGCATCGGCTATGTACCGCAGAAGCTGTATCTCGACACCACGTTGCCACTGACCGTAAACCGTTTTTTACGCTTACGCCCTGGTACACATAAAGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGGGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGCGTGTACTATTAGCGCGAGCATTGTTAAATCGCCCGCAATTATTAGTGCTGGATGAACCCACTCAAGGCGTGG  >  NC_000913/1942861‑1943115
                                                                                                                                    |                                                                                                                          
tGCGCATCGGCTATGTACCGCAGAAGCTGTATCTCGACACCACGTTGCCACTGACCGTAAACCGTTTTTTACGCTTACGCCCTGGTACACATAAAGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCt                                                                                                                      <  2:344344/139‑1 (MQ=255)
  cgcATCGGCTATGTACCGCAGAAGCTGTATCTCGACACCACGTTGCCACTGACCGTAAACCGTTTTTTACGCTTACGCCCTGGTACACATAAAGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCAc                                                                                                                        <  2:141710/135‑1 (MQ=255)
  cgcATCGGCTATGTACCGCAGAAGCTGTATCTCGACACCACGTTGCCACTGACCGTAAACCGTTTTTTACGCTTACGCCCTGGTACACATAAAGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCAc                                                                                                                        >  1:141710/1‑135 (MQ=255)
  cgcATCGGCTATGTACCGCAGAAGCTGTATCTCGACACCACGTTGCCACTGACCGTAAACCGTTTTTTACGCTTACGCCCTGGTACACATAAAGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGa                                                                                                                    >  2:368403/1‑139 (MQ=255)
                    cAGAAGCTGTATCTCGACACCACGTTGCCACTGACCGTAAACCGTTTTTTACGCTTACGCCCTGGTACACATAAAGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGa                                                                                                        <  1:10029/133‑1 (MQ=255)
                    cAGAAGCTGTATCTCGACACCACGTTGCCACTGACCGTAAACCGTTTTTTACGCTTACGCCCTGGTACACATAAAGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGa                                                                                                        >  2:10029/1‑133 (MQ=255)
                                     caccacGTTGCCACTGACCGTAAACCGTTTTTTACGCTTACGCCCTGGTACACATAAAGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGaa                                                                                 >  1:865973/1‑139 (MQ=255)
                                         acGTTGCCACTGACCGTAAACCGTTTTTTACGCTTACGCCCTGGTACACATAAAGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGc                                                                             <  2:865973/139‑1 (MQ=255)
                                          cGTTGCCACTGACCGTAAACCGTTTTTTACGCTTACGCCCTGGTACACATAAAGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCAc                                                                                                           >  2:803661/1‑108 (MQ=255)
                                          cGTTGCCACTGACCGTAAACCGTTTTTTACGCTTACGCCCTGGTACACATAAAGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCAc                                                                                                           <  1:803661/108‑1 (MQ=255)
                                          cGTTGCCACTGACCGTAAACCGTTTTTTACGCTTACGCCCTGGTACACATAAAGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCa                                                                            >  2:780302/1‑139 (MQ=255)
                                          cGTTGCCACTGACCGTAAACCGTTTTTTACGCTTACGCCCTGGTACACATAAAGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCa                                                                            <  2:917337/139‑1 (MQ=255)
                                          cGTTGCCACTGACCGTAAACCGTTTTTTACGCTTACGCCCTGGTACACATAAAGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCa                                                                            >  2:252132/1‑139 (MQ=255)
                                               ccACTGACCGTAAACCGTTTTTTACGCTTACGCCCTGGTACACATAAAGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGCgtg                                                                       <  1:252132/139‑1 (MQ=255)
                                               ccACTGACCGTAAACCGTTTTTTACGCTTACGCCCTGGTACACATAAAGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGCgtg                                                                       <  1:780302/139‑1 (MQ=255)
                                                       cGTAAACCGTTTTTTACGCTTACGCCCTGGTACACATAAAGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGGGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGCGTGTACTATTa                                                               >  1:639757/1‑139 (MQ=255)
                                                              cGTTTTTTACGCTTACGCCCTGGTACACATAAAGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGCGTGTACTATTAGCGCGAg                                                        <  1:368403/139‑1 (MQ=255)
                                                                   tttACGCTTACGCCCTGGTACACATAAAGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGCGTGTACTATTAGCGCGAGCATTg                                                   >  1:249470/1‑139 (MQ=255)
                                                                               cccTGGTACACATAAAGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGCGTGTACTATTAGCGCGAGCATTGTTAAATCGCCCg                                       <  2:639757/139‑1 (MQ=255)
                                                                                  tGGTACACATAAAGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGCGTGTACTATTAGCGCGAGCATTGTTAAATCGCCCGCaa                                    >  1:704905/1‑139 (MQ=255)
                                                                                  tGGTACACATAAAGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGCGTGTACTATTAGCGCGAGCATTGTTAAATCGCCCGCaa                                    >  2:869157/1‑139 (MQ=255)
                                                                                          aTAAAGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGCGTGTACTATTAGCGCGAg                                                        <  2:915665/111‑1 (MQ=255)
                                                                                          aTAAAGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGCGTGTACTATTAGCGCGAg                                                        >  1:915665/1‑111 (MQ=255)
                                                                                          aTAAAGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGCGTGTACTATTAGCGCGAGCATTGTTAAATCGccc                                        <  1:138009/127‑1 (MQ=255)
                                                                                          aTAAAGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGCGTGTACTATTAGCGCGAGCATTGTTAAATCGccc                                        >  2:138009/1‑127 (MQ=255)
                                                                                            aaaGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGCGTGTACTATTAGCGCGAGCATTGTTAAATCGCCCGCAATTATTAGTGc                          <  2:249470/139‑1 (MQ=255)
                                                                                              agaagaTATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAg                                                                           >  2:180203/1‑88 (MQ=255)
                                                                                              agaagaTATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAg                                                                           <  1:180203/88‑1 (MQ=255)
                                                                                              agaagaTATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGCg                                                                         >  2:12895/1‑90 (MQ=255)
                                                                                              agaagaTATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGCg                                                                         <  1:12895/90‑1 (MQ=255)
                                                                                                   atatTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGCGTGTACTATTAGCGCGAGCATTGTTAAATCGCCCGCAATTATTAGTGCTGGATGa                   >  2:652935/1‑139 (MQ=255)
                                                                                                    tatTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGctctc                                                                                           >  2:224357/1‑66 (MQ=255)
                                                                                                    tatTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGctctc                                                                                           <  1:224357/66‑1 (MQ=255)
                                                                                                      ttttGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGCGTGTACTATTAGCGCGAGCATTGTTAAATCGCCCGCAATTATTAGTGCTGGATGAAcc                <  2:704905/139‑1 (MQ=255)
                                                                                                        ttGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGCGTGTACTATTAGCGCGAGCATTGTTAAATCGCCCGCAattat                                >  2:786625/1‑121 (MQ=255)
                                                                                                        ttGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGCGTGTACTATTAGCGCGAGCATTGTTAAATCGCCCGCAattat                                <  1:786625/121‑1 (MQ=255)
                                                                                                          gCCTGCACTGAAACGTGTCCAGGCCGAGCTCCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGCTTGTACTATTAGCGCGAGCATTGTTAAATCGCCCGCAATTATTAGTGCTGGATGAACCCACt            <  1:652935/139‑1 (MQ=255)
                                                                                                          gCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGCGTGTACTATTAGCGCGAGCATTGTTAAATCGCCCGCAATTATTAGTGCTGGATGAACCCACt            <  1:869157/139‑1 (MQ=255)
                                                                                                              gCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGCGTGTACTATTAGCGCGAGCATTGTTAAATCGCCCGCAATTATTAGTGCTGGATGAACCCACTCAAg        >  2:183173/1‑139 (MQ=255)
                                                                                                                    aaaCGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGCGTGTACTATTAGCGCGAGCATTGTTAAATCGCCCGCAATTATTAGTGCTGGATGAACCCACTCAAGGCGTgg  <  2:426624/139‑1 (MQ=255)
                                                                                                                    aaaCGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGCGTGTACTATTAGCGCGAGCATTGTTAAATCGCCCGCAATTATTAGTGCTGGATGAACCCACTCAAGGCGTgg  <  1:183173/139‑1 (MQ=255)
                                                                                                                         tgtCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAg                                                                           <  2:105332/61‑1 (MQ=255)
                                                                                                                         tgtCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAg                                                                           >  1:105332/1‑61 (MQ=255)
                                                                                                                            ccAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGc                                                                          <  2:262275/59‑1 (MQ=255)
                                                                                                                            ccAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGc                                                                          >  1:262275/1‑59 (MQ=255)
                                                                                                                                    |                                                                                                                          
TGCGCATCGGCTATGTACCGCAGAAGCTGTATCTCGACACCACGTTGCCACTGACCGTAAACCGTTTTTTACGCTTACGCCCTGGTACACATAAAGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGGGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGCGTGTACTATTAGCGCGAGCATTGTTAAATCGCCCGCAATTATTAGTGCTGGATGAACCCACTCAAGGCGTGG  >  NC_000913/1942861‑1943115

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: