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.9% 147.1 / NA 48G111E (GGG→GAG) znuCzinc ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (27/20);  minor base C (1/0);  total (28/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

CAGAAGCTGTATCTCGACACCACGTTGCCACTGACCGTAAACCGTTTTTTACGCTTACGCCCTGGTACACATAAAGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGGGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGCGTGTACTATTAGCGCGAGCATTGTTAAATCGCCCGCAATTATTAGTGCTGGATGAACCCACTCAAGGCGTGGATGTA  >  NC_000913/1942881‑1943120
                                                                                                                |                                                                                                                               
cAGAAGCTGTATCTCGACACCACGTTGCCACTGACCGTAAACCGTTTTTTACGCTTACGCCCTGGTACACATAAAGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCaaa                                                                                                       >  1:962968/1‑139 (MQ=255)
   aaGCTGTATCTCGACACCACGTTGCCACTGACCGTAAACCGTTTTTTACGCTTACGCCCTGGTACACATAAAGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGc                                                                                                    >  1:514714/1‑139 (MQ=255)
     gCTGTATCTCGACACCACGTTGCCACTGACCGTAAACCGTTTTTTACGCTTACGCCCTGGTACACATAAAGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGctc                                                                                                  >  1:154118/1‑139 (MQ=255)
        gTATCTCGACACCACGTTGCCACTGACCGTAAACCGTTTTTTACGCTTACGCCCTGGTACACATAAAGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGc                                                                                                                               >  2:435132/1‑107 (MQ=255)
        gTATCTCGACACCACGTTGCCACTGACCGTAAACCGTTTTTTACGCTTACGCCCTGGTACACATAAAGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGc                                                                                                                               <  1:435132/107‑1 (MQ=255)
               gACACCACGTTGCCACTGACCGTAAACCGTTTTTTACGCTTACGCCCTGGTACACATAAAGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCg                                                                                        >  1:699943/1‑139 (MQ=255)
                   ccacGTTGCCACTGACCGTAAACCGTTTTTTACGCTTACGCCCTGGTACACATAAAGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAAc                                                                                    >  2:565407/1‑139 (MQ=255)
                    cacGTTGCCACTGACCGTAAACCGTTTTTTACGCTTACGCCCTGGTACACATAAAGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACg                                                                                   >  2:943874/1‑139 (MQ=255)
                        ttGCCACTGACCGTAAACCGTTTTTTACGCTTACGCCCTGGTACACATAAAGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGc                                                                               <  1:943874/139‑1 (MQ=255)
                           ccACTGACCGTAAACCGTTTTTTACGCTTACGCCCTGGTACACATAAAGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGCgtg                                                                            <  2:154118/139‑1 (MQ=255)
                           ccAATGACCGTAAACCGTTTTTTACGCTTACGCCCTGGTACCCAAAAAGAAAATATTTTGCCTGCCCTGAAACGTTTCCAGGCCGAGCCCCTGATTAACGCCCCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGCgtg                                                                            <  2:962968/139‑1 (MQ=255)
                                     tAAACCGTTTTTTACGCTTACGCCCTGGTCCACATAAAGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTACCGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGc                                                                                  >  1:404445/1‑123 (MQ=255)
                                     tAAACCGTTTTTTACGCTTACGCCCTGGTACACATAAAGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGc                                                                                  <  2:404445/123‑1 (MQ=255)
                                     tAAACCGTTTTTTACGCTTACGCCCTGGTACACATAAAGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGCGTGTACTATTAgc                                                                  >  1:495009/1‑139 (MQ=255)
                                           gTTTTTTACGCTTACGCCCTGGTACACATAAAGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCa                                                                                 <  2:1054211/118‑1 (MQ=255)
                                           gTTTTTTACGCTTACGCCCTGGTACACATAAAGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCa                                                                                 >  1:1054211/1‑118 (MQ=255)
                                                      ttacgcCCTGGTACACATAAAGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGc                                                                               <  1:725719/109‑1 (MQ=255)
                                                      ttacgcCCTGGTACACATAAAGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGc                                                                               >  2:725719/1‑109 (MQ=255)
                                                      ttacgcCCTGGTACACATAAAGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGCGTGTACTATTAGCGCGAGCATTGTTAAATc                                                 <  2:514714/139‑1 (MQ=255)
                                                             cTGGTACACATAAAGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGCGTGTACTATTAGCGCGAGCATTGTTAAATCGCCCGCa                                          <  1:866739/139‑1 (MQ=255)
                                                                gTACACATAAAGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGCGTGTACTATTAGCGCGAGCATTGTTAAATCGCCCGCAatt                                       >  2:701334/1‑139 (MQ=255)
                                                                  acacATAAAGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGCGTGTACTATTAGCGCGAGCATTGTTAAATCGCCCGCAattat                                     <  2:495009/139‑1 (MQ=255)
                                                                  acacATAAAGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGCGTGTACTATTAGCGCGAGCATTGTTAAATCGCCCGCAattat                                     <  2:236178/139‑1 (MQ=255)
                                                                      aTAAAGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGCGTGTACTATTAGCGCGAGc                                                            >  1:303949/1‑112 (MQ=255)
                                                                      aTAAAGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGCGTGTACTATTAGCGCGAGc                                                            <  2:303949/112‑1 (MQ=255)
                                                                      aTAAAGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGCGTGTACTATTAGCGCGAGCATTGTTAAATCGCCCGCAATTATTAGt                                 >  2:878877/1‑139 (MQ=255)
                                                                          agaagaTATTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGCGTGTACTATTAGCGCGAGCATTGTTAAATCGCCCGCAATTATTAGTGCTg                             >  2:673585/1‑139 (MQ=255)
                                                                               atatTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCaa                                                                                                        <  1:1054039/59‑1 (MQ=255)
                                                                               atatTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCaa                                                                                                        >  2:1054039/1‑59 (MQ=255)
                                                                                tatTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGCGTGTACTATTAGCGCGAGCATTGTTAAATCGCCCGCAATTATTAGTGCTGGATGaa                       >  2:1048183/1‑139 (MQ=255)
                                                                                 atTTTGCCTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGCGTGTACTATTAGCGCGAGCATTGTTAAATCGCCCGCAATTATTAGTGCTGGATGAAc                      <  1:1048183/139‑1 (MQ=255)
                                                                                       ccTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGCGTGTACTATTAGCGCGAGCATTGTTAAATCGccc                                             <  1:477568/110‑1 (MQ=255)
                                                                                       ccTGCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGCGTGTACTATTAGCGCGAGCATTGTTAAATCGccc                                             >  2:477568/1‑110 (MQ=255)
                                                                                          gCACTGAAACGTGTCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGCGTGTACTATTAGCGCGAGCATTGTTAAATCGCCCGCAATTATTAGTGCTGGATGAACCCACTCAAg             >  1:823291/1‑139 (MQ=255)
                                                                                           cACTGAAACGTGTCCAGGTCGAGCACCTGATTAATGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGCgtg                                                                            >  2:618342/1‑75 (MQ=255)
                                                                                           cACTGAAACGTGTCCAGGCCGAGCACCTGATTAATGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGCgtg                                                                            <  1:618342/75‑1 (MQ=255)
                                                                                                     tgtCCAGGCCGCGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGCGTGTACTATTAGCGCGAGCATTGTTAAATCGCCCGCAATTATTAGTGCTGGATGAACCCACTCAAGGCGTGGATGTa  >  2:877877/1‑139 (MQ=255)
                                                                                                     tgtCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGCgtg                                                                            <  2:383773/65‑1 (MQ=255)
                                                                                                     tgtCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGCgtg                                                                            >  1:383773/1‑65 (MQ=255)
                                                                                                     tgtCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGCg                                                                              <  2:310032/63‑1 (MQ=255)
                                                                                                     tgtCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGCg                                                                              >  1:310032/1‑63 (MQ=255)
                                                                                                     tgtCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGCGTGTACTATTAGCGCGAGCATTGTTAAAt                                                  >  2:46045/1‑91 (MQ=255)
                                                                                                     tgtCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGCGTGTACTATTAGCGCGAGCATTGTTAAAt                                                  <  1:46045/91‑1 (MQ=255)
                                                                                                     tgtCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGCGTGTACTATTAGCGCGAGCATTGTTAAATCGCCCGCAATTATTAGTGCTGGGTGAACCAACTCAAGGCGTGGCTGTa  >  2:491162/1‑139 (MQ=255)
                                                                                                     tgtCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGCGTGTACTATTAGCGCGAGCATTGTTAAATCGCCCGCAATTATTAGTGCTGGATGAACCCACTCAAGGCGTGGATGTa  >  2:50298/1‑139 (MQ=255)
                                                                                                     tgtCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGCGTGTACTATTAGCGCGAGCATTGTTAAATCGCCCGCAATTATTAGTGCTGGATGAACCCACTCAAGGCGTGGATGTa  >  2:423655/1‑139 (MQ=255)
                                                                                                      gtCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAg                                                                                <  2:277677/60‑1 (MQ=255)
                                                                                                      gtCCAGGCCGAGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAg                                                                                >  1:277677/1‑60 (MQ=255)
                                                                                                                |                                                                                                                               
CAGAAGCTGTATCTCGACACCACGTTGCCACTGACCGTAAACCGTTTTTTACGCTTACGCCCTGGTACACATAAAGAAGATATTTTGCCTGCACTGAAACGTGTCCAGGCCGGGCACCTGATTAACGCACCGATGCAAAAGCTCTCTGGTGGCGAAACGCAGCGTGTACTATTAGCGCGAGCATTGTTAAATCGCCCGCAATTATTAGTGCTGGATGAACCCACTCAAGGCGTGGATGTA  >  NC_000913/1942881‑1943120

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