Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,942,523 T→C T21A (ACA→GCA)  znuA ← zinc ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,942,5230TC92.4% 198.9 / 4.3 66T21A (ACA→GCA) znuAzinc ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base T (3/2);  new base C (28/33);  total (31/35)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-01

GTAAAACCTCTGTTTCTGTTACCCCATCAGCAATGGCAGAAGCGATGAACCCAACGGGTTTAAGCGAAGCGACAACGGCGGCATCTGCGGCCTGTGTTGCACCTCCCCAGAGAGCGGCGGATAATGCTGCGAAAAGAAGCGTTTTTTTATGTAACATAATGCGACCAATAATCGTAATGAATATGAGAAGTGT  >  NC_000913/1942427‑1942619
                                                                                                |                                                                                                
gTAAAACCTCTGTTTCTGTTACCCCATCAGCAATGGCAGAAGCGATGAACCCAACGGGTTTAAGCGAAGCGACAACGGCGGCATCTGCGGCCTGtgctgc                                                                                               <  2:720546/100‑1 (MQ=255)
gTAAAACCTCCGTTTCTGTTACCCCATCAGCAATGGCAGAAGCGATGAACCCAACGGGATTAAGCGAAGCGACAACGGCAGCATCTGCGGCCTGTGTTGc                                                                                               >  2:81959/1‑100 (MQ=255)
  aaaaCCTCTGTTTCTGTTACCCCATCAGCAATGGCAGAAGCGATGAACCCAACGGGTTTAAGCGAAGCGACAACGGCGGCATCTGCGGCCTGtgctgc                                                                                               >  1:86295/1‑98 (MQ=255)
  aaaaCCTCTGTTTCTGTTACCCCATCAGCAATGGCAGAAGCGATGAACCCAACGGGTTTAAGCGAAGCGACAACGGCGGCATCTGCGGCCTGTGCTGCAc                                                                                             >  1:1081029/1‑100 (MQ=255)
       ctctGTTTCTGTTACCCCATCAGCAATGGCAGAAGCGATGAACCCAACGGGTTTAAGCGAAGCGACAACGGCGGCATCTGCGGCCTGTGCTGCACCTccc                                                                                        >  2:1072419/1‑100 (MQ=255)
         ctgtttctgttACCCCATCAGCAATGGCAGAAGCGATGAACCCAACGGGTTTAAGCGAAGCGACAACGGCGGCATCTGCGGCCTGTGCTGCACCTCCCCa                                                                                      >  2:1265642/1‑100 (MQ=255)
          tgtttctgttACCCCATCAGCAATGGCAGAAGCGATGAACCCAACGGGTTTAAGCGAAGCGACAACGGCGGCATCTGCGGCCTGTGCTGCACCTCCCCa                                                                                      <  2:384278/99‑1 (MQ=255)
           gtttctgttACCCCATCAGCAATGGCAGAAGCGATGAACCCAACGGGTTTAAGCGAAGCGACAACGGCGGCATCTGCGGCCTGTGCTGCACCTCCCCaga                                                                                    >  2:1384117/1‑100 (MQ=255)
           gtttctgttACCCCATCAGCAATGGCAGAAGCGATGAACCCAACGGGTTTAAGCGAAGCGACAACGGCGGCATCTGCGGCCTGTGCTGCACCTCCCCaga                                                                                    >  2:1567037/1‑100 (MQ=255)
             ttctgttACCCCATCAGCAATGGCAGAAGCGATGAACCCAACGGGTTTAAGCGAAGCGACAACGGCGGCATCTGCGGCCTGTGCTGCACCTCCCCagaga                                                                                  <  2:1405729/100‑1 (MQ=255)
                   tACCCCATCAGCAATGGCAGAAGCGATGAACCCAACGGGTTTAAGCGAAGCGACAACGGCGGCATCTGCGGCCTGTGCTGCACCTCCCCAGAGAgcggcg                                                                            <  1:1116900/100‑1 (MQ=255)
                   tACCCCATCAGCAATGGCAGAAGCGATGAACCCAACGGGTTTAAGCGAAGCGACAACGGCGGCATCTGCGGCCTGTGCTGCACCTCCCCAGAGAgcggcg                                                                            <  1:846106/100‑1 (MQ=255)
                   tACCCCATCAGCAATGGCAGAAGCGATGAACCCAACGGGTTTAAGCGAAGCGACAACGGCGGCATCTGCGGCCTGTGCTGCACCTCCCCAGAGAgcggcg                                                                            <  1:1623403/100‑1 (MQ=255)
                     ccccATCAGCAATGGCAGAAGCGATGAACCCAACGGGTTTAATCGAAGCGACAACGGCGGCATCTGCGGCCTGTGCTGCACCTCCCCAGAGAGCGGCGGa                                                                          <  1:361413/100‑1 (MQ=255)
                        cATCAGCAATGGCAGAAGCGATGAACCCAACGGGTTTAAGCGAAGCGACAACGGCGGCATCTGCGGCCTGTGCTGCACCTCCCCAGAGAGCGGCGGATaa                                                                       >  2:295524/1‑100 (MQ=255)
                         aTCAGCAATGGCAGAAGCGATGAACCCAACGGGTTTAAGCGAAGCGACAACGGCGGCATCTGCGGCCTGTGCTGCACCTCCCCAGAGAGCGGCGGATAAt                                                                      <  2:1150240/100‑1 (MQ=255)
                          tCAGCAATGGCAGAAGCGATGAACCCAACGGGTTTAAGCGAAGCGACAACGGCGGCATCTGCGGCCTGTGCTGCACCTCCCCAGAGAGCGGCGGATAAtg                                                                     >  2:105276/1‑100 (MQ=255)
                            aGCAATGGCAGAAGCGATGAACCCAACGGGTTTAAGCGAAGCGACAACGGCGGCATCTGCGGCCTGTGCTGCACCTCCCCAGAGAGCGGCGGATAAtgct                                                                   <  1:564436/100‑1 (MQ=255)
                              cAATGGCAGAAGCGATGAACCCAACGGGTTTAAGCGAAGCGACAACGGCGGCATCTGCGGCCTGTGCTGCACCTCCCCAGAGAGCGGCGGATAAtgctgc                                                                 >  2:1760791/1‑100 (MQ=255)
                               aaTTGCAGAAGCGATGAACCCAACGGGTTTAAGCGAAGCGACAACGGCGGCATCTGCGGCCTGTGCTGCACCTCCCCAGAGAGCGGCGGATAATGCTGCg                                                                <  1:1604471/100‑1 (MQ=255)
                               aaTGGCAGAAGCGATGAACCCAACGGGTTTAAGCGAAGCGACAACGGCGGCATCTGCGGCCTGTGCTGCACCTCCCCAGAGAGCGGCGGATAATGCTGCg                                                                <  2:224245/100‑1 (MQ=255)
                                  ggCAGAAGCGATGAACCCAACGGGTTTAAGCGAAGCGACAACGGCGGCATCTGCGGCCTGTGCTGCACCTCCCCAGAGAGCGGCGGATAATGCTGCGaaa                                                             >  2:689338/1‑100 (MQ=255)
                                    cAGAAGCGATGAACCCAACGGGTTTAAGCGAAGCGACAACGGCGGCATCTGCGGCCAGTGCTGCACCTCCCCAGAGAGCGGCGGATAATGCTGCGAAaag                                                           >  1:1158866/1‑100 (MQ=255)
                                     aGAAGCGATGAACCCAACGGGTTTAAGCGAAGCGACAACGGCGGCATCTGCGGCCTGTGCTGCACCTCCCCAGAGAGCGGCGGATAATGCTGCGAAgaga                                                          <  2:1819659/100‑1 (MQ=255)
                                       aaGCGATGAACCCAACGGGTTTAAGCGAAGCGACAACGGCGGCATCTGCGGCCTGTGCTGCACCTCCCCAGAGAGCGGCGGATAATGCTGCGAAaagaag                                                        <  2:558459/100‑1 (MQ=255)
                                         gtgATGAACCCAGATGTTTTAAGTGAAGCGACAACGGCGGTATCTGCGGCCTGTGCTGAACCTCCCCAGGGAGCGGCGGATAATGCTGCGAAaagaag                                                        <  2:725885/96‑1 (MQ=255)
                                         gCGATGAACCCAACGGGTTTAAGCGAAGCGACAACGGCGGCATCTGCGGCCTGTGCTGCACCTCCCCAGAGAGCGGCGGATAATGCTGCGAAAAGAAGCg                                                      <  2:268978/100‑1 (MQ=255)
                                          cGATGAACCCAACGGGTTTAAGCGAAGCGACAACGGCGGCATCTGCAGCCTGTGTTGCACCTCCCCAGAGAGCAACAAATAATGCAGCGAAAAGAAGCg                                                      >  2:1536669/1‑99 (MQ=255)
                                           gATGAACCCAACGGGTTTAAGCGAAGCGACAACGGCGGCATCTGCGGCCTGTGCTGCACCTCCCCAGAGAGCGGCGGATAATGCTGCGAAAAGAAGCGt                                                     >  2:1542604/1‑99 (MQ=255)
                                            aTGAACCCAACGGGTTTAAGCGAAGCGACAACGGCGGCATCTGCGGCCTGTGCTGCACCTCCCCAGAGAGCGGCGGATAATGCTGCGAAAAGAAGCGttt                                                   >  1:1044417/1‑100 (MQ=255)
                                             tGAACCCAACGGGTTTAAGCGAAGCGACAACGGCGGCATCTGCGGCCTGTGCTGCACCTCCCCAGAGAGCGGCGGATAATGCTGCGAAAAGAAGCGtttt                                                  >  2:1620586/1‑100 (MQ=255)
                                               aaCCCAACGGGTTTAAGCGAAGCGACAACGGCGGCATCTGCGGCCTGTGCTGCACCTCCCCAGAGAGCGGCGGATAATGCTGCGAAAAGAAGCGtttttt                                                >  2:1148554/1‑100 (MQ=255)
                                                aCCCAACGGGTTTAAGCGAAGCGACAACGGCGGCATCTGCGGCCTGTGCTGCACCTCCCCAGATAGCGGCGGATAATGCTGCGAAAAGAAGCGttttttt                                               >  2:1732422/1‑100 (MQ=255)
                                                aCCCAACGGGTTTAAGCGAAGCGACAACGGCGGCATCTGCGGCCTGTGCTGCACCTCCCCAGAGAGCGGCGGATAATGCTGCGAAAAGAAGCGttttttt                                               >  2:690427/1‑100 (MQ=255)
                                                aCCCAACGGGTTTAAGCGAAGCGACAACGGCGGCATCTGCGGCCTGTGCTGCACCGCCCCAGAGAGCGGCGGATAATGCTGCGAAAAGAAGCGtttttt                                                >  2:989973/1‑99 (MQ=255)
                                                 cccAACGGGTTTAAGCGAAGCGACAACGGCGGCATCTGCGGCCTGTGCTGCCCCTCCCCAGAGAGCGGCGGATAATGCTGCGAAAAGAAGCGTTTTTTTa                                              >  2:1681852/1‑100 (MQ=255)
                                                 cccAACGGGTTTAAGCGAAGCGACAACGGCGGCATCTGCGGCCTGTGCTGCACCTCCCCAGAGAGCGGCGGATAATGCTGCGAAAAGAAGCGTTTTTTTa                                              <  1:1298411/100‑1 (MQ=255)
                                                   cAACGGGTTTAAGCGAAGCGACAACGGCGGCATCTGCGGCCTGTGCTGCACCTCCCCAGAGAGCGGCGGATAATGCTGCGAAAAGAAGCGTTTTTTTATg                                            <  2:1093095/100‑1 (MQ=255)
                                                   cAACGGGTTTAAGCGAAGCGACAACGGCGGCATCTGCGGCCTGTGCTGCACCTCCCCAGAGAGCGGCGGATAATGCTGCGAAAAGAAGCGTTTTTTTATg                                            <  1:295522/100‑1 (MQ=255)
                                                   cAACGGGTTTAAGCGAAGCGACAACGGCGGCATCTGCGGCCTGTGCTGAACCTCCCCAGAGAGCGGCGGATAATGCTGCGAAGAGAAGCGTTTTTTTATg                                            >  1:1566667/1‑100 (MQ=255)
                                                   cAACGGGTTTAAGCGAAACGACAACGGCGGCATCTGCGGCCTGTGCTGCACCTCCCCAGAGAGCGGCCGATAATGCTGCGAAAAGAAGCGTTTTTTTATg                                            >  1:581417/1‑100 (MQ=255)
                                                       gggTTTAAGCGAAGCGACAACGGCGGCATCTGCGGCCTGTGCTGCACCTCCCCAGAGAGCGGCGGATAATGCTGCGAAAAGAAGCGTTTTTTTATGTAAc                                        <  1:619476/100‑1 (MQ=255)
                                                        ggTTTAAGCGAAGCGACAACGACGGCATCTGCGGCCTGTGCTGCACCTCACCAGAGAGCGGCGGAAAATGCTGCGAAAAGAAGCGTTTTTTTATGTAAc                                        >  1:779280/1‑99 (MQ=255)
                                                         gTTTAAGCGAAGCGACAACGGCGGCATCTGCGGCCTGTGCTGCACCTCCCCAGAGAGCGGCGGATAATGCTGCGAAAAGAAGCGTTTTTTTATGTAACAt                                      >  2:893501/1‑100 (MQ=255)
                                                          tttAAGCGAAGCGACAACGGCGGCATCTGCGGCCTGTGCTGCACCTCCCCAGAGAGCGGCGGATAATGCTGCGAAAAGAAGCGTTTTTTTATGTAACATa                                     >  2:78407/1‑100 (MQ=255)
                                                             aagcgaagcgACAACGGCGGCATCTGCGGCCTGTGCTGCACCTCCCCAGAGAGCGGCGGATAATGCTGCGAAAAGAAGCGTTTTTTTATGTAACATAATg                                  <  1:486272/100‑1 (MQ=255)
                                                              agcgaagcgaCAACGGCGGCATCTGCGGCCTGTGCTGCACCTCCCCAGAGAGCGGCGGATAATGCTGCGAAAAGAAGCGTTTTTTTATGTAACATAATGc                                 <  1:1713380/100‑1 (MQ=255)
                                                               gcgaagcgaCAACGGCGACATCTGCGGCCTGTGCTGCACCTCCCCCGAGAGCGGCGGATAATGCTGCGAAAAGAAGCGTTTTTTTATGTAACATAATACg                                >  1:300804/1‑100 (MQ=255)
                                                                cgaagcgaCAACGGCGGCATCTGCGGCCTGTGCTGCACCTCCCCAGAGAGCGGCGGATAATGCTGCGAAAAGAAGCGTTTTTTTATGCAACATAATGCGa                               <  2:1191578/100‑1 (MQ=255)
                                                                 gaagcgaCAACGGCGGCATCTGCGGCCTGTGCTGCACCTCCCCAGAGAGCGGCGGATAATGCTGCGAAAAGAAGCGTTTTTTTATGTAACATAATGCGAc                              >  1:1339385/1‑100 (MQ=255)
                                                                  aagcgaCAACGGTGGCATCTGCGGCCTGTGCTGCACATCCCCAGAGAGCGGCGGATAATGCTGCGAAAAGAAGCGTTTTTTTATGTAACATAATGCGAcc                             <  2:365751/100‑1 (MQ=255)
                                                                   agcgaCAACGGCGGCATCTGCGGCCTGTGCTGCACCTCCCCAGAGAGCGGCGGATAATGCTGCGAAAAGAAGCGTTTTTTTATGTAACATAATGCGACCa                            >  2:1082303/1‑100 (MQ=255)
                                                                     cgaCAACGGCGGCATCTGCGGCCTGTGCTGCACCTCCCCAGAGAGCGGCGGATAATGCTGCGAAAAGAAGCGTTTTTTTATGTAACATAATGCGACCaat                          <  2:107652/100‑1 (MQ=255)
                                                                      gaCAACGGCGGCATCTGCGGCCTGTGCTGCACCTCCCCAGAGAGCGGCGGATAATGCTGCGAAAAGAAGCGTTTTTTTATGTAACATAATGCGACCaata                         >  2:137835/1‑100 (MQ=255)
                                                                          aCGGCGGCATCTGCGGCCTGTGCTGCACCTCCCCAGAGAGCGGCGGATAATGCTGCGAAAAGAAGCGTTTTTTTATGTAACATAATGCGACCAATAATCg                     <  2:926682/100‑1 (MQ=255)
                                                                            ggcggcATCTGCGGCCTGTGCTGCACCTCCCCAGAGAGCGGCGGATAATGCTGCGAAAAGAAGCGTTTTTTTATGTAACATAATGCGACCAATAATCGTa                   <  1:1265631/100‑1 (MQ=255)
                                                                             gcggcATCTGCGGCCTGTGCTGCATCTCCCCAGAGAGCGGCGGATAATGCTGCGAAAAGAAGCGTTTTTTTATGTAACATAATGCGACCAATAATCGTaa                  <  2:1093113/100‑1 (MQ=255)
                                                                               ggcATCTGCGGCCTGTGCTGCACCTCCCCAGAGAGTGGCGGATAATGCTGCGAAAAGATGCGTTTTTTTATGTAACATAATGCGACCAATAAGCGTaaaa                <  1:560933/100‑3 (MQ=255)
                                                                               ggcATCTGCGGCCTGTGCTGCACCTCCCCAGAGAGCGGCGGATAATGCTGCGAAAAGAAGCGTTTTTTTATGTAACATAATGCGACCAATAATCGTAATg                <  1:818342/100‑1 (MQ=255)
                                                                               ggcATCTGCGGCCTGTGCTGCACCTCCCCAGAGAGCGGCGGATAATGCTGCGAAAAGAAGCGTTTTTTTATGTAACATAATGCGACCAATAATCGTAATg                <  1:1107061/100‑1 (MQ=255)
                                                                                gcATCTGCGGCCTGTGTTGCACCGCCCAAGAGAGCGGCGGATAATGCTGCGAAAAGAAGCGTTTTTTTCTGTAACATAATGCGACCAATCATCGTAATGa               <  1:81959/100‑1 (MQ=255)
                                                                                  aTCTGCGGCCTGTGCTGCACCTCCCCAGAGAGCGGCGGATAATGCTGCGAAAAGAAGCGTTTTTTTATGTAACATAATGCGACCAATAATCGTAATGAat             <  2:1339395/100‑1 (MQ=255)
                                                                                     tGCGGCGTGTGTTGCATCTCCCCAGTTAGCGGCGGATAATGCTGCGAAAAGAAGCGTTTTTTTATGTAACATAATGCGACCAATAATCGTAATGAATATg          <  1:774374/100‑1 (MQ=255)
                                                                                      gCGGCCTGTGTTGCACCTCCCCAGAGAGCGGCGGATAATGCTGCGAAAAGAAGCGTTTTTTTATGTAACATAATGCGACCAATAATCGTAATGAATATga         >  2:854567/1‑100 (MQ=255)
                                                                                      gCGGCCTGTGCTGCACCTCCCCAGAGAGCGGCGGATAATGCTGCGAAAAGAAGCGTTTTTTTATGTAACATAATGCGACCAATTATCGTAATGAata            >  1:1585044/1‑97 (MQ=255)
                                                                                       cGGCCTGTGCTGCACCTCCCCAGAGAGCGGCGGATAATGCTGCGAAAAGAAGCGTTTTTTTATGTAACATAATGCGACCAATAATCGTAATGAATATgag        <  1:1344099/100‑1 (MQ=255)
                                                                                         gCCTGTGCTGCACCTCCCCAGAGAGCGGCGGATAATGCTGCGAAAAGAAGCGTTTTTTTATGTAACATAATGCGACCAATAATCGTAATGAATATGAGaa      <  2:1081039/100‑1 (MQ=255)
                                                                                            tgtgCTGCACCTCCCCAGAGAGCGGCGGATAATGCTGCGAAAAGAAGCGTTTTTTTATGTAACATAATGCGACCAATAATCGTAATGAATATGAGAAg     <  1:787728/98‑1 (MQ=255)
                                                                                             gtgcTGCACCTCCCCAGAGAGCGGCGGATAATGCTGCGAAAAGAAGCGTTTTTTTATGTAACATAATGCGACCAATAATCGTAATGAATATGAGAAgtgt  >  1:1362521/1‑100 (MQ=255)
                                                                                                |                                                                                                
GTAAAACCTCTGTTTCTGTTACCCCATCAGCAATGGCAGAAGCGATGAACCCAACGGGTTTAAGCGAAGCGACAACGGCGGCATCTGCGGCCTGTGTTGCACCTCCCCAGAGAGCGGCGGATAATGCTGCGAAAAGAAGCGTTTTTTTATGTAACATAATGCGACCAATAATCGTAATGAATATGAGAAGTGT  >  NC_000913/1942427‑1942619

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: